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ABSTRACT 


The  development  of  a  unified  picture  of  short-pulse  high-intensity  multiphoton  processes, 
embracing  atoms,  molecules,  and  solids,  appears  possible  through  the  study  of  clusters.  Of 
particular  significance  are  possible  intra-cluster  processes  that  can  influence  the  mechanism  of 
ionization  and  lead  to  the  production  of  inner-shell  vacancies.  Inner-shell  excitation  leading  to 
prompt  X-ray  emission  is  specifically  considered  and  the  treatment  leads  to  the  definition  of  a 
critical  cluster  size  n^,  representing  the  achievement  of  maximal  X-ray  emission  from  the 
ensemble.  These  results  suggest  the  possibility  of  designing  a  new  class  of  molecular  materials 
optimized  for  the  efficient  production  and  amplification  of  X-rays.  Recent  experiments, 
specifically  designed  to  evaluate  certain  essential  predictions  of  the  cluster  picture,  have 
furnished  considerable  supporting  evidence  for  the  main  features  of  the  proposed  model.  These 
studies,  which  involved  Kr  and  Xe,  (1)  established  the  crucial  role  of  cluster  formation  in  the 
generation  of  soft  X-ray  (—  80  -  150  eV)  emission,  (2)  verified  the  scaling  of  this  new 
phenomenon  in  to  the  kilovolt  (~  800  -  2100  eV)  region,  (3)  provided  important  insights 
concerning  the  dynamical  character  of  the  strong-field  coupling  to  the  clusters  and  (4)  combined 
the  cluster  excitation  with  a  recently  discovered  mode  of  channeled  propagation.  These  results 
strongly  indicate  that  the  predicted  amplification  can  be  achieved. 
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INTRODUCTION 


Our  recent  work  has  been  aimed  at  the  study  of  the  properties  of  X-ray  emission  from 
multiphoton  excited  molecules  and  clusters.  As  discussed  in  this  report,  the  basic  scaling 
relations  stemming  from  a  general  analysis  [1]  of  prompt  X-ray  emission  from  ensembles  of 
closely  grouped  atoms  indicate  that  clusters  of  heavy  atoms  and/or  suitably  designed  molecules 
containing  heavy  atoms  have  unusually  favorable  properties  for  the  production  and  amplification 
of  X-rays  in  the  kilovolt  range.  The  scaling  laws  found  indicate  a  very  strong  dependence  on 
certain  atomic  and  molecular  properties  which  favors  heavy  materials  in  a  potent  and  obvious 
way.  In  this  analysis,  explicit  contact  is  made  with  six  independent  experimental  observations 
involving  three  materials  (Ar,  Kr,  and  Xe),  all  either  unexplained  or  anomalous  findings 
concerning  the  generation  of  short  wavelength  (X  <  200  A)  radiation.  Consistent  agreement  of 
these  data  with  the  analysis  is  found  in  every  case.  It  is  considered  significant  that  this  synthesis 
can  be  achieved  with  a  treatment  that  only  involves  known  atomic  and  molecular  properties 
without  the  incorporation  of  additional  free  parameters. 

Four  new  experimental  results  [2,3]  concerning  multiphoton  induced  emissions  (1)  from 
Xe  in  the  145  -  152  A  range,  (2)  from  Kr  in  the  72  -  140  A  region,  (3)  from  Xe  in  the  11  -  17 
A  range,  and  (4)  from  Kr  in  the  7  A  region,  arising  from  experiments  explicitly  designed  to  test 
the  cluster  model,  provide  strong  confirming  evidence  for  the  physical  picture  presented.  These 
results  support  the  concept  of  a  molecular  design  of  specific  systems  for  the  production  and 
amplification  of  X-rays. 

The  prospects  of  X-ray  amplification  are  enormously  enhanced  by  knowledge  stemming 
from  previous  studies  [4,5]  of  confined  propagation  in  which  channel  formation  has  been 
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observed  [4],  This  physical  mechanism  enables  the  spatial  distribution  of  the  excited  material 
to  be  efficiently  produced  in  a  geometry  suitable  for  amplification  [6-8]. 

The  overall  goal  of  our  program  is  the  union  of  these  two  physical  processes,  namely, 
(1)  the  cluster  excitation  and  (2)  the  channeled  propagation,  in  order  to  produce  X-ray  sources 
of  exceptional  brightness  for  3-dimensional  microimaging, 
n.  DISCUSSION  OF  RESEARCH 

A.  Multiphoton  X-Ray  Generation  from  Molecules  and  Clusters 

A  unified  picture  of  short-pulse  high-intensity  multiphoton  processes,  embracing  atoms, 
molecules,  and  solids  has  been  proposed  in  relation  to  the  study  of  the  behavior  of  clusters  [1]. 
This  model  is  described  in  Appendix  A.  Specific  experimental  tests  of  the  cluster  concept  have 
also  been  performed,  the  most  informative  of  which  were  conducted  in  the  kilovolt  region  and 
involved  the  observation  of  Kr  L-shell  and  Xe  M-shell  emissions.  These  findings  are  presented 
in  Appendices  B  and  C. 

Important  scaling  relationships  have  been  established  by  these  studies,  particularly 
concerning  the  energy  transfer  rate.  These  results  are  given  in  Appendix  D.  The  main 
conclusion  is  that  energy  deposition  rates  exceeding  ~  1  W/atom  are  feasible  in  appropriately 
designed  clusters  or  molecules  incorporating  heavy  atoms. 

B.  Channeled  Propagation 

The  spatial  control  of  the  energy  deposited  for  excitation  plays  a  crucial  role  in 
establishing  the  conditions  for  X-ray  amplification.  The  existence  of  a  confined  mode  of 
propagation,  which  occurs  under  conditions  identical  with  those  suitable  for  excitation  of  the 
clusters,  is  now  supported  by  both  theoretical  [6,9]  analysis  and  experimental  [10]  evidence,  as 
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described  in  Appendices  E  and  F,  respectively.  This  phenomenon  has  an  enormous  influence 
on  the  scaling  relationships  [8]  governing  amplification,  the  analysis  of  which  is  presented  in 
Appendix  G.  The  main  conclusion  is  that  ~  1  J  of  excitation  energy  is  sufficient  to  produce 
strong  amplification  at  quantum  energies  well  into  the  kilovolt  range, 
m.  RESEARCH  FINDINGS 

The  goal  of  our  work  is  the  union  of  the  material  excitation,  discussed  in  Section  II.  A, 
with  the  confined  propagation,  described  in  Section  II.B.  This  step  would  bring  into  conjunction 
the  two  basic  requirements  for  efficient  X-ray  amplification,  namely,  (1)  a  highly  excited 
nonequilibrium  medium  and  (2)  a  spatial  distribution  configured  for  directional  amplification. 

Preliminary  experimental  findings  obtained  by  imaging  the  X-ray  emitting  zone  with  a 
pin-hole  camera  indicate  that  this  conjunction  is  possible.  The  results  involve  Xe  M-shell 
emission  (~  1  keV)  produced  under  the  conditions  described  in  Appendix  C.  The  pin-hole 
camera  was  positioned  so  that  a  transverse  view  of  the  excited  region  was  observed.  The 
camera  was  equipped  with  either  a  100  ^m  or  5  fim  pin-hole  and  a  filter  that  eliminated  radiation 
with  a  quantum  energy  less  than  ~  500  eV.  The  signal  was  recorded  on  a  two-dimension  array 
and  stored  in  a  microprocessor.  The  exposures  were  sufficiently  strong  that  single-pulse  images 
were  recorded. 

The  image  of  the  Xe  M-shell  kilovolt  radiation  recorded  with  the  100  ^m  pin-hole  is 
illustrated  in  Fig,  (1).  The  z-coordinate  represents  the  direction  of  propagation  (right  to  left) 
and  the  transverse  direction  is  designated  by  the  r-coordinate.  The  spatial  resolution  of  the 
image  is  ~  100  /im,  the  dimension  of  the  pin-hole.  A  line-out  along  the  axis  [(r,z)  =  (0,z)] 
is  shown  in  Fig.  (2).  A  clear  asymmetry  about  the  focal  position  of  the  248  nm  radiation  is 
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Direction  of  248  nm  Propagation 


Fig.  (1):  Pin-hole  X-ray  camera  image  in  transverse  direction  indicating  the  spatial 

distribution  of  Xe  M-shell  radiation  (-  1  keV).  The  pin-hole  had  a 
diameter  of  —  100  fim.  The  stagnation  pressure  of  the  pulsed-valve  was 
30  psi  and  the  temperature  was  293  K. 
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Fig.  (2);  Line-out  [(r,z)  =  (0,z)]  of  the  data  shown  in  Fig.  (1).  The  position  of  the 
248  nm  focal  point  (z  =  625  /xm)  is  indicated  along  with  a  marker 
indicating  the  Rayleigh  range  (—30  ^m)  corresponding  to  the  focusing 
system  used  for  the  248  nm  radiation.  The  length  5  represents  estimated 
length  of  the  region  in  which  material  from  the  pulse-valve  is  present. 
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present.  Significantly,  the  brighter  zone  is  on  the  fax  side  of  the  focus  and  extends  for  a  distance 
(200  or  300  a  value  many  times  the  Rayleigh  range  (~  30  ^m)  of  the  optical  system  used 
for  focusing.  Figure  (3)  illustrates  the  image  obtained  with  a  5  /im  pin-hole  while  Fig.  (4) 
shows  the  corresponding  axial  line-out.  The  asymmetry  is  clearly  present  in  the  latter  as  in  Fig. 
(2).  Careful  inspection  of  the  data  from  other  exposures  reveals  the  existence  of  a  narrow 
filament  with  a  length  of  ~  200  /nm  and  a  radial  extent  limited  by  the  resolution  of  the  imaging 
system  (~  20  /urn).  This  radial  profile  is  shown  in  Fig.  (5),  which  illustrates  a  transverse  line- 
out  of  the  filamentary  region. 

The  images  show  that  the  power  density  of  the  incident  248  nm  radiation  is  confined  in 
the  region  beyond  the  focus.  This  type  of  behavior  would  be  expected  if  the  channeling 
mechanism  discussed  in  Appendices  E  and  F  were  influencing  the  propagation.  All  the 
experimental  data  indicate  that  this  is  the  case. 

IV.  CONCLUSIONS 

The  union  of  a  (1)  multiphoton  cluster/ molecule  mechanism  for  the  production  of  X-rays 
and  (2)  a  relativistic/charge-displacement  process  for  the  establishment  of  confined  propagation 
appears  possible.  This  combination  intrinsically  has  exceptionally  favorable  characteristics  for 
the  efficient  amplification  of  X-rays  in  the  kilovolt  spectral  region. 
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Direction  of  248  nm  Radiation 


Fig.  (3):  Pin-hole  X-ray  camera  image  in  transverse  direction  indicating  the  spatial 

distribution  of  Xe  M-shell  radiation  (—  1  keV).  The  pin-hole  had  a 
diameter  of  5  ^m.  The  stagnation  pressure  of  the  pulsed-valve  was  30  psi 
and  the  temperature  was  273  K. 


Intensity  (orb.  units') 


Rayleigh  Range 


-  30  um  c  < _  Direction  of 

•2  248  nm  Radiation 


Fig.  (4):  Line-out  [(r,z)  =  (0,z)]  of  the  data  shown  in  Fig.  (3).  The  approximate 

position  of  the  248  nm  focal  point  (z  =  625  inn)  is  indicated  along  with 
a  marker  indicating  the  Rayleigh  range  ( ~  30  /^m)  corresponding  to  the 
focusing  system  used  for  the  248  nm  radiation. 
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Fig-.  (5);  Radial  line-out  of  observed  narrow  filamentary  region  for  an  Xe  pressure 
of  30  psi  and  a  temperature  of  273  K.  The  measured  FWHM  of  20  fun 
corresponds  to  the  limiting  spatial  resolution  of  the  X-ray  imaging  system 
used. 
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Abstract.  The  development  of  a  unified  picture  of  short- 
pulse  high-intensity  multiphoton  processes,  embracing 
atoms,  molecules,  and  solids,  appears  possible  through  the 
study  of  clusters.  Of  particular  significance  axe  ptossible 
intra-cluster  processes  that  can  influence  the  mechanism  of 
ionization  and  lead  to  the  production  of  inner-shell  vacan¬ 
cies.  Inner-shell  excitation  leading  to  prompt  X-ray  emission 
is  specifically  considered  and  the  treatment  leads  to  the  def¬ 
inition  of  a  critical  cluster  size  representing  the  achieve¬ 
ment  of  maximal  X-ray  emission  from  the  ensemble.  These 
results  suggest  the  possibility  of  designing  a  new  class  of 
molecular  materials  optimized  for  the  efficient  production 
«  and  amplification  of  X-rays. 

PACS:  97.85.  33.(K),  36.00,  42.55.Vc 


This  article  introduces  and  explores  the  possibility  of  de¬ 
veloping  a  unified  picture  of  short-pulse  « 100  fs)  high- 
intensity  (>  10‘^  W/cm*)  multiphoton  processes  embracing 
atoms,  molecules,  and  solids,  through  the  study  of  clusters 
[1],  Rare-gas  clusters  on  account  of  their  simplic¬ 

ity  and  close  connection  with  known  atomic  and  molecu¬ 
lar  properties,  appear  particularly  suitable  for  this  purpose. 
Of  particular  significance  are  possible  intra-cluster  processes 
that  can  influence  the  mechanism  of  ionization.  This  can  (1) 
lead  to  the  production  of  inner-shell  vacancies  and  (2)  the 
generation  of  prompt  X-rays  with  a  potential  for  strong  am¬ 
plification.  The  analysis  developed  lilow  gives  simple  rela¬ 
tions  concerning  inner-shell  excitation  and  X-ray  emission 
applicable  over  a  large  range  of  cluster  sizes  (1  <  n  <  10^). 
The  treatment  also  leads  to  the  definition  of  a  critical  clus¬ 
ter  size  representing  the  achievement  of  maximal  X-ray 
emission  from  the  ensemble,  a  parameter  that  is  found  to  be 
directly  analogous  to  the  concept  of  nuclear  critical  mass. 

Explicit  contact  is  made  with  six  independent  experimen¬ 
tal  observations  [2-4]  involving  three  materials  (Ar,  Kr,  Xe), 
all  either  unexplained  or  anomalous  findings  concerning  the 
generation  of  short  wavelength  (A  <  2(X)  A)  radiation.  Con¬ 


sistent  agreement  with  the  analysis  is  found  in  each  case. 
Significantly,  this  synthesis  is  achieved  with  the  exclusive 
use  of  previously  established  atomic  and  molecular  prop¬ 
erties.  A  connection  with  studies  of  X-ray  generation  from 
solid  surfaces  (BaFj)  is  also  made  in  conjunction  with  this 
analysis. 


1  Coupling  to  Clusters 
l.l  Simple  Initial  Formulation 

The  rare-gas  cluster  is  considered  as  a  spherical  group  of 
(n)  atoms  having  mutual  interatomic  spacings  Tq  and  an 
overall  radius  R  for  the  ensemble  given  by  R  = 

>  3),  and  R  (n  =  1,  2),  in  analogy  with 

simple  descriptions  [5, 6]  of  nuclei.  Consider  irradiation  of 
the  cluster  at  wavelength  A  with  the  stipulation  (a)  that 
the  characteristic  dimension  i?  of  a  fully  ionized  cluster 
be  less  than  the  skin  depth  6  and  (b)  the  assumption  that 
the  atoms,  which  are  weakly  bound  by  van  der  Waals 
forces,  respond  to  the  external  optical  field  as  free  atoms 
whose  ionization  can  be  estimated  with  a  tunneling  picture. 
Condition  (a)  enables  all  atoms  in  the  cluster  simultaneously 
to  experience  the  same  field  strength,  since  R  <  6  X. 
Consistent  with  these  assumptions,  we  will  set  the  maximum 
dimension  R^  ~  20-30  A,  or  equivalently,  a  maximum 
cluster  size  of  ~  lO*.  Assumption  (b)  enables  a  simple 
experimentally  evaluated  procedure  to  be  used  to  estimate 
the  level  of  atomic  ionization  produced.  This  gives  the 
threshold  intensity  for  ionization  [7, 8]  as 

4  =  c£^/1287re«Z^  (1) 

where  is  the  ionization  potential,  Z  the  resulting  ionic 
charge,  c  the  speed  of  light,  and  e  the  electronic  charge. 

1.2  Intra-Cluster  Processes 

The  model  of  a  rare-gas  cluster  as  a  sum  of  independent 
atoms,  held  weakly  together  by  van  der  Waals  forces  and 
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responding  as  individual  systems  to  the  external  optical 
field,  is  expected  to  become  invalid  as  the  cluster  becomes 
sufficiently  large  and/or  the  external  electromagnetic  field 
becomes  sufficiently  strong.  This  limit  could  arise  through 
the  action  of  inelastic  scattering  processes  produced  by 
ionized  electrons  originating  from  atoms  in  the  cluster  that 
are  accelerated  by  the  external  optical  field  and  collide  with 
other  atoms  or  ions  in  the  system.  Roughly,  such  intra-cluster 
processes  can  be  considered  as 


e-  +  X’+^e-+(X-'+)*,  (2) 

which  produces  excitation,  or  the  reaction 

(3) 

which  causes  additional  ionization. 

An  approximate  estimate  of  the  number  of  such 
excitations  or  ionizations  can  be  written  as  the  product  of  the 
number  of  electrons  produced  in  the  cluster  (nZ),  the  atomic 
density  in  the  cluster  the  cross  section  for  the  inelastic 
channel  under  consideration  (e.g.  and  the  characteristic 
scale  length  of  the  cluster  (iZ).  Therefore,  for  reaction  (3), 

N,  ^  ^  ,  (n  >  3) ,  (4) 

with  Z  determined  through  (I)  for  the  corresponding  inten¬ 
sity.  Such  a  mechanism  would  cause  the  pattern  of  ionization 
to  depend  upon  the  cluster  size  (n),  thereby  yielding  an  ob¬ 
served  spectrum  of  ionic  charge  states  deviating  from  that 
characteristic  of  the  free  atom  (71=  1),  an  easily  measured 
feature  in  direct  conflict  with  the  picture  of  atomic  indepen¬ 
dence  assumed  above  in  the  simplest  initial  formulation. 

An  important  immediate  consequence  of  (4)  is  the  exis¬ 
tence  of  a  critical  size 


for  inner-shell  excitation 

^2 


4  = 


STcan^/^  /hr^ 


(n  >  3) , 


(6) 


with  Q  =  e^jhc.  Since  this  simple  picture  assumes  that  the 
acceleration  occurs  in  a  time  comparable  to  or  less  than  a 
period  of  the  wave,  validity  of  this  approximation  demands 
that  the  cluster  neither  be  too  large  nor  the  external  field  too 
small;  namely,  it  requires  that  the  amplitude  of  the  driven 
excursion  of  the  electron  caused  by  the  field  be  comparable 
to  or  larger  than  the  cluster  dimension  R.  If  we  set  the 
limiting  situations  zs  =  R,  assume  that  the  electron  is 
free,  and  put  x^  =  eElmufl,  with  m  denoting  the  electron 
mass  and  u  the  angular  frequency  of  the  external  optical 
field,  we  can  define  another  limiting  intensity  as 


(7) 


with  A  denoting  the  wavelength  of  the  optical  field  and  \ 
being  the  electron  Compton  wavelength.  In  this  picture,  an 
inner-shell  electron  with  binding  energy  can  be  ionized  if 
the  intensity  exceeds  a  lower  bound  given  by  /g  or  /„(A), 
whichever  is  greatest 


2  Prompt  X-Ray  Emission 

Excitation  of  the  type  discussed  above  can  lead  to  the 
emission  of  prompt  inner  0  -1)  shell  radiation  in  the  X- 
ray  range.  Generally,  in  addition  to  the  production  of  a  hole 
in  the  0  -  1)  shell,  the  possibility  of  prompt  j  -*  j  -  I 
transitions  requires  retention  of  at  least  one  electron  in  the 
j-shell  during  the  course  of  irradiation.  This  determines  an 
upper  bound  on  the  intensity  of  irradiation  I^(j)  through 
(1)  which  is,  in  the  simple  model  considered  independent 
of  the  cluster  size. 


defined  by  =  n.  Therefore,  if  corresponds  to  an  itmer- 
shell  ionization,  with  n  >  an  inner-shell  vacancy  would 
be  produced  in  every  atom.  Furthermore,  since  the  critical 
size  of  the  cluster  depends  inversely  on  the  square  of  the 
atomic  density  {q~'^  =  r§)  and  inversely  on  the  cube  of  the 
cross  section  for  excitation  is  a  direct  analogue  of 

the  nuclear  concept  of  a  criticaii  mass,  with  cr^  replacing  the 
neutron  fission  cross  section  cr^  and  Z  playing  the  role  of 
the  average  neutron  yield  P. 

Certain  conditions  must  be  met  in  order  for  intra-cluster 
processes  of  the  type  expressed  by  reactions  (2)  and  (3) 
to  occur  for  an  inner-shell.  Specifically,  such  an  excitation 
can  only  take  place  if  sufficient  kinetic  energy  is  possessed 
by  the  colliding  electron.  To  estimate  this  energy,  and  the 
conditions  necessary  for  its  transfer,  we  will  assume  that  the 
electrons  experience  an  acceleration  over  the  characteristic 
dimension  of  (iZ)  of  the  cluster  by  the  external  optical  field 
(E)  in  a  time  short  compared  to  the  period  of  the  oscillating 
wave.  Furthermore,  once  the  electrons  leave  the  vicinity  of 
the  cluster  we  will  assume  that  they  represent  ionization 
and  are  lost  from  the  system.  This  simple  procedure  yields  a 
characteristic  energy  given  by  Si  eER,  a  result  that  can 
be  expressed  as  a  characteristic  minimum  intensity  needed 


2.1  Krypton 


Considering  Kr  Af-shell  emission  as  an  example,  from 
the  threshold  ionization  curve  [7-9]  derived  from  (1)  and 
exhibited  in  Fig.  1,  we  find  that  4ux(iV)  S  2  x  10'*  W/cm^. 
By  combining  the  conditions  stemming  from  (6)  and  (7) 
4uu(^)’  ^8-  2  illustrates  the  conditions  on  intensity 
and  cluster  size  for  which  the  model  allows  prompt  M-shell 
emissions.  The  region  corresponding  to  prompt  L-shell  [10] 
emissions  =  1.67  keV)  is  also  represented  in  Hg.  2. 

Several  important  results  are  manifest  in  the  form  of 
Fig.  2.  First,  the  crossing  points  n^{P)  and  ti,(Q),  found 
by  equating  (1)  with  (6)  and  (7),  respectively,  yield  the 
condition  for  the  existence  of  an  allowed  zone  for  prompt 
emission  as 


(8) 


an  expression  which  we  observe  is  independent  of  Tq. 
Furthermore,  (8)  strongly  favors  the  situation  in  which  the 
magnitudes  of  the  characteristic  energies,  and  hence,  the 
corresponding  levels  of  ionization  Z,  are  high.  This  is 
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Classical  MPI  of  KR 


Fig.  1.  Threshold  intensity  /,^  for  ionization  shown  as  a  function  of 
resulting  charge  state  Z  calculated  for  Kr  with  the  classical  method 
described  in  [7,8].  The  respective  atomic  shells  (K.  L,  M,  N)  are 
designated.  The  maximum  intensities  and  I^(N)  associated 

with  leaving  some  lU-sheil  and  iV-shell  populations,  respectively,  are 
indicated.  Also  shown  is  the  intensity  (~5x  10'*  W/cm'^)  at  which  Kr’"'' 
transitions  at  ~  100  A  were  observed  in  fluorescence  in  experiments 
reported  in  [4, 1 1] 


Cluster  Size  (n) 


□  Kr  K-^ll  (-100  @V)  ©LAM.  Expgrlmant 
^  Kr  L— ^lall  (l,S-s*2  ksV}  Zons 

Fig.  2.  Representation  of  limiting  conditions  applicable  for  production 
of  prompt  Kr  L-sheil  and  Af-shell  emission  from  Kr  clusters  as  a 
function  of  intensity  of  irradiation  and  cluster  size.  The  experimental 
datum  (LAND  corresponding  to  Lee  et  ai.  [2]  and  the  crossing  points 
P  and  Q  are  shown.  The  additional  condition  /„(248  nm),  stemming 
from  (7),  is  also  represented.  See  text  for  other  details 

readily  seen  by  letting  be  approximated  by  hydrogenic 
binding  {E^o:Z^),  keeping  E’p/eg  ~  constant,  and  noting 
that  Hq/rtp  ~  follows.  As  a  natural  consequence  of  the 
systematics  of  atomic  binding  energies,  we  observe  that  large 


Fig.  3a.  b.  Spectra  of  Kr  in  the  region  near  100  A.  a  Multiphoton 
produced  spectrum  of  Kr  produced  in  a  gas  jet,  as  described  in 
(11),  in  the  —85-100  A  region  with  a  peak  intensity  of  irradiation 
of  —5  X  10'*W/cm^.  The  4p  -►  3d  features  from  Kr*^  appear 
prominently,  b  Spectrum  of  Kr  observed  in  a  theta  pinch  plasma  by 
Jones  and  Kailne  [12],  The  4p  —►  3d  transitions  are  indicated.  The 
correspondence  between  the  multiphoton  and  the  theta  pinch  produced 
spectra  involving  Kr’"'’  in  the  — 100  A  region  is  marked 


allowed  zones  are  associated  with  level  stmctures  in  the 
kilovolt  range,  a  feature  apparent  from  comparison  of  the 
regions  of  M-shell  and  L-shell  emission  shown  in  Fig.  2. 

The  information  in  Fig.  2  can  be  directly  related  to  two  in¬ 
dependent  experimental  studies,  one  [2]  involving  the  search 
for  prompt  multiphoton-induced  X-ray  emission  from  free 
Kr  atoms  (n  =  1)  and  another  concerning  the  measurement 
of  fluorescence  produced  in  a  pulsed  g^  target  undergoing 
multiphoton  excitation  [4, 11],  The  datum  corresponding  to 
the  former  [2]  (/  S  3  x  10*^  W/cm^,  n  =  1),  designated 
as  LANL,  appears  on  lug.  2.  Although  the  experiment  ac¬ 
tually  represents  all  intensities  I  <  10*’  W/cm’  due  to  the 
distribution  of  intensity  occurring  in  the  focal  region,  we  see 
from  Fig.  2  that  there  is  no  overlap  of  the  experimental  con¬ 
ditions  with  either  of  the  zones  for  prompt  L-  or  M-shell 
emissions.  The  experimental  region  lies  well  below  the  L- 
shell  region  and  is  excluded  from  the  M-shell  zone  by  the 
crossing  P.  Hence,  prompt  X-ray  emission  is  forbidden,  a 
prediction  fully  in  agreement  with  the  ex|«rimental  obser¬ 
vation  of  none  [2]. 

The  second  experiment  involves  the  anomalous  detec¬ 
tion  [4, 11]  of  intense  Kr  M-shell  fluorescence  from  Kr^ 
4p  — ►  3d  transitions  at  ~  100  A,  as  illustrated  in  Fig.  3a, 
arising  from  material  produced  in  a  pulsed  gas  jet  and  irra- 
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dialed  at  a  peak  intensity  of  ~5  x  10'*  W/cm-  with  subpi¬ 
cosecond  248  nm  radiation.  The  corresponding  spectral  re¬ 
gion  observed  from  a  theta  pinch  by  Jones  and  Kallne  [12] 
is  shown  in  Fig.  3b. 

The  observation  of  any  Kr®'*"  radiation  by  the  multiphoton 
excitation,  either  prompt  or  delayed,  immediately  presents 
a  paradox.  As  shown  in  Fig.  1,  with  a  peak  intensity  of 
~5  X  10‘*  W/cm-,  no  production  of  Kr^'^  is  possible,  a 
result  that  is  experimentally  based  on  earlier  measurements 
of  ion  production  [8]  involving  free  krypton  atoms  under 
identical  conditions  of  irradiation.  Yet  the  Kr®'*'  array  at  ~ 
103  A  was  among  the  brightest  features  observed  throughout 
the  entire  region  studied  [4, 1 1].  However,  Fig.  2  shows  that 
the  formation  of  small  clusters  (n  >  2—3)  would  permit  both 
anomalous  ionization  (Kr®'*'  production)  and  the  emission 
of  Kr^'*'  4  — ►  3  radiation  of  the  kind  illustrated  in  Fig.  3, 
thereby  reconciling  the  paradox.  Furthermore,  a  review  of 
the  hydrodynamic  properries  and  operating  conditions  [4, 11] 
of  the  pulsed  valve,  particularly  in  comparison  to  other 
work  [13-15]  in  which  Ktj  and  Ktj  species  were  formed 
and  directly  measured  [14],  showed  that  the  hydrodynamic 
conditions  were  essentially  the  same  in  both  studies.  Indeed, 
the  observed  pressure  dependence  of  the  Kr®'*'  signal  was 
close  to  that  seen  in  the  earlier  studies  [14]  involving  Ktj 
and  Ktj  formation.  Both  of  these  considerations  lead  firmly 
to  the  conclusion  that  significant  production  of  small  Kr„ 
clusters  was  inevitable  in  the  multiphoton  experiments  in 
which  the  Kr^"^  transitions  were  observed. 


Fig.  4.  Representation  of  limiting  conditions  for  excitation  of  prompt 
At  Z,-sheil  radiation  horn  Ar  clusters  as  a  function  of  intensity  of 
irradiation  and  cluster  size.  The  Af-shell  limit  is  taken  from  Ftg.  7. 
The  allowed  zone  is  shaded.  Points  P  and  Q  designate  the  crossings 
described  in  the  text  The  conditions  /„(248  nm)  and  I„06S  nm), 
stemming  from  (7),  are  also  indicated.  The  experimental  datum  (LANL) 
of  Lee  et  al.  [2]  is  shown  along  with  the  maximum  intensity  (^5  x 
10**  W/cm-)  used  in  fluorescence  studies  of  Ar  discussed  further  in  [3] 


2.2  Argon  and  Xenon 

Additional  observations  [2^]  with  Ar  and  Xe  can  be  ex¬ 
plained  mutatis  mutandis  and  the  corresponding  diagrams 
appear  for  these  materials  in  Figs.  4  and  5,  respectively.  The 
results  [2]  of  Lee  et  al.  (LANL  data),  which  show  no  prompt 
emission,  are  consistent  in  all  cases  illustrated  for  both  ma¬ 
terials.  Note  the  exclusion  of  the  atomic  region  (n  =  1)  for 
Ar  L-shell  and  Xe  N-shell  emissions  by  the  crossings. 

Again,  however,  the  experiments  examining  the  muiti- 
photon-induced  fluorescence  with  the  pulsed  jet  indicate  very 
different  behavior.  Fluorescence  [3]  ftom  Ar,  observed  with 
a  maximum  intensity  of  irradiation  of  ~5  x  10‘*  W/cm^, 
detected  the  clear  presence  of  radiation  at  ~  165  A  on  the 
well-known  doublet  2s2p*  — >  transition  [16]  of  Ar^'*'. 

This  finding,  as  shown  in  Fig.  6,  is  another  paradoxical  re¬ 
sult,  since  Fig.  7  indicates  that  the  threshold  intensity  of 
Ar’^  production  is  ~10‘*W/cm^,  approximately  20-fold 
higher  than  that  used  in  the  experiment.  Furthermore,  spec¬ 
troscopic  analysis  of  the  plasma  [3]  led  to  the  conclusion 
that  electron  collisions  could  not  have  produced  the  Ar^'*' 
species.  In  this  connection,  we  note  that  the  cross  section 
for  electron  collisional  ionization  of  Ar^'*'  is  quite  small 
[17].  However,  as  shown  in  Fig.  4,  small  cluster  forma¬ 
tion  (n  =  2, 3)  would  allow  the  excitation  of  an  L-sheil 
electron.  Since  the  maximum  experimental  intensity  is  also 
above  the  M-shell  limit,  this  interaction  could  readily  gen¬ 
erate  the  Ar’"*"  species  seen.  Therefore,  small  cluster  forma¬ 
tion  removes  the  anomaly  of  the  Ar®'*’  production  and  the 
observed  radiation  on  the  165  A  doublet  could  then  be  gen¬ 
erated  in  the  manner  previously  described  [3]. 


Cluster  Size  (n) 

■  L-flMI  (5-6  k«V) 

(Oa-IX)  k*V)  0  UVM.  Omtm 

■  x«  N-«MI  (-60^60  .V)  □Forttddm  Zoom 

Ay,  By,  Cy.  6aF2  studM) 

Fig.  5.  Rq)resentation  of  limiting  conditions  applicable  for  excitation  of 
prompt  Xe  N-shell.  Af-shell  and  £-shell  radiation  from  Xe  clusters  as 
a  function  of  intensity  of  iiradiation  and  cluster  size.  The  experimental 
datum  (LANL)  corresponding  to  Lee  et  aL  [2],  which  falls  in  the 
forbidden  zone,  is  shown.  The  boundaries  are  established  as  previously 
described  for  Kr  and  Ar  in  Figs.  2  and  4,  respectively.  The  points  A,, 
and  Cj.  correspond  approximately  to  the  conditions  pertaining  to 
the  data  in  Figs.  1 1, 12.  and  13,  respectively.  See  text  for  details 
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153  159  163  167  171  175 

WAVELENGTH  (A) 

Fig.  6.  Transitions  observed  in  Ar  with  a  pulsed-gas  jet  at  an  intensity 
of  ~5  X  10'*  W/cm*  at  248  nm  as  reponed  in  [3].  The  Ar^'''  2a^2p*- 
2a^2p*  doublet  is  clearly  visible  along  with  two  unidentified  (UI)  lines. 
The  latter  may  be  associated  with  satellite  transitions  in  Ar*'*' 


Classical  MPI  of  AR 


Fig.  7.  Threshold  intensity  for  ionization  shown  as  a  function  of 
resulting  charge  state  Z  calculated  for  Ar  with  the  classical  method  de¬ 
scribed  in  [7],  The  respective  atomic  shells  (L  and  A/)  are  designated. 
The  maximum  intensity  associated  with  leaving  some  A/-shell 

population  is  indicated.  Also  stown  is  the  threshold  for  Ar’"'"  produc¬ 
tion  10'*  W/cm-)  and  the  maximum  intensity  (~5  x  10'*  W/cm*) 
utilized  in  fluorescence  snidies  [3]  of  Ar  using  a  pulsed-gas  Jet  in  which 
the  165  A  2a2p*  —  2a^2p*  transition  in  Ar’"''  was  observed 


A  similar  analysis  can  be  applied  to  Xe  L-shell,  M-sheil, 
and  .V-shell  emissions  which  yields  conclusions  analogous 
to  those  found  for  Ar  and  Kr.  The  corresponding  allowed 
zones  for  prompt  emission  from  those  shells  are  indicated 
in  Fig.  5. 

The  multiphoton-induced  fluorescence  observed  [11,18] 
from  Xe,  particularly  from  levels  involving  4d  excitation. 


Xe 


ItBfl  Xe^*  4<j*535p-»’4d’®5s 


Waveler^tft  (A) 

Fig.  8.  Emission  spectrum  of  Xe  in  the  144-174  A  region  observed 
in  a  pulsed-gas  Jet  at  a  peak  intensity  of  ~5  x  10‘*W/cm^  with 
subpicosecond  248  nm  irradiation.  Feanires  associated  with  Xe^"*", 
Xe*'*',  and  Xe*"*"  are  identified 


Classical  MPI  of  XE 


Fig.  9.  Threshold  intensity  shown  as  a  function  of  resulting  charge 
state  Z  calculated  for  Xe  with  the  classical  method  described  in  [7]. 
The  respective  atomic  shells  (fC,  L,  M,  N,  O)  are  designan-vi  The 
maximum  intensities  /^(Af),  ami  I^[0)  associated  with 

leaving  some  Af-shell.  iV-shell  and  O-shell  populations,  respectively, 
are  indicated.  The  experimental  intensity  Si  5  x  10'*  W/cm-) 
used  in  studies  of  fluorescence  from  Xe  using  a  pulsed-gas  Jet  is 
indicated.  This  intensity  is  Just  sufficient  to  lead  to  observable  Xe’"'’  in 
ion  measurements 


suggests  additional  important  insights.  Transitions  [19-21] 
originating  from  states  involving  a  4d  hole  were  observed 
in  three  species,  Xe^+,  Xe*''’,  and  Xe^.  The  identified 
configurations  of  the  lines  shown  in  Fig.  8  were  4d®5s5p  -♦ 
4d‘05a(Xe^+),  dd’Sp  ^  4<i‘°(Xe*+),  and  4<i«5p  ^ 
4d^(Xe^'*’).  The  maximum  intensity  used  at  248  nm  was 
~5  X  10‘*  W/cm^  which  is  just  sufficient,  as  seen  from  Fig.  9, 
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to  lead  to  observable  Xe®"*"  production  in  ion  charge  state 
measurements  [18]. 

Experimentally,  these  transitions  presented  two  obviously 
anomious  features,  ( 1 )  their  exceptional  strength  and  (2)  the 
nature  of  the  radiating  configurations  detected.  The  Xe^''" 
line  at  165.3  A  was  by  far  the  most  intense  feature  seen  in 
any  of  our  studies  of  fluorescence,  indicating  that  copious 
4(f  excitation  was  certainly  being  produced.  Furthermore, 
consider  the  configuration  [4d*5p]  of  the  excited  state  in 
Xe^"*"  producing  the  array  in  the  ~150A  region.  This  level 
has  two  4d  holes,  but  retains  a  relatively  weakly  bound  5p 
electron.  Production  of  this  level  by  recombination  ftom 
Xe'°'''(4d®)  is  not  possible,  since  the  intensity  was  too  low 
for  its  production  [18].  Collisional  excitation  by  sufficiently 
energetic  electrons  in  the  plasma  formed  is  possible,  but 
calls  for  an  unusually  high  efficiency  given  the  very  high 
observed  intensity  of  the  xenon  transitions.  Furthermore,  it 
was  not  possible  to  reconcile  in  any  reasonable  way  the 
observed  density  dependence  of  the  emission  with  collisions 
arising  from  plasma  electrons.  In  fact,  as  in  the  case  of  Kr®'*’ 
radiation  discussed  above,  the  observed  pressure  dependence 
was  in  far  closer  accord  with  that  associated  with  the 
formation  of  small  rare  gas  clusters  [14,22]  arising  from 
the  complex  hydrodynamics  associated  with  the  flow. 

The  apparently  anomalous  4<i  excitation  implied  by  the 
Xe  fluorescence  would  be  expected,  however,  if  the  cluster 
process  discussed  above  in  relation  to  Ar  and  Kr  is  operative. 
In  Fig.  5  the  triangular  zone  representing  the  allowed  region 
for  prompt  iV-shell  emission,  as  shown  enlarged  in  Fig.  10, 
indicates  that  dimer  (Xcj)  formation  would  be  sufficient 
to  permit  prompt  emission  to  occur  in  the  regions  of  the 
focus  for  which  the  intensity  corresponds  to  the  correct  range 
[/  <  7^(0)].  Emission  from  Xe^"'"  and  Xe*"*"  species  having 
a  4<i  excitation  could  be  expected  to  originate  from  the 


Fig.  10.  Expansion  of  Fig.  S  showing  the  parameters  concerning  4<1- 
subsheil  excitation  in  Xe.  The  experimental  intensity  ~  5  x 
10'*  W/cm-)  used  in  the  studies  of  fluorescence  is  indicated.  The 
limiting  curve  /n(248  nm)  is  shown.  The  O-shell  limit  is  taken  from 
Fig.  9  and  the  allowed  zone  for  prompt  IV-sheil  emission  is  shown. 
The  critical  cluster  size  3t  5)  corresponding  to  4d-shell  excitation 
for  the  prevailing  experimental  conditions  is  indicated 


small  clusters  (2  <  n  <  20)  corresponding  to  the  indicated 
allowed  zone.  Anomalous  ionization  would  also  be  expected, 
which  could  possibly  produce  Xe‘°''‘(4d*),  and  consequently, 
through  recombination  in  the  plasma,  the  Xe®'*'  4<i*5p  — ►  4^ 
lines  at  ~  150  A.  However,  the  information  shown  in  Fig.  10 
suggests  the  possibility  of  another  mechanism  of  production 
for  the  Xe^'''(4<i*5p)  state. 

For  electrons  with  a  kinetic  energy  of  ~100eV,  the 
cross  section  for  4d  ionization  in  Xe  is  unusually  high,  a 
feature  arising  from  the  shape  resonance  involving  the  5p, 
5s,  and  4d  shells  well  known  from  studies  of  photoionization 
and  electron  collisional  ionization  in  that  atom  [23-30]. 
The  measurements  establish  a  maximum  cross  section  of 
(Tj,  =  1.6  X  10“ cm*  for  4<i  electrons  [24].  A  comparably 
high  value  [23,31]  also  holds  for  ions  such  as  Xe^"*".  Since 
the  value  Z  =  7  corresponds  to  the  production  of  Xe*"'',  we 
can  now  evaluate  the  critical  cluster  size  through  (2).  With 
[32]  Tq  =  4.4  A,  we  obtain  =  5,  a  remarkably  small 
value  and  one  that  matches  the  allowed  region  shown  in 
Fig.  10.  Therefore,  since  the  critical  cluster  size  is  so  small 
in  this  case,  the  production  of  single  4d  vacancies  can  be 
nearly  saturated  even  in  relatively  small  clusters,  and  the 
production  of  double  vacancies  could  become  significant. 
Hence,  direct  ^uble  vacancy  production  induced  by  the 
intra-cluster  mechanism  becomes  a  possibility  for  the  origin 
of  the  Xe*'''(4<i*5p)  species. 


2  J  Solid-State  Studies  fBaF2) 

The  results  for  Xe  presented  in  Fig.  5,  with  (a)  the  ne¬ 
glect  of  the  small  difference  in  atomic  number  separating 
Xe(Z  =  54)  and  Ba(Z  =  56),  (b)  the  further  neglect  of  the 
relatively  small  difference  in  atomic  spacings,  and  (c)  the  as¬ 
sumption  that  n  ~  10^  corresponds  to  a  cluster  having  prop¬ 
erties  approximately  those  of  bulk  solid  matter,  enable  us 
to  make  contact  with  data  recently  obtained  in  experiments 
studying  X-ray  generation  from  BaFj  surfaces.  The  spectra 
of  interest  are  shown  in  Figs.  11, 12  and  13.  The  conditions 
pertaining  to  the  data  shown  in  Fig.  11  are  represented  ap¬ 
proximately  by  the  point  in  Fig.  5.  Since  this  point  lies 
at  the  border  of  the  M-shell  allowed  zone,  it  carries  the 
strong  implication  that  anomalous  features  stemming  from 
the  type  of  direct  coupling  discussed  above  may  actually  be 
observable  in  the  spectrum  obtained  with  solid  BaF2. 

Transitions  involving  the  Ba  M-sheil  arising  from  this 
form  of  interaction  would  be  expected  (1)  to  appear  as 
relatively  narrow  features,  (2)  exhibit  observable  satellite 
strucmre  due  to  the  distribution  of  electrons  presumably 
produced  in  the  iV-shell  of  the  system  during  the  interaction, 
and  (3)  display  a  considerably  different  dependence  on 
the  conditions  of  irradiation  than  the  broad  quasicontinuum 
discussed  in  [33]  and  [34],  Indeed,  an  anomalous  feature 
(Aj.  =  13.40±0.05  A)  exhibiting  all  three  of  these  expected 
characteristics  is  seen  in  Fig.  11.  The  transition  appears  to 
be  associated  with  weak  satellite  structure  and  it  was  found 
experimentally  to  be  quite  sensitive  to  the  conditions  of 
irradiation,  a  fcature  illustrated  by  the  comparison  of  Fig.  1 1 
with  Figs.  12  and  13.  Higher  intensity,  corresponding  to 
point  in  Fig.  5,  gives  relatively  less  radiation  at  A^,,  and 
irradiation  with  a  long  (~2  ns)  low  intensity  (<  10*^  W/cm*) 
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WAVELENGTH  (A) 

Fig.  11.  Spectrum  of  BaF,  irradiated  with  subpicosecond  (~270fs) 
ultraviolet  (248  nm)  radiation  at  a  peak  intensity  of  ~8  x  10'*  W/cm- 
observed  with  a  von  Hamos  curved  crystal  spectrometer  on  a  single 
exposure.  Occurring  along  with  the  Ba  iVf-shell  continuum,  the 
f.a(I4.98A)  line  and  an  unidentified  feature  at  13.40  ±0.05  A  are 
seen.  These  data  were  recorded  in  Chicago  by  A.  McPherson  and 
A.P.  Shevel'ko.  We  acknowledge  the  assistance  of  P.G.  Burkhalter 
and  D.A.  Newman  of  the  Naval  Research  Laboratory  in  obtaining  the 
densitometer  trace  of  these  data.  Assuming  n  S'  Ifr*.  the  conditions  of 
these  data  are  represented  on  Fig.  5  approximately  by  the  point 


11. (X)  Wavelength  [A]  16.00 


Fig.  12.  Spectrum  of  BaF,  irradiated  as  in  Fig.  1 1  except  with  an 
incident  intensity  of  ~  10'^  W/cm-.  The  F*''"  transition  appears 
stronger,  but  the  unknown  feature  has  become  relatively  weakened. 
These  data  correspond  approximately  to  point  in  Fig.  5 


a  10  II  12  12  14  15 


WAVELENGTH  (A) 

Fig.  13.  Spectrum  of  BaF;  irradiated  by  1.06  pm  radiation  with 
a  2  ns  pulse  length  at  an  intensity  <  10‘*  W/cm*.  The  M-sheil 
quasicontinuum  is  present,  but  the  unknown  A^  feature  is  absent. 
These  data,  which  were  recorded  at  the  Lebedev  Physical  Institute. 
Moscow,  by  A.  Shevei’ko  and  A.  McPherson  in  April  1993,  correspond 
approximately  to  the  region  below  point  in  Fig.  5 


pulse,  corresponding  to  a  position  far  removed  from  the 
Af -shell  allowed  zone  in  Fig.  5  (region  below  point  C^), 
leads  to  a  vanishing  of  the  line.  Finally,  theoretical 
calculations  performed  by  the  Moscow  group  (A.  Shevel’ko 
and  CO  workers)  indicate  that  corresponds  best  to  the 
3ti‘°4s-  -  3(f4s-4f  transition,  a  transition  whose  nature 
could  be  expected  with  the  proposed  method  of  coupling 
by  3<i  -»  4/  excitation  through  reaction  (2). 

Two  additional  comments  relevant  to  the  data  shown 
in  Fig.  II  can  be  made.  First,  the  unknown  feature  A^  is 
definitely  not  weak,  since  the  9-13  A  continuum  [33,34] 
against  which  it  must  compete,  is  very  strong.  Second,  if  the 
transition  is  associated  with  a  direct  process  of  excitation 
of  the  nature  discussed  above,  then  it  can  only  occur  in 
the  early  stages  of  the  development  of  the  hot  plasma, 
since  the  high  temperature  plasma  would  quickly  lead  to 
full  ionization  of  the  N-shell.  On  the  basis  of  LASNEX 
simulation  [34]  of  the  experiments  [33]  on  Ba  in  which  the 
electron  temperature  is  ciculated,  the  relevant  time  appears 
to  be  approximately  the  initial  ~50-l{X)fs  of  irradiation.  It 
is  further  implied  that  the  excitation  would  also  be  confined 
to  a  region  quite  close  to  the  surface,  perhaps  in  the  first 
skin  depth.  Hence,  such  transitions  might  be  able  to  serve 
as  highly  spatially  localized  probes  of  the  interaction  for 
times  sufficiently  short  that  equilibrium  conditions  cannot 
be  established.  Of  course,  these  observations  only  give  a 
hint  at  the  possibilities,  and  much  further  study  is  needed  to 
verify  or  deny  these  hypotheses  and  possibilities. 


3  Implications  for  X*Ray  Amplification 

Overall,  six  previously  unexplained  or  anomalous  experi¬ 
mental  observations  involving  three  materials  (Ar,  Kr,  and 
Xe)  can  be  accounted  for  by  the  formation  of  small  clus¬ 
ters.  In  addition,  a  possible  connection  with  X-ray  generation 
from  solids  has  been  made.  Furthermore,  it  can  be  similarly 
shown  [35]  that  both  the  observed  [36]  production  of  2<Tg 
excitation  in  Nj,  involving  subsequent  emission  at  55.8  nm, 
and  the  anomalous  production  [37]  of  from  UFg  are 
consistent  with  the  picture  developed. 

In  aggregate,  these  results  have  two  important  implica¬ 
tions.  One  involves  the  possibility  of  developing  a  unified 
picture  of  high-intensity  multiphoton  processes  applicable 
for  sufficiently  short  pulses  to  atoms,  molecules,  and  solids. 
The  second  concerns  the  concept  of  a  molecular  design  for 
systems  optimized  for  efficient  X-ray  production  and  am¬ 
plification.  The  latter  results  from  the  potential  confluence 
of  several  favorable  factors.  They  are  (1)  the  inner-shell 
specificity  of  the  excitation  expressed  in  Figs.  2, 4, 5,  and 
10,  (2)  the  high  efficiency  of  coupling  potentially  indicated 
(n  ~  n^),  and  (3)  the  ability  to  combine  simultaneously  the 
propitious  qualities  of  both  high  and  low  density  in  the  dy¬ 
namic  of  the  medium.  The  selectivity  (1)  and  the  efficiency 

(2)  are  consequences  of  the  fact  that  the  electronic  motions 
involved  in  the  coupling  occupy  only  a  relatively  small  frac¬ 
tion  of  the  available  phase  space.  With  respect  to  the  latter 

(3) ,  the  use  of  clusters  achieves  this  paradoxical  combi¬ 
nation  of  opposites  automatically.  Namely,  clusters  enable 
atomic  excitation  rates  to  achieve  the  extraordinarily  high 
values  characteristic  of  condensed  matter  while  simultane- 
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ously  eliminating,  through  local  rapid  expansion  to  low  par¬ 
ticle  density,  undesirable  collisional  interactions  that  would 
establish  equilibrium  conditions  in  the  excited  material. 

The  interaction  also  leads  naturally  to  the  production  of 
inverted  population  densities.  Consider  the  case  of  Ar.  as 
illustrated  in  Fig.  7,  irradiated  at  an  intensity  I  that  is  just 
below  In  this  circumstance,  seven  electrons  are 

removed  frorn^  the  M-shell,  producing  Ar^^  with  the  config¬ 
uration  ls"2s"2p^3s,  but  no  Ar^'*’  having  the  configuration 
ls-2s-2p®  can  be  formed.  However,  as  shown  in  Fig.  4,  in 
(Ar)„  clusters  with  n  >  3,  an  I-shell  electron  can  be  addi¬ 
tionally  liberated  at  the  same  intensity  with  the  outcome 

that  Ar*'^  can  be  produced,  but  only  with  the  excited  con¬ 
figuration  ls"2s*2p^3s.  Furthermore,  if  we  have  n>n 
every  atom  in  the  cluster  is  produced  in  the  excited  (Ar^'*’)* 
state  and  no  ground  configuration  Ar^"''  species  is  formed. 
Basically,  the  generation  of  the  ground  Ar®'*'  state  is  totally 
blocked,  but  the  channel  leading  to  the  excited  level  is  fully 
open.  Therefore,  in  addition  to  the  selective  and  efficient 
production  of  radiation  on  the  3s  — ►  2p  transition  of  Ar*'*', 
these  dynamics  also  automatically  generate  inverted  popula¬ 
tion  distributions  with  the  capacity  for  amplification. 

Intra-cluster  collisional  ionization  of  the  A/-shell  elec¬ 
trons  by  reaction  (3)  can  potentially  destroy  the  inversion 
by  enabling  the  generation  of  Ar*"^  ls^2s^2p^  ground  state 
species  in  the  example  considered  above.  However,  the 
use  of  heavy  materials,  which  exhibit  shape  resonances  for 
inner-shell  excitation  (e.g.  Xe  4d)  and  which  have  large  cor¬ 
responding  cross  sections  (<7^^  ~  10' ‘^cm^).  appears  to 
sharply  reduce  the  significance  of  such  losses.  Estimates 
of  the  cross  section  for  the  loss  channel  noted  above  (M- 
shell  collisional  ionization),  computed  for  relevant  condi¬ 
tions  using  the  Bethe  approximation  [38],  give  an  estimate 
of  ~  10“  - 10“  '*  cm^,  so  the  shape  resonances  correspond¬ 
ing  to  inner-shell  excitation  can  still  dominate.  Basically,  this 
question  reduces  to  the  issue  of  whether  ionization  cross 
sections,  as  a  function  of  energy,  can  exhibit  significant 
edges  with  respect  to  inner-shell  excitation.  As  shown  in 
Figs.  14  and  15,  the  experimental  evidence  demonstrates  that 
when  the  inner-shell  interaction  involves  a  shape  resonance 
[3 1, 39],  such  edges  can  be  prominent.  However,  the  dynam- 


Fig.  14.  Single  ionization  cross  sections  of  Xe^^  as  a  function  of 
electron  energy  {solid  circles).  Also  shown  are  experimental  results 
of  Griffin  et  al.  {X)  and  Achenbach  et  ui.  (-t-).  The  curve  is  the 
semiempirical  calculation  using  the  Lotz  formula.  Figure  and  caption 
taken  from  (3 1  ]  which  should  be  consulted  for  further  details.  Figure 
used  with  permission 


Fig.  IS.  Absolute  cross  sections  for  single  ionization  of  Ba-^  ions. 
Open  squares  represent  measurements  employing  a  “cold”  Penning  ion 
source  (low  content  of  metastable  ions  in  the  parent  beam);  the  open 
circles  are  measurements  employing  a  “hot”  ECR  ion  source  (producing 
a  substantial  amount  of  metastable  ions  in  the  parent  beam).  The 
statistical  uncertainties  are  smaller  than  the  size  of  the  symbols.  The 
solid  line  is  a  distoned-wave  calculation  for  direct  ionization  of  the  5p 
and  5s  subshells  of  ions  in  the  5s-5p®  ground-state  configuration.  The 
chain  curve  is  a  distoned-wave  calculation  for  direct  ionization  of  the 
5p  and  5s  subshells  of  ions  in  the  5s-5p^5d  metastabie  configuration. 
Figure  and  caption  taken  from  [39]  which  should  be  consulted  for 
further  details.  Figure  used  with  permission 


ics  of  these  collisional  amplitudes,  which  will  occur  when 
the  strong  radiative  field  is  present,  can  only  be  reliably 
tested  by  experiment. 

Auger  processes  can  also  generate  losses.  In  this  case,  the 
Auger  rate  depends  quite  sensitively  on  the  configuration 
[40],  a  property  that  can  be  used  to  suppress  strongly 
the  Auger  decay.  For  the  short  pulse  excitation  (~  100  fs) 
being  considered,  Ar  ions  with  the  excited  state  configura¬ 
tions  Ar^+  (la-2s^2^3s-),  Ar^+ 

(ls-2s-2p^3s3p),  and  Ar*+  (ls^2s-2p*3s)  appear  to  have 
sufficiently  reduced  Auger  rates.  Overall,  the  figures  on 
both  collisional  and  Auger  losses  suggest  that,  in  properly 
selected  systems  undergoing  irradiation  at  an  appropriate 
intensity,  the  loss  channels  may  be  a  small  fraction  (<  10%) 
of  the  total  excitation. 

This  situation  could  lead  to  very  significant  levels  of 
amplification.  If  we  assume  that  the  amplifiying  X-ray 
transition  is  Doppler  broadened,  on  account  of  the  ionic 
velocity  arising  from  the  coulomb  explosion,  ignore  the 
collisional  and  Auger  loss  channels,  and  take  the  cluster  size 
n  =  n^,  it  can  be  shown  [35]  that  the  gain  constant  per 
unit  electron  density  at  wavelength  is  given  by 


9A  _  A, 

n*  SnZ\/2  c 


for  homonuclear  clusters  composed  of  atoms  with  mass  M 
having  an  X-ray  spontaneous  emission  rate  In  (9)  we 
have  assumed,  in  order  to  determine  the  velocity  associated 
with  the  Doppler  width,  that  the  ion  kinetic  energy  is  given 
by  the  expression 


72p2_2/3  p2  - 

=  £J-1—  =  02/3^2  f_  %  _ 
''o  ^0 


(10) 


a  value  that  can  be  shown  to  be  the  upper  bound  of  possibil¬ 
ities  for  a  uniformly  ionized  spherical  cluster  [35].  Table  1 
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Table  1.  Estimated  quMtities  associated  with  the  dynamics  of  muitiphoton  excited  clusters  and  soft  X-ray  amplification.  The  dominance  of  Doppler 
broadening  is  assumed  and  n  =  n^.  See  text  for  definitions  and  details  uoppier 


Matenal 

'"o 

[Al 

Transition 

A,  [A] 

[eV] 

Z 

Anaa 

[W/cm“] 

[cm*] 

[s-'j 

[cm-1 

At 

L»J 

bo 

L(2-2) 

~  50 

250 

7 

3  X  lO'* 

~  1.2  X  10-'* 

~5  X  10* 

1.3  X  10" 

1.9  X  10-1’ 

Kr 

Br, 

Xe 

4.0 

2.3 

4.4 

M(4-3) 

~100 

Af(4-3) 

~100 

.Vf5-i) 

~i50 

210 

184 

~140 

7 

6 

7 

~  2  X  10'* 

~  2  X  10'* 

~10'* 

1-3  X  10-'^ 

1-3  X  10-'^ 

~  1.6  X  10“'* 

4.4  X  ICF 
-1.2  X  lO’ 

~8 

~2.2  X  10^ 

~5 

~10" 

~10" 

5  X  10'" 

1.3-3.9  X  10-'’ 

0.3-1.0  X  10-'* 

4.7  X  10-'* 

E 

2.7 

Ar(5^) 

~150 

~130 

6 

~10'* 

~  1.6  X  10-'* 

~  I 

5  X  10'° 

~io-'* 

Xe 

4.4 

M(4-3) 

~  30 

~700 

25 

6  X  10'* 

~  4x  10-'’ 

~7  X  10* 

~10'2 

4.6  X  10-2° 

U 

2.8 

0(6-5) 

~100 

~200 

10 

~  2  X  10'* 

~  2  X  10-'* 

~1 

~10" 

1-3  X  10-'* 

~4 

0(6-5) 

~I00 

~200 

~I0 

~  2  X  10'* 

~  2  X  10-'* 

~1 

~10" 

~io-'* 

UF6 

--4 

iV(5-4) 

~  30 

~500 

~12 

~  5  X  10'* 

~io-'* 

~4 

o 

7 

~  4  X  10-'’ 

presents  the  estimated  quantities  for  a  selected  group  of  ma¬ 
terials.  It  is  immediately  seen  that  very  favorable  conditions 
for  amplification  {g^  ~  10^cm~')  may  be  achievable  for 
electron  densities  <  10‘®  cm"^  at  wavelengths  as  short 
as  S  100  A  with  the  use  of  subpicosecond  excitation. 


4  Concept  for  Molecular  Design 

A  constraint  governing  the  molecular  design  of  a  system 
optimized  for  X-ray  production  can  be  derived  from  the 
conditions  expressed  in  (1),  (5),  (6),  and  (7).  In  general,  the 
allowed  zone  for  prompt  X-ray  emission  from  the  y-shell 
of  atoms  in  a  homonuclear  cluster  appears  as  a  triangular 
region  of  the  kind  illustrated  in  Fig.  16.  If  we  identify  the 


emission  from  a  homonuclear  cluster  as  a  function  of  intensity  (/) 
and  cluster  size  n.  The  three  constraints  arising  from  (1),  (6).  and  (7) 
are  shown  along  with  the  vertex  3  having  coordinates  (n^,  I^).  The 
diagram  is  similar  to  that  shown  for  the  Xe  iV-shell  in  Fig.  10 


optimization  of  the. system  with  (1)  the  utilization  of  the 
lowest  possible  intensity  and  (2)  excitation  of  the  cluster 
under  conditions  of  maximal  amplification  per  unit  density, 
then  the  vertex  0  defines  that  optimum  subject  to  the 
subsidiary  condition  rig  =  n^.  An  explicit  statement  for  this 
situation  can  be  easily  found.  From  equality  of  (6)  and  (7), 
we  find 


an  expression,  which,  when  combined  with  the  condition 
Ug  =  with  the  use  of  (5),  gives  the  result 


Finally,  if  the  factor  Z  is  represented  by  sening  Ig  = 
from  (1),  then  the  overall  expression 


IZSttq-*  \mc^J  >,lnig 

is  found,  a  statement  which  relates  ( 1 )  the  atomic  properties 
(e^,  Ep,  (7^),  (2)  the  molecular  structure  constant  (rg),  and  (3) 
the  parameters  describing  the  irradiation  (Ig,  A),  Although 
not  fully  establishing  any  particular  molecular  configuration, 
(13)  explicitly  illustrates  the  relative  dependencies  of  these 
important  variables  with  respect  to  the  optimum  condition. 


iTtrl) 


(13) 


5  Conclusions 

A  unified  picmre  of  high-intensity  multiphoton  processes, 
involving  atoms,  molecules,  and  solids,  may  be  developed 
through  the  study  of  clusters.  Of  particular  significance  are 
possible  intra-cluster  processes  that  can  influence  the  level 
of  ionization  produced  and  lead  to  the  excitation  of  inner- 
shell  states  with  the  subsequent  generation  of  prompt  X-rays. 
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The  behavior  of  shape  resonances  [41]  characteristic  of  in¬ 
elastic  electron  scattering  processes,  particularly  for  ionized 
systems  [42]  for  which  rather  little  is  known,  feature  promi¬ 
nently  in  the  dynamics.  Overall,  at  least  six  previously  un¬ 
explained  or  anomalous  experimental  observations  involving 
three  materials  find  interpretations  consistent  with  the  phys¬ 
ical  picture  developed.  It  is  considered  significant  that  this 
synthesis  can  be  achieved  with  the  exclusive  use  of  known 
atomic  and  molecular  properties. 

The  definition  of  a  critical  cluster  size  =  {r\jZa^f 
for  inner-shell  excitation  highlights  the  complexity  of  the 
relationships  governing  the  coupling  by  combining  (1)  a 
static  structural  factor  (rg),  (2)  a  dynamical  consideration 
involving  inelastic  electron  collisional  processes  and 
(3)  the  nonlinearity  of  the  interaction  with  the  strength  of 
the  radiation  field,  which  is  implicitly  represented  by  the 
factor  Z. 

The  general  result  is  that  molecular  systems  and  clusters 
intrinsically  have  properties  making  them  considerably  more 
suitable  than  atoms  for  the  producvan  of  X-ray  emission 
through  the  mechanism  of  multiphoton  coupling  at  high  field 
strengths.  This  outcome  creates  the  possibility  of  designing 
molecules  [43]  optimized  for  efficient  X-ray  production  at 
selected  wavelengths  and  provides  a  general  method  for  the 
generation  of  nonequilibrium  distributions  of  excited  states 
suitable  for  X-ray  amplification.  If  this  can  be  accomplished, 
a  class  of  presently  unknown  molecular  materials  is  likely 
to  emerge.  Presently,  since  the  available  clues  indicate 
that  heavy  atoms  are  strongly  favored,  such  materials  may 
possibly  be  related  to  metallic  clusters  [1],  giant  molecules 
[44]  involving  atoms  such  .  Au.  Pb,  Th,  and  U,  and  systems 
incorporating  rare  eanh  ...oms  [45, 46]  or  transition  metals 
[47]. 
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and  £-Sheil  ( —  6  A)  Transitions 
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Exf>eriments  demonstrating  the  role  of  cluster  formation  on  multiphoton-induced  x-ray  emission  and 
the  scaling  of  this  phenomenon  into  the  kilovolt  range  have  been  performed  on  Kr.  For  the  Kr  M  shell, 
augmentation  of  Kr,  formation  leads  to  a  large  increase  in  Kr’*  (4p—  3d)  emission  (-100  A)  and  the 
appearance  of  a  strong  band  at  -90  A.  The  observation  of  /.-shell  transitions  (-5-7.5  A)  demon¬ 
strates  the  scaling  of  this  phenomenon  into  the  kilovolt  region  and  leads  to  the  conclusion  that  the  in¬ 
teraction  produces  direct  inner-shell  excitation  with  the  emission  of  prompt  x  rays. 

PACS  numbers:  32.80.  Wr,  32.30.Rj 


It  has  been  postulated  [l]  that  a  unified  picture  of 
high-intensity  short-pulse  multiphoton  processes,  embrac¬ 
ing  atoms,  molecules,  and  solids,  appears  possible  through 
the  study  of  atomic  and  molecular  clusters.  Of  particular 
significance  in  this  analysis  are  intracluster  processes  in¬ 
volving  inelastic  electron  collisions  of  the  form 

e-+X^e~-)-e~  +  (X*)^  (I) 

that  can  lead  to  the  production  of  inner-shell  excited 
species  capable  of  prompt  x-ray  emission. 

Essential  predictions  of  the  proposed  model  have  been 
explicitly  tested  in  new  experimental  studies  involving  five 
atomic  shells  [Kr(M),  Kr(Z,),  Xe(N),  Xe(Af),  and 
Xe(Z,)l.  This  Letter  establishes  the  crucial  role  of  cluster 
formation  in  the  generation  of  the  anomalous  emissions 
previously  observed  [2-4J  from  the  Kr  M  shell  (—110 
eV)  and  demonstrates  the  scaling  of  this  new  phenome¬ 
non  into  the  kilovolt  range  (— 800-2 1 00  eV)  through 
measurements  of  Kr  /.-shell  radiation. 

The  detection  of  intense  Kr’'*'  (4p— “3</)  emission  at 
—  100  A  in  earlier  spectroscopic  studies  [2,3]  using 
pulsed-gas  targets  presented  an  obvious  paradox.  Mea¬ 
surements  of  ion  production  [5]  using  tenuous  gas  targets, 
conducted  under  identical  conditions  of  irradiation,  had 
demonstrated  that  no  production  of  Kr’'*'  was  possible 
from  free  Kr  atoms.  Hence,  the  identification  [6,7]  of  the 
Kr’'*'  emissions,  prtxiuced  by  either  a  prompt  or  a  de¬ 
layed  mechanism,  was  fully  anomalous.  However,  if 
sufficient  cluster  formation  were  produced  in  the  hydro- 
dynamic  flow  of  the  jet  furnishing  the  target  [2-4],  the 
cluster  model  showed  that  this  paradox  could  be  entirely 
resolved;  the  anomalous  emissions  found  a  simple  ex¬ 
planation  which  applied  generally  to  all  cases  considered 
[IJ. 

The  multiphoton-induced  production  of  prompt  short 
wavelength  radiation  from  clusters  requires  that  certain 
conditions  be  met  [i].  These  requirements  establish  rela¬ 
tionships  between  the  intensity  of  irradiation  {/)  and  the 
cluster  size  in)  which  define  allowed  zones  for  the  emis¬ 


sion  to  occur.  Such  emission  can  only  take  place  if 
sufficient  kinetic  energy  is  possessed  by  the  colliding  elec¬ 
tron  in  reaction  (1)  to  excite  an  inner  electron.  This  con¬ 
dition,  based  on  a  simple  formulation  of  the  work  done  by 
the  external  field  on  a  free  electron  in  the  cluster,  leads  to 
the  specification  of  two  limiting  intensities  given  by 


/o- 


Sttan^^hri 


/-(A) 


(2) 


in  which  e„  r©,  a.  A,  m,  Xf,  and  c  denote  the  inner-shell 
binding  energy  [8],  the  interatomic  spacing,  the  fine 
structure  constant,  the  wavelength  of  irradiation,  the 
electron  mass,  the  electron  Compton  wavelength,  and  the 
speed  of  light,  respectively.  Hence,  an  inner-shell  j  —  1 
electron  with  binding  energy  e,  can  be  ionized  if  the  in¬ 
tensity  exceeds  a  lower  bound  given  by  either  /©  or  /,(A), 
whichever  is  greatest.  In  addition  to  the  production  of  a 
vacancy  in  the  j—\  shell,  prompt  y— » j  —  \  emission  re¬ 
quires  the  retention  of  at  least  one  electron  in  the  outer  j 
shell  during  the  course  of  irradiation.  This  determines  an 
upper  bound  on  the  intensity  which,  with  the  use  of  the 
tunneling  ionization  model  [9],  leads  to  a  limiting  intensi¬ 
ty  given  by 


/max(;)“c£/(y)/128«^*Z^  (3) 


where  Ep{j)  is  the  ionization  potential  of  the  most  tightly 
bound  y-shell  electron,  Z  the  resulting  ionic  charge,  and  e 
the  electronic  charge.  Hence,  the  allowed  zones  for 
prompt  x-ray  emission  from  atoms  of  a  homonuclear  clus¬ 
ter  generally  apf^r  in  the  form  shown  for  the  Kr  M  and 
L  shells  in  Fig.  1 . 

This  picture  leads  to  the  prediction  that  an  increase  in 
the  cluster  density  and  average  cluster  size  should  lead 
to  both  a  corresponding  intensification  and  a  spectral 
modification  of  the  detected  emi^ion.  Since  Kro  cluster 
formation  can  be  strongly  augmented  by  c(x>iing  the  flow 
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in  the  jet  [10,1 1],  in  light  of  the  form  of  the  allowed  zone 
for  the  t\4  shell  in  Fig.  I,  a  comparison  of  spectra  ob¬ 
served  at  comparable  pressures  (densities)  and  substan¬ 
tially  different  temperatures  should  indicate  the  impor¬ 
tance  of  the  role  of  cluster  species  in  the  production  of  the 
radiating  ions. 

The  experimental  apparatus  used  has  been  described 
elsewhere  [3].  The  laser  [12]  employed  in  the  present 
study,  which  had  a  wavelength  of  248  nm.  a  pulse  width 
of  300  fs,  and  a  power  of  ~0.7  TW,  was  focused  with 
an  //lO  CaF:  lens  into  the  gas  target  and  provided  a 
maximum  focal  intensity  of  (0.5-I.0)  x  io'’  W/cm^  con¬ 
ditions  essentially  the  same  as  in  the  previous  work  [2-4l. 
The  gas  target  was  produced  by  a  high-pressure  pulsed 
valve  fitted  with  a  circular  sonic  nozzle  having  a  diameter 
of  0.5  mm.  The  gas  flow  was  cooled  with  free-flowing  dry 
nitrogen  which  had  passed  through  a  liquid  nitrogen  bath. 
A  thermocouple  attached  to  the  valve  body  measured  the 
temperature  of  the  gas  and  a  grazing  incidence  spectrom¬ 
eter  was  used  to  record  the  spectra  [2-4l. 

The  relevant  spectral  comparison  for  Kr  in  the  72-112 
A  region  is  presented  in  Fig.  2.  The  Kr  spectrum  illus¬ 
trated  in  Fig.  2(a)  was  obtained  with  a  stagnation  pres¬ 
sure  of  1 1 5  psia  at  293  K,  conditions  very  similar  to  those 
of  the  earlier  work  [2,3j.  Figure  2(b)  presents  the  Kr 
spectrum  corresponding  to  a  significantly  reduced  nozzle 
temperature  (238  K)  and  a  stagnation  pressure  (130 
psia).  nearly  the  same  as  that  pertaining  to  the  data  in 
Fig,  2(a). 

Three  principal  differences  distinguish  the  spectra 
shown  in  Figs.  2(a)  and  2(b).  They  are  (I)  a  large  in¬ 
crease  in  the  observed  Kr^"*"  (4/7— »3</)  signal  upon 
reduction  of  the  temperature  (e.g.,  from  —30  to  —450  at 
—  103  A),  (2)  the  appearance  of  a  strong  highly  struc- 


Cluster  Size  (n) 


FIG.  I.  Representation  of  the  allowed  zones  for  production 
of  prompt  Kr  M-  and  L-shell  emission  from  Kr  clusters  as  a 
^notion  of  the  intensity  of  248  nm  irradiation  and  cluster  size. 
The  b^ndanes  outlining  the  allowed  region  are  defined  in  the 
tnaximum  experimental  intensities  used  for  the  M~ 
shell  and  L-sheil  studies  correspond  to  I.OxiO'^  W/cm^  and 
— 8x  10“  W/cm^,  respectively. 


tured  band  at  90  A  at  the  lower  temperature,  and  (3)  the 
large  relative  reduction  in  Fig.  2(b)  of  the  narrow 
features  appearing  in  the  72-83  A  region  of  Fig.  2(a). 
Spectral  data  [6.7]  indicate  that  the  features  observed  in 
Fig.  2(b)  in  the  86-97  A  region  arise  from  Kr‘°''' 
(4/7—  3t/)  transitions  in  the  -86-91  A  range  and  Kr’"'' 
(4j4/7  — 4j3<y)  lines  in  the  -90-97  A  interval.  A  con¬ 
tribution  from  3/7  excitation  in  both  Kr’'*’  and  Kr'®'*'  may 
occur  in  this  same  region.  The  spectrum  also  gives  evi¬ 
dence  for  the  presence  of  a  broad  unresolved  band  in  the 
—  76-85  A  region  composed  of  dense  arrays  l6]  originat¬ 
ing  from  Kr*®'*'  (4^4/7—4^34/)  and  Kr"'*'  (4/7— 3</) 
transitions. 

These  two  signatures,  a  sharply  increased  signal 
strength  and  the  appearance  of  new  categories  of  transi¬ 
tions  and  higher  charge  states,  are  precisely  the  outcomes 
anticipated  from  enhanced  cluster  formation  [1].  One 
expects  that  a  higher  cluster  (Kr„)  density  for  fixed  n 
would  cause  a  corresponding  increase  in  the  emission  and 
that  the  formation  of  larger  clusters  would  lead  to  a  fur¬ 
ther  intensification  of'the  emission  and  to  new  stages  of 
ionization  and  classes  of  transitions.  Since  the  90  A  band 
only  appeared  in  the  cooled  spectrum,  we  believe  that  it  is 
mainly  associated  with  the  formation  of  the  larger  clus¬ 
ters  produced  at  the  lower  temperature  (238  K). 

Figure  2(a)  also  exhibits  the  clear  visibility  of  relative¬ 
ly  weak  narrow  transitions  in  the  -72-83  A  region 
which  can  be  associated  with  4/  —  Id  transitions  [7]  in 
Kr*"*".  In  distinct  contrast  to  the  lines  in  the  85-104  A 
range,  these  features  do  not  grow  dramatically  as  the 
temperature  is  decreased.  They  remain  weak  and  become 
relatively  insignificant  as  the  temperature  is  reduced. 
This  sharply  contrasting  behavior  indicates  the  action  of 
a  different  mechanism  for  their  excitation.  Since  the  two 
spectra  shown  in  Fig.  2  were  generated  under  essentially 
the  same  average  plasma  conditions,  conventional  col- 
lisional  mechanisms,  such  as  recombination,  should  occur 
in  largely  the  same  fashion  in  both  cases.  This  expecta¬ 
tion  matches  the  observed  behavior  of  the  transitions  as¬ 
sociated  with  4/  orbital  excitation  and  we  attribute  those 
weak  features  to  recombination. 

The  cluster  model  predicu  the  existence  of  large  al¬ 
lowed  regions  associated  with  level  structures  in  the  kilo¬ 
volt  range,  a  property  manifest  in  Fig.  I  through  compar¬ 
ison  of  the  zones  associated  with  the  M  and  L  shells. 
Specifically,  the  model  [l]  predicts  that  the  irradiation  of 
Kr  clusters  with  intensities  above  —10'*  W/cm^  for 
nearly  an  arbitrary  cluster  distribution,  should  result  in 
emission  of  prompt  x  rays  from  L-shell  transitions  [7 
13,14]  in  Kr  ions  (—7  A). 

Experiments  testing  this  scaling  behavior  have  been 
conducted  with  a  peak  intensity  of  — 8x  lo'*  W/cm^.  In 
this  case,  the  focusing  was  provided  by  an  //3  off-axis 
parabolic  mirror  and  a  crystal  spectrometer,  using  Kodak 
Industrex  film  and  PET  for  the  4.5-7.5  range,  was  used. 

The  L-shell  emission  spectrum  of  Kr  shown  in  Fig.  3, 
which  involves  transitions  17,13,14]  from  Kr’^  (q-24- 
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FIG.  2.  Comparison  of  Kr  spectra  in  the  72-112  A  region 
produced  at  different  temperatures  by  multiphoton  excitation 
with  subpicosecond  irradiation  at  248  nm  with  a  maximum  in¬ 
tensity  of  (0.5-l.0)xl0'’  W/cm^.  The  vertical  scales  in  both 
si^ctra  are  absolute  values  giving  a  valid  measure  of  the  com¬ 
parative  signal  strengths.  The  spectral  resolution  is  0.9  ±0.1  A. 
(a)  Nozzle  temperature  293  K  and  stagnation  pressure  115 
psia.  UI  designates  unidentified  transitions.  The  base  line  cor¬ 
responds  to  —17  counts/shot,  (b)  Nozzle  temperature  238  K 
and  stagnation  pressure  130  psia.  Unresolved  arrays  from 
Kr  and  Kr"'*'  may  contribute  to  the  signal  in  the  -76-85  A 
region. 


27)  in  the  4.5-7.5  A  region,  was  recorded  with  720  laser 
pulses.  This  spectrum  exhibits  several  salient  properties; 
(Da  group  of  strong  3/—*  21'  lines  in  the  A  re¬ 

gion,  (2)  a  small  set  of  far  weaker  4/, 5/—  2p  transitions 
in  the  4.8-5.6  A  region,  (3)  lines  with  significant  strength 
having  2s  excitation  (e.g.,  F),  and  (4)  a  substantial  num¬ 
ber  of  transitions  involving  satellite  lines  (I4l  [e.g., 
Rk^^'^(K/2p*3/3/'— 2/>^3zf3/3/')  and  Kr^^*(U/2p^2i 
2p^3d2l)  shown  in  the  inset]. 

Several  considerations  bear  directly  on  the  interpreta¬ 
tion  of  these  data.  First,  the  radiating  charge  states  are 
anomalous.  In  analogy  with  the  Kr’’*'  transitions  dis¬ 
cussed  above,  the  presence  of  lines  from  Kr^'*'(^ 
'==24-27)  is  equally  anomalous,  since  the  tunneling  ion¬ 
ization  picture  [5,9]  predicts  no  charge  states  higher  than 
Kr  .  Second,  the  radiative  transitions  are  fast.  The 
lifetimes  [I5]  of  the  3j— 2/?,  3d^2p,  and  3p— 2j 
transitions  occurring  in  the  -6-7  A  region  fall  largely  in 
the  10-1000  fs  range  and  the  2p—  2j(  — 200  eV)  transi¬ 
tion  [7]  has  a  lifetime  [16,17]  less  than  10  ps.  Third, 
recombination  and  plasma  collisionai  processes  are  rela¬ 
tively  slow.  Considering  the  properties  of  the  flow,  the 
relative  positions  of  the  focal  zone  and  the  nozzle,  and  the 
level  of  ionization  (Zs25),  the  maximum  average  elec¬ 
tron  density  [4]  is  —  5  x  10^°  cm  This  gives  a  rate  less 

1812 


FIG.  3.  Kr  Z.-shell  emission  spectrum  produced  by  multipho¬ 
ton  excitation  with  248  nm  radiation  at  a  peak  intensity  of 
— 8x|0“  W/cmA  The  stagnation  pressure  of  the  pulsed  jet 
was  1 15  psia  and  the  gas  temperature  was  243  K.  Transitions 
from  Kr*'*'  (<7  "24-27)  arc  identified.  The  inset,  presenting 
data  from  another  exposure,  details  the  6.5-7.0  A  region  and 
shows  the  prominent  Kr^^'''(F)  and  satellite  lines  ob¬ 

served  in  these  experiments. 


than  10'®  s~'  for  both  recombination  [18]  and  2s— '2p 
collisionai  excitation  [1 6]. 

In  light  of  these  considerations,  the  experimental  obser¬ 
vations  lead  to  the  following  interpretations.  The  five 
weak  transitions  iA-E)  associated  with  Kr^®'*’  all  involve 
transitions  of  the  form  n/—  2p  with  n  >  4.  Since  the  N- 
shcH  limit  [/mBx()V)]  falls  entirely  below  the  allowed  Z.- 
shell  zone  in  Fig.  I,  the  postulated  cluster  mechanism 
cannot  account  for  the  presence  of  these  lines.  In  analogy 
with  the  weak  4/-—  3d  features  occurring  in  Fig.  2(a), 
these  lines  are  attributed  to  recombination  from  Kr^^'*',  a 
species  identified  by  transitions  (//,/,/,/V,0). 

The  strong  lines  observed  in  the  6.3”7.7  A  region  uni¬ 
formly  involve  transitions  with  3s,  3p,  or  3d  electrons. 
Using  the  most  weakly  bound  of  these  (3^),  and  combin¬ 
ing  the  information  in  Fig,  1  with  the  conditions  for  tun¬ 
neling  ionization  [1,9],  we  can  estimate  the  minimum 
cluster  size  necessary  to  observe  Z, -shell  ionization  and 
still  retain  a  3d  electron.  Since  the  maximum  intensity 
consistent  with  the  retention  of  a  34  electron  is  ~2x  I0‘* 
W/cm^  the  Z.-sheil  zone  indicates  that  clusters  with 
n>3  would  suffice.  Evidence  [10,11]  on  Kr^  formation 
indicates  that  an  appreciable  density  with  sizes  at  least  as 
large  as  n  —  1 0  is  formed  under  our  experimental  condi¬ 
tions. 

Lines  F  and  G  involve  2s  excitation  in  Kr^'*’.  Since 
2p—  2j  radiative  decay  is  fast  and  electron  collisionai 
communication  between  the  2s  and  2p  states  is  slow,  the 
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formation  of  Islp^hp  levels  by  recombination 

from  excited  Islp^  Kr"^'*'  ions  would  have  to  be  rapid. 
However,  since  the  plasma  conditions  cannot  support  a 
recombination  rate  comparable  to  the  radiative  decay,  it 
follows  that  the  excited  states  are  directly  pro¬ 

duced  in  the  excitation  of  the  duster.  The  same  reason¬ 
ing  and  conclusion  hold  in  relation  to  the  2d~*  2p  satel¬ 
lite  lines  [e.g.,  Kr^^-*-(|/)  and 
The  comparison  of  the  spectrum  illustrated  in  Fig.  3 
with  corresponding  data  obtained  with  a  Z  pinch  [I3j  re¬ 
veals  a  relative  abundance  of  satellite  lines  in  the  multi¬ 
photon  spectrum  and  significant  reversals  in  relative  line 
intensities.  Both  observations  are  in  natural  accord  with 
the  cluster  model.  The  strong  appearance  of  satellite 
transitions  in  the  multiphoton  spectrum  [e.g.,  Kr^^‘*’(K)], 
at  the  relatively  low  average  electron  density  characteris¬ 
tic  of  this  experiment,  is  again  indicative  of  a  direct 
mechanism  of  production  [l]  in  the  cluster  which  can  be 
associated  with  the  existence  of  the  allowed  L-shell  re¬ 
gion.  In  the  6.3-7.7  A  region,  the  multiphoton  (my) 
^s—2p  transitions  consistently  exhibit  the  anomalous 
property  of  having  a  greater  intensity  than  the  corre¬ 
sponding  2d—2p  lines,  while  the  reverse  holds  for 
the  Z-pinch  (Z)  spectrum  [13].  Specifically,  for  Kr“'*‘ 
R>L  (my),  R<L  (Z)  and  for  Kr”"^  N>J  (my), 

N  <J  (Z).  Recombining  collisionally  dominated  optical¬ 
ly  thin  plasmas  characteristically  produce  emission  with 
transitions  dominating  [I9l  the  corresponding 
/w—»  2/j  lines. 

The  anomalous  strength  of  the  3j  transitions  congenial¬ 
ly  matches  the  expectations  of  the  cluster  picture.  Elec¬ 
tron  loss  can  occur  by  tunneling  ionization  and  intraclus- 
ter  inelastic  electron  collisions.  With  a  maximum  intensi¬ 
ty  of  8  x  |0'*  W/cm^  3j  electrons  cannot  be  removed  by 
the  tunneling  process  [9],  although  the  full  3</ shell  is  ion¬ 
ized  for  intensities  >  2x  10‘*  W/cm^.  Collisionally,  since 
the  2d  ionization  cross  section  [20,21]  is  approximately 
tenfold  the  value  for  the  2s  level,  retention  of  the  2s  elec¬ 
tron  is  again  favored.  It  follows  that  in  the  region  corre¬ 
sponding  to  L-shell  excitation,  particularly  for  intensities 
above  —  lO'*  W/cm^  a  higher  probability  exists  for  re¬ 
tention  of  3j  electrons  than  for  2d.  This  tendency  should 
also  be  enhanced  for  the  higher  charge  states  seen,  as¬ 
suming  a  greater  intensity  for  their  production.  Namely, 
the  2s/2d  ratio  should  be  greater  for  than  for 

Kr  ,  a  feature  that  is  reflected  in  the  multiphoton  data 
(N/J>  R/L). 

In  summary,  experimental  studies  of  multiphoton- 
induced  emissions  in  the  5-7.5  A  region  from  Kr  have  es¬ 
tablished  that  the  cluster  mechanism  demonstrated  in  the 
-80-1 10  eV  range  scales  into  the  kilovolt  zone.  The  ex¬ 
perimental  evidence  supports  the  conclusion  that  multi¬ 
photon  coupling  to  clusters  can  directly  generate  inner- 
shell  vacancies  and  the  production  of  prompt  x  rays.  Fi¬ 
nally,  recent  measurements  [22]  demonstrating  the  multi- 
photon  production  of  Xe(L)  emission  further  corroborate 
this  conclusion. 
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"Multiphoton  X-Ray  Generation  from  Clusters:  Analysis  of  Energy  Transfer 
Scaling" 


29 
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ABSTRACT 

A  preliminary  analysis  is  given  concerning  the  scaling  relations  for  multiphoton 
production  of  inner-shell  vacancies  and  prompt  X-ray  emission  from  molecules  and  clusters. 
The  discussion,  which  is  based  on  recent  studies  of  X-ray  emission  (5-16  A)  observed  from  Kr 
and  Xe  clusters,  indicates  that  energy  deposition  rates  exceeding  -  1  W/atom  are  feasible  in 
appropriately  designed  molecules  incorporating  heavy  atoms. 

1.  Introduction 

A  cluster/molecule  mechanism^  has  been  postulated  for  the  multiphoton  production  of 
X-rays.  Of  panicular  significance  in  this  analysis  are  intra-cluster  processes  of  the  form 

e-  >  X  —  e-  i.  g-  ^  (xy  (1) 

which  can  lead  to  the  excitation  of  inner-shell  electrons  and  prompt  X-ray  emission.  This 
concept  represents  the  cluster/molwule  analog  of  an  interaction  previously  considered  for 
multiphoton  excitation  of  free  atoms”. 

Recent  experiments,^’^  specifically  designed  to  evaluate  certain  essential  predictions  of 
the  cluster  picture,  have  furnished  considerable  supporting  evidence  for  the  main  features  of 
the  proposed  model  These  studies,  which  involved  Krand  Xe,  (1)  established  the  crucial  role 
of  cluster  formation^  in  the  generation  of  soft  X-ray  (-  80-150  eV)  emission,  (2)  verified  the 
scaling  of  diis  new  phenomenon  into  the  kilovolt  (-  800-2100  eV)  region,'^  and  (3)  provided 
important  insights  concerning  the  dynamical  character  of  the  strong-field  coupling  to  the 
clusters.  Although  these  findings  are  only  the  initial  studies  of  these  interactions  and  much 
funher  examination  is  necessary ,  there  currently  exists  a  sufficient  body  of  information  on  which 
to  base  an  analysis  of  the  consequences  of  these  results.  TTierefore,  a  preliminary  assessment 
of  certain  key  scaling  relations  is  the  goal  of  the  discussion  presented  below.  Particular 
emphasis  is  given  to  an  exploration  of  the  upper  limits  of  power  and  energy  transfer,  since  these 
limits  will  govern  the  efficacy  of  multiphoton  processes  in  the  generation  of  bright  X-ray 
sources.  ^ 


Submitted  to  the  ICOMP  VI  Proceedings,  9  September  1993 


30 


2.  Discussion 


2.1.  Backgroumi 

coupiing^o^e'^m?°  Th^m?  concenung  limiting  cross  secdons  for  multiphoton 

coupling  to  ^  atoms.^  The  main  interest  was  the  evaluation  of  the  effective  cross  section  Z 
energy  transfer  in  the  high  intensity  limit  (>  10^9  L  1  e«ecnve  cross  section  for 

considered  only  individual  atoms  and  incoroorated  ^“vy  ^naten^s.  The  analysis 

induced  by  the  incident  wave  consistent  estafaUshin?  the  electronic  motions 

The  outcome  was  the  prediction  that  an  enn^gy  transfer  rate  in  section, 

feasible  with  an  intensity  of  irradiation  lO^^-lO^O  w/cm^  of  0. 1-1.0  W/atom  was 
rate  of  energy  transfer  is  roughly  comparable  to  that  devpion>J^  enormous  value  for  the 
irradiation  of  soHd  sm^  ar  intenS^  developed  during  the  subpicosecond 

encrsy  “?  "•'If 

prop^y  <toip,ed  are  Iwe 

experimental  findings  suggest  that  molecules  can  ^  i^mmde.  Specifically,  the 

behave  ttansiendy  like  solid  matter.  Hence  a  zas  of  thes»"m  “  ^  ^  diem  to 

average  density  L  be  independently  responds  as  a  solid  whose 

molecules  subsequently  undergo  a  rapid  coulomb  explosion,  the  systtnT^Sckly  com!^°iS 
mng  a  tdaavdy  low  avaiagt  panicle  densS-  with  an  nnisoally  eleva^ 


z. 


Conditions  for  Ouster  Excitation  and  X-Ray  Emission 


^  muldphoton  inducfid  production 

Se^fnZrdTr  “  thrmtensirofSaurin" 

size  (n)  which  define  allowed  zones  for  the  emission  to  occur.  Certainly  such  emission  can 
only  take  piace  if  sufficient  kinetic  energy  is  possessed  by  the  colliding  electron  in  Eq  (1)  to 
^v^hT  decnon.  This  condition  Icada  »  rhe  spificanon  o/rwo  linnung  Scndto 


and 


(2) 


®  Uy 


(3) 


m  which  Je.  Tq-  ^  K  m,  ?fg,  and  c  denote  the  inner-sheil  binding  energy,*  the  interatomic 
spacing  of  the  duster,  the  fine  structure  constant,  the  wavdength  of  iiSkdon  the  di^n 
mas,  the  d^n  Compton  wavategih,  and  dm  spend  of  light,f^=ctivdy.  In’  lUs  ^emra, 
m  mnd-sh^  0-1)  deemn  with  binding  dlHgy  ,,  can  be  ionized  if  the  intenaiiy  ~1L<.  i 
la_gl±mnd  given  by  either  Ig  or  Ijj(X),  whichever  is  greatest.  In  addition  to  the  production  of 
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a  vacancy  in  the  O'-O-sheil,  j  —  j-1  emission  requires  the  retention  of  at  least  one 

eiecffon  in  the  outer  j-sheil  during  the  course  of  irradiatiou.  This  Hi.t^ines  an  unner  hmmri 
on  the  intensity  which,  wi±  the  use  of  the  tunneiing  ionization  modei,^  leads  to  a  limiting 
given  by 


4uxW  -  c£^0yi2&i:e‘Z\ 


(4) 


where  ^(,)  th=  .omzaDon  poientiallO  of  the  most  tighUy  bound  j-shell  electron,  Z  the  resuldng 

f  V  f  “T  “““■  ““  I”“P'  X-ny  enntsion 

tarn  atome  of  a  homcnuclear  clutter  appeart  generally  at  a  triangular  re^on  of  the  form  thown 


Fig.  1:  Geneni  representation  of  the  allowed  zone  for  prompt  j-sheil  emiaian  fimn  n  hnn,»«nri^r  , 

faction  of  mteasity  (I)  and  cluster  size  n.  Hie  throe  consaamu  uumg  ih>m  Eqi.  (2).  (3).  and  (4)  .» 
shown  along  with  the  corresponding  mtersecuons  P,  Q,  and  0. 

It  c^  be  readily  shown ^  ±at  an  approximate  estimate  of  the  number  N,  of  i»Tr,>f.H»nT 
occurring  through  reaction  (1)  in  the  cluster  can  be  written  as  the  product  of  the  number  of 
electrons  produced  in  the  cluster  (nZ),  the  atomic  density  in  the  cluster  (ro*^,  the  cross  section 
for  the  inelastic  channel  under  consideration  (<T^j),  and  the  characteristic  seaiw  length  R  of  the 
cluster  given  approximately  by 
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Therefore,  for  reaction  (I), 


(5) 


■  (nZ)c7^r^~’r^n^  =  n^Z—, 


wi±  Z  determined  through  Eq.  (4)  for  the  corresponding*  intensity.  Equation  (6)  is  a 

^  the  extem^fieid  4erience 

only  a  transit  through  the  duster  and  are  subsequently  lost  to  ionization. 

An  important  immediate  consequence  of  Eq.  (6)  is  the  existence  of  a  critical  size 


derin^  by  condiuon  ®  n,  "a^fore,  corresponds  to  an  inner-shell  ionization,  wi± 

J  ^°tn.  Furthermore,  since  the  cridcai 

size  of  the  duster  deprads  inversdy  on  the  square  of  the  atomic  density  0 -2  r^  and 
inversely  on  the  cmbe  of  the  cross  section  for  excitation  (<r^-3),  is  a  direct  ^alogue  of  the 
ear  concept  o  a  macal  mass,  with  replacing  the  neutron  fission  cross  secnon  <r^  and 
Z  p  aymg  the  role  of  the  average  neutron  ^eld  p.  We  note  that  the  definition  of  a  cSrical 
us  er  si^  n  or  inn»-^eil  ^citation  highlights  the  complexity  of  the  relationships  governing 
e  coup  ng  y  combining  (1)  a  static  smictural  factor  (Fq),  (2)  a  dynamical  consideration 
datron  coilisional  processes  and  (3)  the  nonlinearity  of  the  interaction 
ith  the  strength  of  the  radianon  fidd,  which  is  implidtly  represented  by  the  factor  Z. 

2.J.  Synopsis  of  Experimemal  Findings 

nie  generd  picture  presented  above  has  been  spedfically  tested  in  a  r^ent  series  of 
Kpenments.  -  Specifically,  these  studies  have  examined  the  Xe(N)  and  Kr(M)  shells^  in  the 
80-150  eV  region  and  the  Xe(M)  and  Kr(L)  shells^  in  the  800-2100  eV  range.  The  spectroscop¬ 
ic  results  obtMed  in  all  four  casM  were  found  to  be  generally  in  agreement^’**  with  the 
expiations  of  the  duster  picture.*  The  studies  of  the  Xe(N)  and  Kr(M)  shdl  emissions^ 
established  the  crucial  role  of  duster  formation  in  the  generation  of  the  observed  spectra  The 
corresponding  examinations  of  the  Xe<M)  and  Kr(L)  shell  radiation  established  imponant  aspects 
of  the  scaling  of  this  new  phenomenon  into  the  kilovolt  spectral  range. 

3.  PRELIMINARY  ANALYSIS  OF  SCALING  RELATIONSHIPS 

Important  dues  concerning  Ae  mechanism  of  interaction  and  the  scaling  of  the  enerry 
transfer  rate  ^se  from  the  of  the  spectra  in  the  5-16  A  range  involving  the  Kr  L-sheU 

Xe  M-shdl  emissions.  Spouficaily,  Uiese  data  provided  detailed  information  concerning  the 
dynamics  of  the  energy  transfer  and  indicated  (I)  the  significance  of  multiple  transits  of  the 
dnven  electrons  through  the  duster  and  (2)  the  possible  influence  of  coherent  dectron  motions 
on  the  magnitude  of  the  multiphoton  coupling. 
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3.1.  Esrimaie  of  Multiple  Transits 

The  q«iroscepic  evidence*  indicaies  dm  the  endtaiion  is  refined  to  a  shon  time 

of  tSftime  integrity  ofX^uster,  An  esri.^ 

ranee  antitf  nated  from  rh  ^  dme  r  —  10  fs,  a  value  which  fits  leasonahly  well  with  the 
?  dynamics  of  a  coulomb  explosion  of  the  ionized  cluster.  Since  the 
radiauon  is  significantly  shoner  than  this  time,  the  electrons  driven  by  the 

“ve“SSSS  each  half-cycle 

mcaoent  wave.  TTie  number  m  of  such  transits  in  the  time  r  is  given  by 

_  ,  2rc 

T-  (8) 

3.2.  Scaling  o/Energy  Transfer  Rate  and  Possible  Role  of  Coherent  Motions 

It  was  possible  to  establish  a  lower  bound  on  the  eAiu«i  iimni  n  j 

thte  p„  a»m  fan.  the  Kr  darn*  The  value  S 


rate  per  atom  P„  from  the  Kr  data.^^  IHe 

excepuonaily  high  magnitude.  Moreover,  this  exoerimentai  “ 

h“Sf  orute- -- 

ct“ 

In  contmt,  acceptable  agr^ment  with  the  experimental  fijrure  P  was  found  in 
^mp^inson  with  an  e^lier  estimate^  of  the  upper  bomd  of  the  eff^  mSti^ton 

m  fonnulated  m  the  high  mtensity  (>  10^®  wW)  limiL  This  boundlJI 
the  maximum  cross  secuon  <r^  was  the  uncomplicated  <tafemew» 

<T.  •  8«*e  *  3-fi  X  (9) 

In  combination  with  the  mtensity  used  (T^  -  8  x  10«  W/cm2)  in  the  Idlovoll  studies  ^  this 
value  gave  a  corresponding  theoretical  figure  P„  -  <r  I  a  290  V^7’  ■  ^ 

which  appropriately  caceeda  the  eaperimeJml  vaSe.  “ 

The  importance  of  this  finding  is  that  the  d^vation  of  <r_  given  in  Eq.  (9)  included  a 

dectSitiTmotion^  iSuSbylie 
external  field.  Without  this  feature  m  the  interaction,  the  upper  bound  derived  would  have  been 

sTemro^vtiSTv^S!  indicated  value.  The  enhancement  in  the  coupling 

ste.ms  physically  from  construcave  mterference  of  the  driven  eJectrons  rZi  in  the 

.catt^g  pTOcas  Icmiug  to  alonric  eaciiadon.  Roughly,  the  ua4S  SL 
mto  the  analysis  of  the  clusters  leads  to  a  modification  in  Eq.  (6)  of  the  frjnn 

a.  —  Za^.  (10) 

(8)  Md  (10)  both  lead  to  multiplicative  corrections  of  Eq.  (6),  the  conservative 

^“c”7=,t tr  V  "r  f  ='“=•  S=^ore,  iSTSS 

dS  Sds  combined  factor  T  w  ItzZD.  from  the  product  of  Eqs. 
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AT 


2vcZ^a_ 


(II) 


a  result  which  leads  to  a  correspondingly  changed  criticai  duster  si 


size 


Xr^ 


2tcZ^a 


r-^s 


(12) 


Of:,.. 

IS  sufastanoaJIy  raised,  (2)  the  critical  duster  size  is  ver^  n?:,  — 

wavelength  dependence  is  introduced  into  both  of  these  sigmf  csnt 

heavy  atoms  having  high  levds  of  ionitraaon  (Z)  are  oo^tlv  favnr^?  com^sed  of 

magnitude  of  T  can  be  quite  large  very  Furthermore,  smce  the 

experimentally  produced  ^  be^nL^  ^  ^  ^Stbs 


experimentally  producai  can  be  expected  Far  Mamr.u  u  •  x-  ”  .  - 

concoping  Kr  cgd^noo/  Z  =  i 


concerning  i,-5neu  exaianon,-  z  s  25  and  m  a  20  zivinsr  r 
case  of  ±e  Kr  L-shsll,  the  ratio  of  Eq.  (6)  to  Eq.  (12)  is  ’  ^ 

=  r*  «  L25  10*, 

an  enormous  reduction  in  the  critical  size  which  gives  n' 
pracricaily  achieved.  12’ ® 


(13) 

33,  a  value  in  a  range  that  can  be 


J.J.  Molecular  Design  Concept 


Z^a.xc 


bV”“  ^ 


i_Y5w)f^y ., 


I28:ca^J\  mc^JI^^ 


'^rfj±Yi£] 


an  expression  which  relates  (1)  the  atomic  properties  fe  E_m  ^  t  n)  fh^  , 
parameters  (r^,  r),  and  (3)  the  intensity  of  madSToS  ^Alth^S' 
particular  molecular  configuration  Eo  riffi  estafalisfaing  any 

th=«  taponan:  wch  n^^of  *=  ton.1mriB  of 

J.^.  Kr  L-Shell  System 

r  . 

Smcc  WC  see  tto  this  factor  (T)  cm  be  large,  it  cm  also  have  m  enonnous  influence  on  the 
^um^ces  may  allow  ampiiflcation.  ‘  We  now  eaplore  the  c^^e^STrf  to 

S^eS  “mSSf  “PUfloi™  on  3  —  2  annsaons  to  Kr  ions  m  the 

I  assembles  the  pa^«ers  relevant  to  die  exdtarion  of  the  Kr  L-shell  system  We 
wiU  specificaily  consider  the  Kr^6^  transirion^-^.-LS  (ls22s22p<S  —  ls22s22D^drrs  70  i 
whsch  has  a  ^iimve  lifennte'S  of  -  10  fs.  An  eshnie^ Doppl^S^g  L 
^cur  as  a  n^lt  of  the  coulomb  explosion  of  the  duster  indict  that  Doppler^ broadening 
exceeds  de  ra^uve  width,  so  we  will  use  a  stimulated  emission  cross  secS  reflecting  thm 
I^er  wdd.  Using  the  results  devdoped  dsewhere,  ^  and  incorporating  the  enhancement  liven 
by  Eq.  (17).  we  can  wnte  the  .X-ray  gain  constant  g;,  per  unit  dectron  density  n,  ^  ^ 


"*  ^  j  1 7” I 


toni°"a^‘'t“o^“  M'‘T‘'““  rfSSIx.  S 

dat  the  greatest  perm^ble  dectron  density  of  the  medium  is  the  critical  dectron  density 

cor^pondm^  to  wavdength  X,  then  the  maximum  value  of  the  gain  constant 

scales  as  X  ^  a  dependence  heavily  favoring  shorter  wavelengths.  acmevaOle 
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Talile  I: 


Pa«mst«s  to  ths  Kr  (Ishshp^  _  i^hp^3d)  iratmrioa  at  6.7 

^  coiama  sasa^s  that  the  full  ionic  popukdon  reside 


4.  ConcIimoBs 

Recent  erperimenoi  findings  have  indicated  that  Kr  and  Xe  dusters  undersoing 
muldphoton  exatsaon  in  strong  fields  produce  intense  X-tay  emission  from  directly  etdtSonif 

Sg“StS  of 
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A  simple  derivation  in  the  Coulomb  gauge  of  the  nonlinear  Schrodlnger  equation  describing  propaga* 
tion  of  powerful  ultrashort  circularly  polarized  laser  pulses  in  underdense  cold  inhomogeneous  plasmas 
is  presented.  Numerical  solutions  are  given  for  the  two-dimensional  axisymmetric  case  for  both  initially 
homogeneous  plasmas  and  static  preformed  plasma  columns.  These  solutions  account  for  (i)  diffraction, 

(ii)  refraction  arising  from  variations  in  the  refractive  index  due  to  the  spatial  profile  of  the  electron  den¬ 
sity  distribution,  (iii)  the  relativistic  electronic  mass  shift,  and  (iv)  the  charge  displacement  resulting  from 
the  transverse  ponderomotive  force.  The  most  important  spatial  modes  of  propagation  corresponding  to 
(1)  purely  relativistic  focusing  and  (2)  the  combined  action  of  both  the  relativistic  and  charge- 
displacement  mechanisms  are  described.  The  latter  leads  to  the  formation  of  stable  confined  modes  of 
propagation  having  paraxially  localized  regions  of  high  intensity  and  corresponding  paraxially  situated 
cavitated  channels  in  the  electron  density.  It  is  further  demonstrated  that  the  dynamical  solutions  of  the 
propagation  tend  asymptotically  to  the  lowest  eigenmodes  of  the  governing  nonlinear  Schrodinger  equa¬ 
tion.  Finally,  the  calculations  illustrate  the  dynamics  of  the  propagation  and  show  that  the  relativistic 
mechanism  promotes  the  initial  concentration  of  the  radiative  energy  and  that  the  subsequent  charge 
displacement  stabilizes  this  confinement  and  produces  waveguidelike  chaimels. 

PACS  numberis):  52.40.Nk,  42.65 Jx,  52.35.Mw,  52.40.Db 


I.  INTRODUCTION 

The  interaction  of  relativistically  intense  subpi¬ 
cosecond  laser  pulses  with  gaseous  media  has  been  an 
area  of  vigorous  research  for  the  past  several  years.  For 
ultraviolet  wavelengths  on  the  order  of  200-300  nm,  the 
intensity  region  of  interest,  in  which  relativistic  effects  be¬ 
come  important,  lies  above  ~  10‘*  W/cm^.  The  propaga¬ 
tion  of  radiation  in  such  media,  for  intensities  greater 
than  ~  10**  W/cm^,  naturally  causes  strong  nonlinear 
ionization  in  all  matter.  Hence,  the  pulse  itself,  even  in 
regions  where  the  intensity  is  relatively  low  compared  to 
the  peak  value,  removes  many  electrons  [1,2]  from  the 
atomic  or  molecular  constituents  creating  a  plasma 
column  in  which  the  main  high-intensity  component  of 
the  pulse  propagates.  Therefore,  in  a  reasonable  first  ap¬ 
proximation,  the  investigation  of  the  resulting  propaga¬ 
tion  can  be  divided  into  two  separate  and  distinct  areas. 
They  are  (1)  the  atomic  and  plasma  physics  occurring  in 
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the  field  of  the  intense  electromagnetic  wave  leading  to 
the  ionization  and  (2)  the  subsequent  nonlinear  propaga¬ 
tion  of  the  radiation  in  the  plasma  that  is  generated.  The 
work  described  below  concerns  the  latter  issue. 

To  our  knowledge,  Akhiezer  and  Polovin  published  the 
first  treatment  of  high-intensity  electromagnetic  waves  in 
a  cold  plasma  [3].  They  derived  the  equations  describing 
the  propagation  as  a  function  of  the  single  canonical  ar¬ 
gument  (at—kz)  appropriate  for  one-dimensional  wave 
motion,  reduced  the  problem  to  Lagrangian  form  with 
two  integrals  of  motion,  and  presented  either  exact  or  ap¬ 
proximate  solutions  corresponding  to  particular  sim¬ 
plified  special  cases. 

Several  subsequent  treatments  have  been  devoted  to 
the  acceleration  of  charged  particles  in  either  relativistic 
beat  waves  or  in  the  tail  of  a  single  relativistic  pulse. 
Specifically,  Noble  [4]  applied  the  equations  derived  by 
Akhiezer  and  Polovin  to  the  study  of  single-  and  double¬ 
wave  propagation.  In  other  work  [5]  beat  waves  in  hot 
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plasmas  are  described  kineticaJly.  The  excitation  of  par¬ 
ticles  in  cold  plasmas  occurring  in  the  tail  of  a  single  in¬ 
tense  laser  pulse  has  also  been  investigated  [6], 
Importantly,  this  previous  work  [6],  which  considered 
the  case  of  linear  polarization,  established  that  harmonics 
of  the  fundamental  frequency  are  not  produced  at  a 
significant  power,  if  the  plasma  is  underdense  and  the 
group  velocity  of  light  is  close  to  the  speed  of  light. 

Utilizing  the  quasistatic  approximation  for  the  fluid 
equations  of  cold  underdense  plasmas,  Sprangle,  Esarey, 
and  Ting  [7]  derived  the  first  nonlinear,  fully  self- 
consistent  set  of  equations  describing  the  propagation  of 
relativistic  laser  pulses  into  plasmas  for  a  one¬ 
dimensional  geometry.  This  formulation  was  used  to  ob¬ 
tain  insight  into  relativistic  self-focusing,  including  the 
self-consistent  electron-density  profiles,  wake-field  gen¬ 
eration,  optical  guiding,  and  second-harmonic  genera¬ 
tion.  However,  since  this  treatment  was  one-dimensional, 
the  multidimensional  case  remained  unsolved. 

An  important  finding  of  Sprangle,  Esarey,  and  Ting  [7] 
concerned  the  diffractional  erosion  of  the  leading  edge  of 
a  pulse  propagating  in  the  plasma.  However,  in  the  work 
described  below,  we  concentrate  on  the  evolution  of  the 
central  portion  of  the  pulse  as  it  propagates  in  the  plasma 
and  neglect  the  associated  erosion  of  the  leading  edge. 

Extensive  literature  exists  on  the  motion  of  electrons  in 
radiation  fields  of  certain  explicitly  given  forms.  The 
solution  for  the  case  of  plane  monochromatic  waves  is 
known  [8],  and  this  result  has  been  considerably  extended 
in  subsequent  work  [9,10].  Of  course,  the  well- 
established  Volkov  solution  for  the  Dirac  equation  also 
exists  [11].  In  panicular,  Bardsley,  Penetrante,  and  Mit- 
*  tleman  [12]  have  numerically  simulated  the  relativistic 
dynamics  of  electrons  in  a  one-electron  picture  that  in¬ 
cludes  the  effects  of  space  charge  and  the  spatial  distribu¬ 
tion  of  the  radiation  field. 

The  first  treatment  of  the  relativistic  self-focusing  in 
plasmas  was  developed  by  Max,  Arons,  and  Langdon 
[13].  In  addition,  the  general  character  of  the  elec¬ 
tromagnetic  propagation  in  plasmas  has  undergone  con¬ 
siderable  analysis.  Schmidt  and  Horton  [14],  Hora  [15], 
and  Sprangle  et  al.  [16,17]  have  evaluated  the  thresholds 
for  relativistic  self-focusing  using  analytic  methods. 
Especially  germane  is  the  work  of  Sun  et  al.  [18],  who, 
for  initially  homogeneous  plasmas,  derived  the  two- 
dimensional  {r,z)  nonlinear  Schrodinger  equation  govern¬ 
ing  propagation,  including  consideration  of  the  combined 
effect  of  the  relativistic  nonlinearity  and  charge  displace¬ 
ment.  This  work  presented  the  lowest  eigenmode  of  the 
nonlinear  Schrodinger  equation  and  included  numerical 
evaluation  of  the  threshold  of  relativistic  self-focusing, 
the  value  of  which  has  been  approximately  estimated  in 
earlier  work  [14-17].  Finally,  this  analysis  [18]  presented 
the  {r,z)  dynamics  of  the  propagation  for  the  perturbed 
lowest  eigenmode  of  the  nonlinear  Schrodinger  equation 
for  cases  not  involving  spatial  cavitation  of  the  electron 
density.  We  note  that  Kurki-Suonio,  Morrison,  and  Taji- 
ma  [19]  have  also  developed  the  stationary  analytic  solu¬ 
tions  to  this  equation  for  a  one-dimensional  geometry. 

Additional  related  works  can  be  cited  in  this  context. 
Relativistic  self-focusing  and  beat-wave  phase-velocity 


control  in  plasma  accelerators  are  kindred  subjects  [20]. 
Computations  involving  particle  simulations  revealing 
the  initial  process  of  self-focusing  and  subsequent  pon- 
deromotively  driven  electron  motions  have  been  per¬ 
formed  [21].  Other  calculations  analyzing  the  plasma  dy¬ 
namics  and  self-focusing  in  heat-wave  accelerators  [22], 
as  well  as  the  consideration  of  the  nonlinear  focusing  of 
coupled  waves  [23],  have  also  appeared. 

In  prior  publications,  we  have  ( 1 )  investigated  the  gen¬ 
eral  behavior  of  two-dimensional  (r,z)  axisymmetric  rela¬ 
tivistic  self-focusing,  (2)  presented  results  on  the  process 
of  stabilization  of  laser  pulses  in  plasma  columns  [24],  (3) 
described,  with  the  use  of  an  analytical  model  [25],  the 
steady-state  characteristics  in  cavitated  channels  having 
overdense  wails,  (4)  reported  preliminary  results  of  calcu¬ 
lations  that  evaluated  the  combined  action  of  the  relativ¬ 
istic  and  charge-displacement  mechanisms  and  indicated 
the  formation  of  stable  confined  cavitated  modes  of  prop¬ 
agation  [26],  and  (5)  presented  the  experimental  evidence 
of  relativistic  and  charge-displacement  self-channeling  of 
intense  subpicosecond  ultraviolet  radiation  in  plasmas 
[27],  including  specific  comparison  with  the  results  of 
this  computational  model  [27]. 

The  current  work  presents  a  full  description  of  the 
theory  of  nonlinear  propagation  of  intense  axisymmetric 
ultrashort  laser  pulses  in  cold  underdense  plasmas.  In 
this  analysis  we  use  the  term  ultrashort  to  indicate  that 
the  duration  of  the  pulses  r  satisfies  the  inequality 
r,  »r»T,,  with  t,  and  r,  designating  the  response 
times  of  the  ions  and  electrons,  respectively.  We  consider 
both  homogeneous  plasmas  and  preformed  plasma 
columns.  In  this  study,  no  attempt  is  made  to  establish 
consistency  between  the  local  ionization  state  and  the  lo¬ 
cal  laser  intensity,  the  issue  outlined  in  Sec.  I.  For 
sufficiently-low-Z  materials,  such  as  hydrogen  (Hj),  this 
calculation  would  be  unimportant,  since  full  (maximal) 
ionization  would  be  achieved  even  in  rather  low  intensity 
(~10'’-10*®  W/cm^)  regions  [28],  Since  the  main  re¬ 
gime  of  interest  for  this  work  involves  intensities  greater 
than  ~  10‘*  W/cm^,  we  believe  that  the  plasma  condi¬ 
tions  we  have  chosen  for  analysis  are  sufficiently  close  to 
the  true  self-consistent  state  to  be  adequate  for  the  in¬ 
tended  scope  of  this  study. 

The  calculations  discussed  below  have  been  performed 
with  the  specific  goal  of  exploring  the  characteristic  dy¬ 
namics  and  stability  of  the  propagation,  including  partic¬ 
ularly,  the  interplay  of  the  relativistic  and  charge- 
displacement  processes.  Therefore,  in  order  to  illustrate 
this  behavior,  numerical  results  are  presented  that  por¬ 
tray  the  propagation  and  self-focusing  action  as  a  conse¬ 
quence  of  (1)  the  purely  relativistic  nonlinearity  and  (2) 
the  combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms.  Although  the  relativistic 
influence  cannot  be  truly  physically  separated  from  the 
charge-displacement  process,  it  has  been  examined  sepa¬ 
rately  because  this  comparison  provides  insight  on  the 
dynamics  of  the  focusing  action,  specifically  the  process 
by  which  the  charge-displacement  and  the  confined  prop¬ 
agation  in  the  electronically  cavitated  channel  develop. 
Since  both  the  relativistic  effect  and  the  charge  displace- 
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ment  tend  to  locally  produce  a  reduction  in  the  plasma 
frequency  (<u^)  that  is  more  significant  in  the  high- 
intensity  regions,  both  of  these  effects  perturb  the  wave 
fronts  in  a  manner  that  encourages  convergence  of  the 
wave  and  the  formation  of  localized  high-intensity  zones. 
Furthermore,  for  both  the  purely  relativistic  and  the 
combined  cases,  we  have  calculated  the  stationary  eigen- 
modes  of  the  nonlinear  Schrodinger  equation  and  show 
their  relationship  to  the  modes  of  propagation  dynami¬ 
cally  developed  in  the  plasmas.  In  particular,  it  is  shown 
that  the  charge  displacement,  especially  that  resulting  in 
dynamical  cavitation,  has  a  very  strong  effect  on  the  spa¬ 
tial  character  of  the  propagation  [26]. 

The  calculations  indicate  that  the  efficiency  of 
confinement  of  the  propagating  energy  is  potentially 
high,  namely,  that  a  large  fraction  of  the  incident  power 
can  be  trapped  in  a  channeled  mode  by  the  combined  ac¬ 
tion  of  the  relativistic  and  charge-displacement  processes. 
Furthermore,  in  the  asymptotic  regime,  it  is  found  that 
the  channel  is  characterized  spatially  by  an  intensity 
profile  and  electron-density  distribution  corresponding  to 
the  lowest  z-independent  eigenstate  of  the  nonlinear 
Schrodinger  equation.  Overall,  the  principal  conclusion 
of  this  work  is  that,  under  appropriate  conditions,  a  new 
dynamical  mode  of  stable  highly  confined  propagation 
naturally  evolves  for  the  propagation  of  sufficiently  short 
(t,  »r»r, )  pulses  of  coherent  radiation  in  plasmas.  By 
new  dynamical  modes,  we  mean  the  self-channeling  of  the 
radiation  through  the  formation  of  stabilized  electroni¬ 
cally  cavitated  paraxial  modes,  which  result  from  the 
combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms.  Interestingly,  for  ultraviolet 
wavelengths  in  the  (200— 300)-nni  range,  the  power  densi¬ 
ties  naturally  associated  with  these  environments  can  ap¬ 
proach  thermonuclear  values. 

In  Sec.  II  the  underlying  physical  concepts  are  present¬ 
ed  and  the  governing  nonlinear  Schrodinger  equation  is 
derived.  Section  III  presents  the  analysis  of  the  eigen- 
modes  of  this  equation.  Representative  results  of  numeri¬ 
cal  simulations  of  the  two-dimensional  axisymmetric 
propagation  are  given  in  Sec.  IV.  Section  IV  also  con¬ 
tains  a  comparison  with  the  purely  relativistic  case.  Fi¬ 
nally,  the  conclusions  are  summarized  in  Sec.  V. 


trons  from  the  high-intensity  zone.  For  sufficiently  short 
pulses,  only  the  electrons  are  expelled  from  the  laser 
beam  and  the  more  massive  ions,  due  to  their  substantial¬ 
ly  greater  inertia,  are  regarded  as  motionless  [17,18,25], 
A  third  mechanism  is  the  nonlinear  response  arising  from 
the  induced  dipoles  of  the  ions,  but  this  is  generally  small 
and  negligible  [25]. 

(c)  Defocusing  mechanisms,  caused  by  diffraction  from 
the  finite  aperture  and  refraction  by  the  transverse  inho¬ 
mogeneities  in  the  electron  density. 

(d)  Dissipation  of  laser-beam  energy  by  (i)  motion  of 
the  electrons,  (ii)  ionization  of  the  gas  atoms,  (iii)  genera¬ 
tion  of  harmonic  radiation,  (iv)  production  of  inverse 
bremsstrahlung,  (v)  Compton  scattering,  and  (vi)  other 
amplitudes  of  nonlinear  scattering. 

The  present  work  incorporates  four  phenomena:  (1) 
the  nonlinear  response  of  the  refractive  index  of  the  plas¬ 
ma  due  to  the  relativistic  increase  in  the  mass  of  the  free 
electrons,  (2)  the  refractive  index  variation  due  to  the  per¬ 
turbation  of  the  electron  density  by  the  ponderomotive 
force,  (3)  the  diffraction  caused  by  the  finite  aperture  of 
the  propagating  energy,  and  (4)  the  refraction  generated 
by  the  transversely  inhomogeneous  plasma  density  associ¬ 
ated  with  the  formation  of  a  plasma  column.  In  these 
calculations,  preformed  static  plasma  columns  were  used 
in  order  to  approximate  the  radial  distribution  of  ioniza¬ 
tion  that  is  expected,  if  the  incident  laser  pulses  were  pro¬ 
ducing  the  ionization  on  their  temporally  leading  edge. 
Finally,  the  calculations  were  performed  for  a  length  of 
propagation  that  is  much  shorter  than  the  characteristic 
length  for  dissipation  of  the  energy  of  the  pulse  through 
ionization  or  other  modes  of  energy  loss. 

It  should  be  noted  that  the  relativistic  intensities 
characteristic  of  the  phenomena  examined  in  this  work 
can  be  currently  obtained  experimentally  [30-32].  In 
particular,  the  experimental  parameters  characteristic  of 
the  ranges  that  would  apply  to  the  study  of  these  phe¬ 
nomena  are  presented  in  Table  I.  With  the  wavelength 
and  range  of  electron  densities  shown  in  Table  I,  the  plas¬ 
ma  is  always  underdense,  namely,  (Wp/6>)^«1,  where 
Q}=2vc/k  is  the  angular  frequency  of  the  laser  radiation 
and  w^=(4ire^V,/m,  is  the  customary  plasma  fre¬ 
quency. 


U.  GENERAL  CONSIDERATIONS 
A,  Physical  model 

Several  physical  phenomena  [1-7,14-19,24-26,29] 
play  a  role  in  the  nonlinear  dynamics  governing  the  prop¬ 
agation  of  intense  coherent  radiation  under  the  condi¬ 
tions  being  examined  in  this  study.  They  are  the  follow¬ 
ing. 

(a)  The  creation  of  a  plasma  column  by  ionization  in 
the  temporally  early  region  of  the  laser  pulse. 

(b)  The  influence  of  the  spatial  variation  of  the  refrac¬ 
tive  index  arising  from  the  nonlinear  response  of  the 
dielectric  properties  of  the  medium.  Two  mechanisms 
are  related  to  the  electronic  component;  specifically,  the 
relativistic  shift  in  electron  mass  and  the  ponderomotive- 
ly  driven  electron  motion  which  tends  to  displace  elec- 


B.  The  propagation  equation 

Consider  the  propagation  of  an  intense  ultrashort  laser 
pulse  in  a  plasma  with  an  initially  radially  inhomogene¬ 
ous  electron  density,  described  by  the  function  f{r),  so 


TABLE  I.  Experimental  parameters  characteristic  of  the 
ranges  that  would  apply  for  the  study  of  relativistic  and  charge- 
displacement_self-focusing  of  laser  pulses  in  plasmas. 


Peak  intensity 

/=10‘*-10“  W/cm^ 

Pulse  length 

r=»100-l(X)0  fs 

Initial  focal  spot  radius 

To  =1-3  pm 

Wavelength 

A. =0.248  pm  (KrF*  laser  light) 

Target  gas  densities 

p=10“-10“  cm"^ 

Initial  unperturbed 
electron  density 

AT,  =  10*^- 10^'  cm-’ 
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that  =  max/(r)=l.  We  denote  the  vector 

and  the  scalar  potentials  of  the  electromagnetic  fields  as 
A  and  0,  respectively,  and  the  corresponding  electric 
field  as  E.  Let  the  momentum  of  the  electrons,  the 
current  density,  and  the  charge  density  be  denoted  as  P^, 
j,  and  p,  respectively.  We  assume,  for  the  short  timescale 
of  interest,  that  the  ions  are  inertially  frozen  in  space 
[17,18,24,25].  Then, 

□  A  =  c  — (47r/c)j  ,  (1) 

at 


—  —  4trp  , 


(V- A)  =  0  , 


—  +  Vg-V 

at  ® 

Pe=-e 

dA 

at 


-v<i> 


(2) 

(3) 


+c"‘[v,X(VX  A)]  , 
j=-eA^^v,,  , 

r  =  [l  +  |pJV(m,_oc)^]‘^^  . 


(4) 

(5) 

(6) 


In  the  set  of  statements  above,  Eqs.  (1)  and  (2),  are  the 
Maxwell  equations,  Eq.  (3)  is  the  Coulomb  gauge  condi¬ 
tion,  Eq.  (4)  is  the  equation  of  motion  of  the  electrons, 
Eqs.  (5)  are  the  definition  of  the  current  density  and  the 
charge  density  is  the  initial  charge  density,  while 
is  the  dynamical  charge  density),  and  Eqs.  (6)  represent 
the  relativistic  relation  between  the  velocity  v,  and  the 
'  momentum  of  the  electrons.  In  Eqs.  (6),  m^  Q  is  the  rest 
electron  mass,  □  =  V^— c~^a^/3f^. 

It  is  convenient  to  normalize  the  values  in  the  equa¬ 
tions  presented  above  as  follows: 


A  =  (e/m^_oC^)  A,  $={e /m^gC^)<l)  , 

£  =  {e/m^oC^)E  ,  =  p, ,  (7) 

^,=V,/C,  Ne=^e/^,.0> 


both  space  and  time  is  much  greater  than  both  the  plas¬ 
ma  (lire /(Up  q)  and  electromagnetic  (^)  wavelengths. 
This  leads  to  the  validity  of  the  inequalities 


_i  da 

aa 

at 

» 

az 

«lkaUkpa 


where  k^=kl—kj,  kQ=oi/c,  and  kp=cOpQ/c.  We  use 
the  notation  C!}po=4iTe^NgQ/mgQ  specifically  to  denote 
the  unperturbed  plasma  frequency.  With  the  assump¬ 
tions  and  approximations  stated  above,  Eqs.  (l)-(8)  be¬ 
come 


aA  =  k^NgY  'A, 

(10) 

^^<t>  =  k^[Ng-f{r)]  , 

(11) 

o 

1 

> 

(12) 

y=(l^-lp,P)‘^^  Pg  =  A. 

(13) 

The  term  c~'dVif>/dt  is  omitted  in  Eq.  (10),  since  y  and 
(/>,  according  to  Eqs.  (9),  (12),  and  (13),  do  not  have  a 
high-frequency  dependence.  Equations -(11)  and  (12)  re¬ 
sult  in  the  expression  for  the  electron  density, 

iV,=max[0,/(r)-(-A:p“^VV]  •  (14) 

The  logical  function  max(  0,  )  provides  for  the  physical¬ 
ly  obvious  and  necessary  condition  Ng  >  0.  The  analo¬ 
gous  expression  for  the  electron  density  has  been  previ¬ 
ously  derived  by  Sun  et  al.  [18]  for  the  case  /(r)  =  l.  It 
should  be  noted  that  Eq.  (12)  states  the  condition  for  the 
balance  of  the  ponderomotive  and  the  electrostatic  forces 
for  the  relativistic  case.  Through  combining  Eqs.  (9), 
(10),  and  (14),  we  establish  the  equation  for  the  slowly 
varying  complex  amplitude  of  the  vector  potential 
a(r,z,t)  as 

1  a  ^  a 

Ug  ar  az 

+  Vja  +kl{  1  -r~'^m3iz.[0,f(r)+k~'^'7\Y]\a  ) 

=0  .  (15) 


with  the  understanding  that,  henceforth,  the  tilde  sign 
will  be  suppressed.  Using  the  relations  (p^-Vlp, 
=  VtpjV2-p,X(VXpJ  and  Vr  =  V|pJV27,  Eq.  (4) 
becomes 

-|■(p,- A)-v,X[VX(p,- A)]  =  V(<^-r)  .  (8) 

In  the  limit  T»2iT/<Op,  the  expression  p,  =  A  is  approxi¬ 
mately  valid.  As  we  shall  see  below,  this  condition  means 
that  the  electron  response  can  be  regarded  as  adiabatic. 
Furthermore,  assuming  the  vector  potential  to  be  circu¬ 
larly  polarized,  we  write 


In  accord  with  the  previous  assumptions  used  in  deriving 
Eq.  (15),  we  have  neglected  the  second  z  and  t  derivatives 
in  this  expression.  In  Eq.  (15),  is  the  group  ve¬ 

locity  of  light  in  unperturbed  plasma, 
is  its  corresponding  dielectric  constant,  and 
T=(l  -bjap)*^.  It  should  also  be  noted  that  the  electric 
field  vector  E  and  the  vector  potential  A  are  related 
through  the  approximate  relation 


E(r,z,t)^  —ik^  A(r,z,t)  . 


A{r,z,t)  =  j\{e^  +is^  )a{r,z,t)exp[i{o)S  —  kz)]+c.c.  ]  . 

(9) 

Consistent  with  the  statements  made  above,  we  em¬ 
phasize  the  use  of  the  assumption  that  the  pulse  length  in 


The  calculation  of  the  propagating  wave  form  for  the 
central  high-intensity  region  of  the  pulse  is  accomplished 
by  considering  the  solutions  of  Eqs.  (15)  along  its  charac¬ 
teristics.  Changing  the  variable  t  to  g  =  t  -z/Ug,Eq.(15) 
becomes 
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0 

-^a+^(Vla+kj[\-r~'m^[QJ(r)  +  k~'^VlY]\a) 

=0  ,  (16) 

for  which  we  seek  solutions  of  the  form  q  =  const.  We  re¬ 
call  that  certain  solutions  of  Eq.  (16)  for  the  special  case 
/(r)=  1  have  already  been  presented  [18,19,33], 

The  results  of  the  computations  presented  in  this  work 
are  given  with  respect  to  the  variables  (q,r,z). 
Specifically,  the  data  presented  in  the  figures  below  illus¬ 
trate  the  propagation  of  the  radiation  along  the  z  axis  for 
g  =  const.  Equation  (16)  describes  the  two-dimensional 
(r,z)  dynamics  of  propagation  for  coherent  circularly  po¬ 
larized  radiation  in  plasmas  having  an  inhomogeneous 
transverse  (r-dependent)  distribution  of  electron  density. 
The  basic  physical  phenomena  explicitly  embodied  in 
these  equations  have  been  outlined  in  Sec.  II  A. 


III.  STATIONARY  SELF-LOCALIZED  MODES 
OF  PROPAGATION 


Equation  (16),  in  the  case  of  homogeneous  plasmas 
[/(r)=  1  ],  has  axisymmetric  partial  solutions  of  the  form 

a(r,z)=U^^„(kpr)tx^[i(kj/2k)(s —  \)z]  ,  (17) 

where  s  is  a  real-valued  dimensionless  parameter  and  the 
real-valued  function  obeys  the  ordinary  differential 
equation  [18] 

VjCf,„-f-[5+F(£/,^,  )]£/,,  =0.  (18) 


In  this  discussion  the  dimensionless  argument  p=kpr  is 
used  and  Vf  in  Eq.  (18)  denotes  the  Laplacian  with  the 
derivatives  designated  with  respect  to  this  variable  (p). 
The  nonlinear  term  in  Eq.  (18)  is 


^j.n=niax(0,l+Vjyj_J  ,  (19) 


The  natural  boundary  conditions  for  Eq.  (18)  are 


dC/s.„(0) 

dp 


C/„(co)=0. 


(20) 


The  first  condition  assures  axial  symmetry  of  the  solu¬ 
tions,  while  the  second  is  necessary  for  the  property  of 
finite  energy  for  the  solutions.  We  designate  these  solu¬ 
tions  as  “eigenmodes”  and  their  significance  is  explored 
below. 

Consider  initially,  however,  the  eigenmodes  corre¬ 
sponding  to  the  purely  relativistic  case,  functions  which 
are  obtained  from  Eqs.  (18)  and  (19)  by  neglecting  the 
term  Since,  we  are  explicitly  eliminating  the  pon- 

deromotive  potential  and  its  influence  on  the  motion  of 
the  electrons,  in  this  situation  the  electron  density  will 
not  be  self-consistent.  In  this  restricted  case,  Eq.  ( 1 8)  can 
be  rewritten  in  the  form 


_d_ 

dp 


dp 


+  RC/,,„.s) 


-0-1 


dp 


where  f^(  ,5 )  =  (s  /2 ) - ( 1  +  ) ‘ ^  -h  1 . 

The  second  boundary  condition,  stated  in  Eq.  (20),  en¬ 
sures  finiteness  of  energy  only  in  the  situation  where  the 
dynamical  system  defined  by  the  last  two  expressions  has 
three  rest  points,  one  of  them  being  zero,  that  is  for  the 
case  0  <  j  <  1 .  Exactly  as  in  the  classic  theory  of  cubic 
media  [34,35]  for  values  of  the  parameter  s  belonging  to 
this  interval,  the  localized  eigenmodes  constitute  a  count¬ 
able  set  that  is  ordered  by  the  number  of  zeros-n  in  each 
of  its  members  for  finite  values  of  p.  Figure  1  illustrates 
the  first  four  eigenmodes  for  s  =0.95. 

Consider  now  the  eigenmode  corresponding  to  the 
unrestricted  situation  that  includes  both  the  relativistic 
and  charge-displacement  mechanisms,  a  case  which  is 
found  to  be  fully  analogous.  The  eigenmodes  of  the  equa¬ 
tion  describing  the  combined  relativistic  and  charge- 
displacement  processes  for  the  interval  0  <  s  <  1  also 
make  a  countable  set  that  can  be  ordered  in  the  same  way 
as  the  one  associated  with  the  purely  relativistic  case. 
The  zeroth  (lowest)  eigenmode  U^^ip)  with  s=0.95, 
along  with  the  corresponding  electron  density  iV,o(p), 
given  by  the  second  formula  in  Eq.  (19),  are  depicted  in 
Fig.  2.  These  two  functions  have  been  developed  in  ear¬ 
lier  work  [18].  The  first  and  second  relativistic  and 
charge-displacement  eigenmodes  for  the  value  of  the  pa¬ 
rameter  s  =0.95  and  the  corresponding  electron-density 
eigenmodes  N^^ip)  and  N^^iip)  are  depicted  in  Figs.  3  and 
4,  respectively.  These  higher  eigenmodes  exhibit  the  im¬ 
portant  feature  that  cavitation  may  occur  in  the  electron¬ 
ic  component  of  plasma  [see  Fig.  (4)],  even  if  it  does  not 
occur  in  the  case  of  the  lower  eigenmodes  [see  Figs.  2  and 

3]. 

A  prior  study  [18]  treated  the  (r,z)  dynamics  of  the 
evolution  of  the  perturbed  lowest  eigenmodes  for  the  case 
in  which  no  cavitation  of  the  electronic  component  of 
plasma  occurs.  The  lowest  eigenmodes  were  perturbed 
by  multiplying  them  by  constants  in  the  neighborhood  of 
unity.  In  this  case,  the  known  regime  involving  oscillato¬ 
ry  propagation  was  observed.  However,  the  instability  of 
the  numerical  method  applied  in  that  work  [18]  prevent¬ 
ed  the  performance  of  similar  computations  for  cases  in¬ 
volving  cavitation.  Our  calculations  show  that  the  analo¬ 
gously  perturbed  lowest  eigenmodes  V^^ip)  also  propa¬ 
gate  in  the  same  oscillatory  fashion  when  electronic  cavi¬ 
tation  occurs,  the  corresponding  density  of  which  is  given 

^s,o^P^-  It  is  ^so  found  that  the  same  behavior  devel- 


FIG.  1.  Stationary  axially  symmetric  eigenmodes  with 
I  =0.95  corresponding  to  purely  relativistic  self-focusing. 
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FIG.  2.  The  zeroth  (lowest)  axially  symmetric  stationary 
eigenmode  with  s=0.95  corresponding  to  relativistic  and 
charge-displacement  self-focusing:  the  normalized  field  ampli¬ 
tude  distribution  U,_o(p)  and  the  normalized  electron  density 
distribution  N.^olp). 


ops  for  the  perturbed  first  eigenmodes,  which  have  one  or 
two  cavitated  channels  in  the  electronic  density  described 
by  ^spip).  The  results  of  these  simulations  can  be  re¬ 
garded  as  evidence  of  the  (r,z)  stability  of  axially  sym¬ 
metric  lowest  and  first  eigenmodes  against  small  pertur¬ 
bations  in  the  amplitude.  It  should  be  noted,  however, 
that  these  higher  eigenmodes  are  presumably  unstable 
against  small  azimuthal  perturbations,  the  nature  of 
which  violates  the  assumed  axial  symmetry  of  the  distri¬ 
butions. 

The  important  consequence  of  the  character  of  the 
solutions  discussed  above  is  the  fact  that  the  power 


FIG.  4.  The  second  axially  symmetric  stationary  eigenmode 
with  s=0.95  corresponding  to  relativistic  and  charge- 
displacement  propagation:  the  normalized  amplitude  distribu¬ 
tion  and  the  normalized  electron-density  distribution 

NsJp). 

The  explicit  evaluation  of  this  critical  value  is  presented 
below.  As  U^q(p)—*0,  N^q(p)^  \  for  s— »I— 0;  we  have 
[18]  therefore, 

r7.o^.,o- 1 =^,.o(  1  +  ula  1 

*  +  *  -Ujo/2. 

I- 1-0  s^l-O 

Furthermore, 

U,,Q{p)  »  Uo(p), 
j  — i-o 

where  Uoip)  is  the  positive,  monotonically  decreasing 
(reaching  zero  for  no  finite  value  of  p)  solution  to  the 
boundary-value  problem  [18,36] 


Fj.=2rr  fJ’uloip)pdp 

contained  in  the  intensity  distribution  corresponding  to 
the  lowest  eigenmode  U,  Q(p),  unlike  the  case  of  a  cubic 
medium,  depends  on  the  parameter  s.  Namely,  this 
power  decreases  as  s  increases. 

The  infimum  of  P,  by  s,  in  the  internal  0<s  <  1,  is 
called  the  critical  power  (P„)  of  the  relativistic  and 
charge-displacement  self-focusing.  This  power  equals  the 
critical  power  of  the  purely  relativistic  self-focusing 
[18,36]  and 


FIG.  3.  The  first  axially  symmetric  stationary  eigenmode 
with  r=0.95  corresponding  to  relativistic  and  charge- 
displacement  propagation:  the  normalized  amplitude  distribu¬ 
tion  t/,,i(p)  and  the  normalized  electron-density  distribution 
iV^,(p). 


dUo 

-j^(0)=0,  Cfo(oo)=0, 
dp 

with  e=  1  —s.  A  change  of  the  variables,  which  is  a  stan¬ 
dard  procedure  for  treating  the  case  of  the  cubic  non¬ 
linearity  [37],  shows  that  [36] 

Uo{p)={2e)'^'^go{e^''^p)  , 

where  go  is  the  customary  Townes  mode,  i.e.,  the  posi¬ 
tive,  monotonically  vanishing  solution  to  the  following 
boundary-value  problem: 

^i8o-8o+8o=0  , 

dgo 

^(0)=0,  go(«)=0. 

dp 

From  the  definition  of  P„  and  the  relation  between  Uq{p) 
and  gpip),  it  follows  that 

P„=  inf  P=  lim  P, 

0<i<l  j^l-O  * 

=  2rr U\(p)pdp 
=^fJ‘g^(p)pdp=2P„^^  , 

where  P„  c  =2?r/®g^(p)pdp  is  the  critical  power  of  the 
Kerr  self-focusing  in  cubic  media  [34,37]. 

Using  the  resulting  relation  between  the  normalized 
values  of  the  critical  power  of  the  relativistic  and  charge- 
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dijpiacement  self-focusing  Pj-r  the  critical  power  of 
the  Kerr  self-focusing  P„  ^  we  have  the  elementary  result 
that 


Finally,  by  calculating  go  (p)  numerically,  we  find  the  nor¬ 
malized  value  of  P„  with  high  precision  [36]  to  be 

^cr  =  2P  =47r  r°°g^(p)pti/3  =  23.4018  . 

•'  0 

Hence,  we  establish  the  expression  for  the  critical  power 
[36]  as 

^cr,i  =(>nlQC^/e-)  gl[p)pdp(co/a^Q)^ 

=  1.6198X10‘°(£u/ai^,o)^  W  . 

The  constant  factor  involved  in  this  statement  improves 
on  that  given  previously  [18]. 

The  relation  between  C/o(p)  and  go(p)  enables  the  de¬ 
velopment  of  the  asymptotic  expressions  (in  the  case 
e— »-0)  for  both  the  peak  value  of  the  amplitude  of  the  sta¬ 
tionary  eigenmode  and  for  the  radius  of  the  eigenmode  by 
using  the  appropriate  characteristics  of  the  Townes 
mode.  The  results  are 

£/,o(0)-(2£)‘/'go(0) 

and 


for  e— *0. 

Special  computations  were  performed  in  order  to 
determine  how  close  the  values  for  U^q(0)  and 
£/o(0)  =  (2e)‘'^^go(0)  are  in  the  case  where  e«l.  The 
calculations  show  that  is 

2.026X  10"’  and  2.096X  10"*  for  e=  10"’  and  10“*,  re¬ 
spectively. 

IV.  THE  TWO-DIMENSIONAL  CASE 

It  is  convenient  to  treat  Eq.  (16)  numerically  using  the 
normalized  coordinates- 

z^=z/(2krl)  , 
u(r,,2i  )=a^'a(r,z)  , 

where  is  a  characteristic  radius  of  the  initial  intensity 
profile  and  Aq  =  inax|a  (r,0)|  for  the  selected  value  of  g. 
For  simplicity,  we  put  f\(r^)=f[r)  and  omit  below  the 
subscript  1  of  r,  and  z,  for  brevity.  Thus,  the  mathemat¬ 
ical  statements  can  be  expressed  in  the  following  set  of 
equations; 


A 

^-(-/V^u -f  iF(/,(r),l«|^)M  =0,  z>0 

(21) 

u  (r,0)  =  Ug(r),  maxi«g|  =  l 

r 

(22) 

du  _ 

— (0,z)  =  0,  u(oo,z)=0. 
or 

(23) 

The  nonlinear  term  F  is  the  real-valued  operator 


P{/,,|-)  =  a,{l-(l-f-aj^)-'/2 

Xmax[0,/,(r)-baf'V[(l-baj|-)'/2jj  _ 

(24) 

and  the  dimensionless  parameters  a ],  02  are  defined  as 

(25) 

where  /g  is  the  peak  intensity  at  the  entrance  plane 
(z  =0)  of  the  medium.  The  parameter 
Ir~'n^QCo^c^ /{4'rre^)  is  known  as  the  relativistic  intensity 
[15]. 

The  ratio  of  the  power  of  the  beam  Pq  and  the  critical 
power  of  the  relativistic  and  charge-displacement  self- 
focusing  defined  in  Sec.  Ill,  is  an  important  parame¬ 
ter  characterizing  Eqs.  (21)-(24).  In  the  notations  of  the 
present  section,  the  value  of  Fg/^cr  can  be  expressed  in 
the  following  way  [36]: 

Po/PzT~*'<i\‘Zi/B)  f  |ug(r)|Vdr  , 

•'  0 

with  dimensionless  constant  B  given  by  [36] 

B  =  2  fj‘g^(p)pdp=z3.724  51  . 

When  the  initial  transverse-intensity  distribution  is 
Gaussian,  namely,  lug(r)|^=exp(  -r^),  we  have 

Fg/F„=a,a2/(1B)  . 

We  note  that  in  several  other  studies,  the  critical  power 
of  the  relativistic  self-focusing  is  defined  alternatively,  ba¬ 
sically  with  5  =4.0.  In  the  present  discussion,  the  ex¬ 
pression  for  Pq/Pct  involves  a  value  of  B  <4.0,  so  that 
the  corresponding  ratio  Fg/F^  is  raised.  Therefore,  a 
pulse  having  a  Gaussian  initial  transverse-intensity  distri¬ 
bution  and  a  flat  initial  wave  front  with  the  parameters 
^1  “248.6192  and  02  =0.031,  undergoes  self-focusing.  In 
this  case,  the  ratio  Fg/F„,  as  defined  above,  evaluates  to 
1.0347,  namely,  Fg>F„.  However,  using  the  value 
5=4.0  would  give  Fg<F„,  a  statement  contradicting 
the  results  of  the  computations. 

A.  The  initial  conditions 

In  this  section,  we  examine  the  self-focusing  of 
coherent  radiation  for  pulses  having  Gaussian  or  hyper- 
Gaussian  transverse  and  longitudinal  intensity  distribu¬ 
tions  [24]  of  the  form 

/|j=.o=/o(r,r)=/mexp[— (t/r)^'  — (r/rg)^^]  , 

iV,  >2,  N2>2  (26) 

with  r  and  t  being  dimensional.  We  assume  [24]  that 
10**  Wcm"^  with  I*  designating  the  ap¬ 
proximate  value  of  the  threshold  for  rapid  nonlinear  ion¬ 
ization  [1,2].  The  spatial  amplitude  distribution  of  the  in¬ 
cident  radiation,  defined  by  Eq.  (26)  for  the  case  of  a  flat 
incident  phase  front,  is  of  the  form 
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Uo(^)  =  exp(  —  r'''V2)  ,  .  (27) 

The  dimensionless  parameter  ai  in  Eq.  (25),  correspond¬ 
ing  to  the  incident  pulse  intensity  on  the  axis  (r  =Q)  of 
the  beam  /o(t),  is 

a2=/o(r)//,  =  (/„//,)exp[-(r/r)^']  .  (28) 

The  transverse  profile  of  the  plasma  column,  created  by 
the  temporally  leading  edge  of  the  pulse,  is  simulated  by 
the  hyper-Gaussian  function 

f(r)  =  QX^[-{r/r^f^],  N-^tl  ,  (29) 

where  r  is  the  radial  coordinate.  The  aperture  of  the 
plasma  column  can  be  estimated  [24]  with  the  use  of 
the  relation 

•  (30) 

For  example,  in  the  case  of  the  Gaussian  transverse  inten¬ 
sity  distribution  (fV,  =2),  the  aperture  of  the  plasma 
column  for  /®  =  10‘°  Wcm“^  /,  =0.45 X  10^°  Wcm~^ 
and  /o(ro)  =  (0. 1 )/,,  gives  r^^lAlr^.  In  this  situation, 
the  homogeneous  plasma  approximation  /(r)=l  is  valid 
[24].  In  contrast,  it  follows  from  Eq.  (30)  for  plateaulike 
incident-transverse-intensity  distributions  that  the  aper¬ 
ture  of  the  simulated  plasma  column  tends  to  the 
value  of  the  beam  aperture  rg.  Therefore,  in  the  example 
given  above  for  we  have  r,  =  1.25rQ.  Thus,  de- 

focusing  of  the  beam,  which  is  significant  because  of  the 
near  coincidence  of  the  apertures  of  the  laser  beam  and 
the  plasma  column,  must  be  taken  into  account  when  the 
'  evolution  of  a  beam  with  a  plateaulike  incident- 
transverse-intensity  distribution  is  studied  [24]. 


where 

(f){r,^)=  f^F{r,T])d7] 

=a,||-(2/a2)/,(r)[(l+a2^)‘/^-l]j  .  (33) 

Figure  5  illustrates  the  nature  of  the  calculated  solu¬ 
tions  corresponding  to  this  purely  relativistic  case.  The 
parameters  of  the  incident  radiation  and  plasma,  for  the 
examples  presented  in  Fig.  5,  are  A.=0.248  /im, 
/,  =0.45X10^°  Wcm■^  /o  =  }/,  =  3X  10*’  Wcm■^ 
/•o  =  3  /im,  and  N,q=7.5X10^°  cm“^  The  correspond¬ 
ing  values  of  the  associated  dimensionless  parameters  are 
a,  =2.486  192X  10^  and  02  =  7- 

Figure  5(a)  presents  the  result  for  purely  relativistic 
propagation  in  a  homogeneous  plasma  along  the  z  axis 
for  a  pulse  having  an  incident  Gaussian  transverse  inten¬ 
sity  distribution  and  a  flat  wave  front  [iV2  =2  in  Eq.  (27)] 
for  the  value  of  q  defined  in  Sec.  II B,  corresponding  to 
lQ  =  jIr.  Figure  5(b)  presents  the  analogous  graph  for  a 
plateau  like  incident  transverse  intensity  distribution 
[iV2~8  in  Eq.  (27)].  We  conclude  from  these  results  that 
the  solution  is  critically  dependent  upon  the  initial  condi¬ 
tion  represented  by  the  detailed  character  of  the  incident 
transverse  intensity  distribution. 

The  ratio  of  the  beam  power  {Pq)  to  the  critical  power 
(P„)  of  the  relativistic  self-focusing  for  the  given  values 
of  the  parameters  a,  and  a 2  yields  Pq/P„=22.2S1  in  the 
case  of  the  Gaussian  initial  transverse  intensity  distribu¬ 
tion  [1V2~3  in  Eq.  (27),  Fig.  5(a)]  and  Pq/P„~2Q.  168  in 
the  latter  example  of  the  plateaulike  [A'^2~8  in  Eq.  (27), 
Fig.  5(b)]  initial  transverse  intensity  distribution. 

Note  that  for  the  values  of  the  parameters  01,02  and 
incident  wave  forms  studied,  we  have  from  Eq.  (32) 
Pi  <  0.  In  this  situation,  the  following  inequality  is  valid: 


B.  Reiadyistic  self-focusing 


Consider  the  two-dimensional  (r,z)  solutions  of  the  sys¬ 
tem  of  equations  embodied  in  Eqs.  (21),  (22),  and  (23)  for 
the  purely  relativistic  nonlinear  term,  the  form  of  which 
can  be  obtained  from  Eq.  (24)  by  disregarding  the  term 
involving  Vj.  In  this  case,  we  find 

F(/,,^)  =  o,[l-/,(r)(l+02|)-‘^]  . 


The  relativistic  self-focusing  mechanism  prevails  over  the 
charge-displacement  mechanism  outside  of  the  focal  spot 
under  the  conditions  ai»l,  02  — !•  (Discussion  of  the 
conditions  for  the  prevailing  of  the  relativistic  self- 
focusing  can  also  be  found  in  Ref.  [33]).  The  situation 
represented  by  Eqs.  (21),  (22),  and  (23),  with  the  above 
nonlinear  term,  describes  self-focusing  with  a  nondissipa- 
tive  saturation  of  the  nonlinearity.  The  properties  of  the 
solutions  for  this  case  depend  essentially  on  the  values  of 
two  conserved  integrals  given  specifically  by 

P,  =  f”\u(r)l^rdr  (31) 

•^0 


and 


du 

dr 


2 

-<^(^|ui^) 


rdr  , 


(32) 


max|M(r,z)|^>(4/£t,a2)|F2l/'Pi  »  (34) 

r 


namely,  with  respect  to  the  radial  coordinate  r,  the  max¬ 
imum  beam  intensity  has  a  positive  lower  bound  indepen¬ 
dent  of  z.  This  conclusion  can  be  established  in  the  same 
way  that  Zakharov,  Sobolev,  and  Synakh  [37]  demon¬ 
strated  the  analogous  result  in  their  earlier  work  on  self- 
focusing.  Therefore,  a  powerful  relativistic  beam,  which 
self-focuses  in  a  homogeneous  nonabsorbing  plasma 
{fi—0.  Ref  [24]),  results  in  a  field  distribution  represent¬ 
ing  a  pulsing  waveguide  when  Pi  <p.  Figures  5(a)  and 
5(b)  explicitly  illustrate  the  formation  of  such  a  regime. 
These  pulsing  waveguides  consist  of  alternations  of  ring 
structures  and  focal  spots  on  the  axis  of  the  beam.  The 
power  confined  in  these  complex  modes  represents  ap¬ 
proximately  50%  and  90%  of  the  total  incident  power  for 
the  cases  depicted  in  Figs.  5(a)  and  5(b),  respectively.  It 
should  be  noted  that  oscillating  periodic  solutions  of  this 
type  have  also  been  presented  by  other  workers  [17]. 

In  the  case  P2  -()>  ^  inequality  comparable  to  (34) 
cannot  be  obtained.  Following  the  method  of  Zakharov, 
Sobolev,  and  Synakh  [37],  we  simulate  the  case  Pj  -  by 
considering  a  beam  initially  focused  or  defocused  by  a 
lens  at  the  entrance  plane  (z  =0)  of  the  medium.  Let  the 
focusing  (or  defocusing)  length  of  this  lens  be  P  =  kr^Ff 
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FIG.  5.  Purely  relativistic  propagation  with  /q  =  3  X  10” 
W/cm%  ro  =  3  fim,  A.=0.248  fim,  and  N,_q  =  7.5X  ltf°  cni“\  la) 
The  formation  of  the  pulsing  waveguide  regime  in  the  case  of 
the  relativistic  seif-focusing  of  a  pulse  with  a  flat  incident  wave 
front;  Gaussian  initial  transverse-intensity  distribution  [lVi=2 
in  Eq.  (27)],  homogeneous  plasma,  (b)  The  formation  of  the 
pulsing  waveguide  regime  in  the  case  of  the  relativistic  self- 
focusing  of  a  pulse  with  a  flat  incident  wave  front;  hyper- 
Gaussian  initial  transverse-intensity  distribution  [^2  =  8  in  Eq. 
(27)],  homogeneous  plasma,  (c)  The  single-focus  regime  in  the 
case  of  the  relativistic  self-focusing  of  an  initially  focused  pulse 
[Rf^Rf^o/2  in  Eq.  (35)];  Gaussian  incident-transverse-intensity 
distribution  [^2=2  in  Eq.  (35)],  homogeneous  plasma,  (d)  The 
formation  of  the  quasistabilized  regime  in  the  case  of  the  rela¬ 
tivistic  self-focusing  of  a  pulse  with  a  flat  incident  wave  front  in 
a  plasma  column  [N-j  =8,  r,  =ro  in  Eq.  (29)];  hyper-Gaussian 
initial  transverse-intense  distribution  [N2=i  in  Eq.  (27)]. 


( R f  dimensionless).  The  condition  R^>0  signifies  that 
the  pulse  is  initially  focused  and  <  0  indicates  that  it  is 
initially  defocused.  Then,  the  corresponding  initial  con¬ 
dition  can  be  written  as 

Uo(r)=exp[  — r^V2  +  />‘/(2i?y )],  .  (35) 

Furthermore,  let  Pi  =0  for  |/?y  |  =R f^Q.  We  note  that  the 
case  ?2  >  0  corresponds  physically  to  a  high  degree  of  the 
initial  focusing  or  defocusing:  g  >  I  1  >  0.  Moreover, 

when <  |/{^ I  <  00 ,  <  0. 

Figures  5(a)  and  5(c)  display  the  calculated  intensity 
distributions  for  the  case  of  relativistic  self-focusing  in 
homogeneous  plasma  of  beams  with  Gaussian  incident- 
transverse-intensity  distributions  and  values  of  Rf  given 
by  Rf  =  ’X  (flat  wave  front)  and  R f^o/2,  respectively.  A 
more  detailed  study  of  the  dynamics  of  the  transition 
from  the  pulsing  waveguide  regime  [see  Fig.  5(a)],  to  the 
single-focus  regime  [see  Fig.  5(c)  representing  the  propa¬ 
gation  of  an  initially  sufficiently  sharply  focused  beam)] 
can  be  found  in  Ref.  [38].  Importantly,  the  computations 
show  that  the  value  .P2  =0  is  not  the  threshold  separating 
these  two  regimes  of  relativistic  self-focusing.  This  tran¬ 
sition  occurs  as  the  first  focus  gains  power  and  is  shifted 
closer  to  the  entrance  of  the  medium  (z=0),  while  the 
remaining  foci  are  shifted  in  the  opposite  direction  and, 
in  the  limit  of  large  displacement,  become  diffused.  Note 
that  single-focus  regimes  of  propagation  have  also  been 
observed  by  Sprangle,  Tang,  and  Esarey  [17], 

Defocused  beams  can  evolve  in  a  different  fashion.  Ini¬ 
tially  sufficiently  sharply  defocused  pulses  {R^<0, 
|Ry|  ^Rf  Q/2,  P2>0)  propagating  in  the  relativistic  re¬ 
gime  monotonically  diffiise  on  the  radial  periphery  and 
do  not  exhibit  the  phenomenon  of  self-focusing. 

The  detailed  spatial  character  of  the  plasma  column 
can  have  a  strong  influence  on  the  evolution  of  the  propa¬ 
gation.  Figure  5(d)  iUustrates  the  relativistic  propagation 
in  a  plasma  column  along  the  z  axis  for  a  pulse  corre¬ 
sponding  to  Ni^S  in  Eq.  (27).  The  transverse  profile  of 
the  plasma  column  is  given  by  /(r)  as  defined  by  Eq.  (29) 
with  A^3=8  and  r,  =rQ.  The  comparison  between  Figs. 
5(b)  and  5(d)  demonstrates  that  the  defocusing  of  a 
hyper-Gaussian  beam  in  a  plasma  column,  with  an  aper¬ 
ture  close  to  the  radius  of  the  beam,  fundamentally  alters 
the  spatial  dynamics  of  the  propagating  energy  [24].  The 
defocusing  causes  a  fraction  of  the  beam  to  spread  away 
from  the  column,  but  the  remaining  energy  of  the  beam 
resolves  into  a  state  that  balances  the  relativistic  self- 
focusing,  defocusing,  and  diffraction.  The  power  trapped 
in  this  quasistabilized  state  so  formed  is  approximately 
25%  of  the  incident  power. 

Previous  analytical  estimates  [38]  have  shown  that  the 
relativistic  self-focusing  length  is  minimal  when  the  in¬ 
cident  pulse  intensity  on  the  axis  of  the  beam  satisfies  the 
condition  /git) =2/,.  Moreover,  specific  computations 
also  show  that  this  inference  remains  valid  for  the  case  of 
the  relativistic  self-focusing  of  pulses  with  flat  initial 
phase  fronts  and  Gaussian  incident-transverse-intensity 
distributions  in  homogeneous  plasmas. 

The  locus  of  the  first  focus  for  the  case  of  the  relativis¬ 
tic  self-focusing  has  been  presented  previously  [24].  In 
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that  study,  the  beam  initially  had  both  a  flat  phase  front 
and  a  Gaussian  transverse-intensity  distribution.  The  pa¬ 
rameters  were  the  same  as  those  applying  to  Fig.  5{a),  and 
the  duration  of  the  pulse  corresponded  to  r=0.5X  10~'^ 
s.  The  minimal  z  of  the  locus  is  reached  when  /o(f)=2/,.. 
If  the  maximum  intensity  on  the  axis  of  the  beam  is  such 
that  >  21^,  then  the  trajectory  of  the  first  focus  exe¬ 
cutes  a  path  that  has  three  reversal  points.  Two  of  them, 
corresponding  to  an  identical  value  of  z,  are  due  to  the 
condition  /o(t)  =  2/,..  The  central  one,  corresponding  to 
a  greater  value  of  z,  arises  when  /o(r)=/„.  In  principle, 
this  behavior  makes  it  possible  to  distinguish  the  process 
of  relativistic  self-focusing  from  that  associated  with  the 
Kerr  nonlinearity,  since  the  latter  results  in  a  locus  of  the 
first  focus  having  only  a  single  reversal  point. 

C.  Relativistic  and  charge-dispiacement  seif-channeling 

This  section  is  devoted  to  the  description  of  the  seif¬ 
channeling  occurring  when  both  the  relativistic  and 
charge-displacement  mechanisms  are  included  in  the  in¬ 
teraction.  This  situation  is  described  by  Eqs.  (21)-(24) 
and  includes  the  important  nonlinear  term  involving  Vfy 
appearing  in  Eq.  (24).  The  numerical  results  show  that 
sufficiently  intense  short  duration  (T,»r»rJ  axisym- 
metric  pulses  readily  undergo  self-channeling  in  plasmas 
over  a  rather  broad  range  of  conditions.  Moreover,  it  is 
found  that  a  large  fraction  of  the  total  incident  power  of 
the  beam  can  be  trapped  in  a  stabilized  mode  confined  to 
the  axis  of  propagation. 

A  specific  example  is  informative  in  representing  the 
general  behavior  exhibited  by  the  propagation  in  the  re¬ 
gime  in  which  the  influence  of  relativity  and  charge  dis¬ 
placement  are  both  significant.  These  results  are  exhibit¬ 
ed  in  Fig.  6.  In  the  case  of  propagation  of  an  initially 
Gaussian  transverse  wave  form  incident  in  a  homogene¬ 
ous  plasma,  with  the  parameters  A. =0.248  /tm, 
Wcm~^  rQ  =  3/xm,  and  N,q  =  7.5 
X10^°  cm“^  (a,  =248.6192,  02  =  },  Po/P„=22.252), 
the  numerical  computations  show  that,  as  soon  as  the 
first  focus  on  the  axis  of  propagation  is  formed,  electronic 
cavitation  occurs.  Specifically,  this  leads  to  complete  ex¬ 
pulsion  of  the  electronic  component  of  plasma  from  the 
paraxial  domain  [Fig.  6(b)].  This  process  results  in  a 
quasistabilized  cavitated  channel  in  the  electron  distribu¬ 
tion  which  extends  along  the  entire  axis  of  propagation 
past  the  location  of  the  initial  focus  [26].  We  note  that 
some  of  the  lowest  stationary  solutions  corresponding  to 
the  relativistic  and  charge-displacement  problem,  includ¬ 
ing  cavitation,  were  developed  by  Sun  et  ai.  [18].  The 
analysis  shows  that  the  first  focus  involves  about  45%  of 
the  total  incident  power  of  the  propagating  energy.  A 
fraction  of  the  remaining  power  is  dissipated  through 
diffraction  on  the  periphery,  while  another  component  is 
temporarily  involved  in  the  formation  of  a  pulsing  ring- 
shaped  structure  [Fig.  6(a)].  Subsequent  energy  exchange 
between  the  ringlike  feature  and  the  paraxial  zone  is  ob¬ 
served,  and  as  a  consequence  of  this  interaction,  a  certain 
part  of  the  energy  of  the  pulsing  ring  is  diffracted  away 
while  the  remaining  power  joins  the  paraxial  domain. 

In  finer  detail,  the  following  aspects  of  the  evolution  of 


the  radiative  energy  are  also  revealed  by  the  calculations. 
As  described  above,  after  the  formation  of  the  first  focus, 
considerable  power  is  transferred  from  the  paraxial  focai 
zone  to  the  ring-shaped  feature,  which,  at  this  stage  in 
the  evolution  of  the  pulse,  contains  approximately  68% 
of  its  total  initial  power.  This  intense  ring,  which  spreads 
away  from  the  paraxial  domain,  produces  a  correspond¬ 
ing  ring-shaped  cavity  in  the  electronic  distribution  [see 
Fig.  6(b)].  The  refraction  resulting  from  this  strongly 
perturbed  electron-density  profile,  together  with  the  rela¬ 
tivistic  self-focusing  mechanism,  causes  the  wave  to  re¬ 
turn  energy  to  the  core  of  the  beam.  Thus,  the  charge 
displacement  produces  a  potent  additional  self-focusing 
action,  which  leads  to  the  formation  of  a  confined  paraxi¬ 
al  mode  of  high  intensity  stabilized  along  the  axis  of 
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FIG.  6.  The  seif-channeling  of  a  pulse  with  a  Gaussian  initial 
transverse-intensity  distribution  [^2=2  in  Eq.  (27)]  and  a  flat 
incident  wave  front  in  initially  homogeneous  plasma  in  the  case 
of  the  relativistic  and  charge-displacement  propagation  with 
/o  =  3X10'’  W/cm%  ro  =  3  pm,  A.=0.248  pm,  and 
l^o—l^e.o  —  7.5X  10^°  cm“^  (a)  The  distribution  of  the  normal¬ 
ized  intensity,  (b)  The  distribution  of  the  normalized  electron 
density,  (c)  Radial  dependence  of  the  asymptotic  solutions  for 
the  normalized  amplitude  [/,(r)//o]''^  and  the  normalized  elec¬ 
tron  density  Ns{r)/No  for  s  =0.554. 
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propagation.  This  phenomenon  of  confined  propagation 
is  designated  as  self-channeling.  It  should  be  noted  that 
the  calculations  of  the  electrostatic  energy  associated 
with  the  charge  displacement,  which  is  given  by  the  ex¬ 
pression  ff'=4~'f^‘E^rdr,  with  the  electrostatic  field 
E,  defined  by  the  equation  VE^  =  -4irp  and  r,  designat¬ 
ing  the  radial  extent  of  the  channel,  show  that  this  energy 
can  be  relatively  small.  Specifically,  for  the  conditions 
represented  in  Fig.  6(b)  [z  =95.4  /tm],  the  electrostatic 
energy  (  W)  accounts  for  only  0.18%  of  the  total  energy 
of  the  laser  radiation  per  unit  length. 

The  essential  finding  of  these  calculations  is  that  the 
combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms  produces  a  strong  tendency  for 
the  generation  of  spatially  highly  confined  modes  of  prop¬ 
agation  which  are  stabilized  along  the  axis  of  propaga¬ 
tion.  Furthermore,  the  study  of  a  range  of  other  cases  in¬ 
dicates  that  these  modes  are  exceptionally  stable  and  that 
a  considerable  fraction  of  the  incident  power  can  be 
confined  in  them.  The  result  is  the  controlled  generation 
of  a  very  high  peak  intensity  in  these  channeled  modes 
with  values  reaching  ~  10^’  W/cm^  for  the  range  of  con¬ 
ditions  studied. 

Important  characteristics  of  the  asymptotic  behavior 
of  these  confined  modes  have  also  been  established.  It 
has  been  shown  that  the  distribution  of  the  amplitude 
u(r,z),  for  large  values  of  z,  tends  asymptotically  to  the 
lowest  eigenmode  of  the  nonlinear  Schrodinger  equation 
(see  Sec.  III).  Specifically,  for  the  example  discussed 
above,  the  computations  have  demonstrated  that  the 
asymptotic  radial  amplitude  distribution  corresponds  to 
the  lowest  eigenmode  with  5=0. 554.  In  this  case,  the 
asymptotic  intensity  distribution  I,(r)  =  Ulo(r)  contains 
46%  of  the  total  incident  power. 

normalized  asymptotic  field  amplitude 
[/j(r)//o]  and  the  corresponding  normalized  asymp¬ 
totic  plasma  electron  density  N,(r)/No  are  depicted  in 
Fig.  6(c).  It  should  be  noted  that  the  profiles  of  the  inten¬ 
sity  I{r,z)  and  the  electron  density  N(r,z),  obtained  as 
the  results  of  the  dynamical  calculations  of  the  propaga¬ 
tion  for  z=:900  /xm,  differ  from  I,(r)  and  N,(r)  for 
r  ~0. 554  by  much  less  than  1%.  In  addition,  we  observe 
that  the  energy  of  charge  displacement  for  =2.5  /xm  in 
the  asymptotic  state  accounts  for  a  fraction  of 
9.45X10  *  of  the  total  energy  of  the  beam  per  unit 
length.  We  observe  that  this  tendency  for  a  solution  of  a 
nonlinear  Schrodinger  equation  involving  a  saturable 
nonlinearity  to  converge  to  the  lowest  stationary  solution 
was  originally  discovered  by  Zakharov,  Sobolev,  and 
Synakh  [37]. 

For  the  range  of  parameters  studied,  the  calculations 
clearly  show  that  the  charge  displacement  has  a  very 
strong  influence  on  the  character  of  the  propagation  after 
the  first  focus  is  formed.  The  pulsing  intensity  structure, 
consisting  of  alternating  foci  on  the  axis  of  propagation 
and  peripheral  focal  rings,  which  is  the  usual  behavior  for 
the  purely  relativistic  self-focusing  [see  Figs.  5(a)  and 
5(b)],  is  converted  into  a  stabilized  and  uniform  channel. 
Collaterally,  a  stabilized  cavitated  channel  in  the  electron 
density  is  also  formed.  We  remark  that  the  periodic  in¬ 


tensity  structure  characteristic  of  the  relativistic  self- 
focusing  occurs  for  great  values  of  z  and,  therefore,  may 
be  regarded  as  the  corresponding  asymptotic  solution  of 
the  purely  relativistic  case.  The  charge  displacement 
leads  to  the  formation  of  the  asymptotic  amplitude  distri¬ 
bution  represented  by  the  corresponding  lowest  eigen¬ 
mode  U^  Q(r)  instead  of  this  pulsing  structure. 

Computations  have  been  also  performed  for  the  propa¬ 
gation  of  incident  plateaulike  wave  forms  with  flat  in¬ 
cident  phase  fronts  in  both  homogeneous  plasmas  and 
plasma  columns.  Additional  calculations,  have  also  ex¬ 
amined  the  behavior  of  focused  Gaussian  and  plateaulike 
incident  wave  forms  as  well  as  defocused  Gaussian  in¬ 
cident  wave  forms  in  homogeneous  plasmas  (see  Figs. 
7-11).  It  is  found  that  the  main  features  described 


FIG.  7.  The  self-channeling  of  a  pulse  with  a  hyper-Gaussian 
incident-transverse-intensity  distribution  [1V2  =  8  in  Eq.  (27)] 
and  a  flat  initial  wave  front  in  initially  homogeneous  plasma  for 
the  case  of  a  relativistic  and  charge-displacement  propagation 
with  /o  =  3X10'’  W/cm^  ro=3  fim,  A.=0.248  ^m,  and 
lfo-N,,o  =  7.5X\0‘°  cm“^  (a)  The  distribution  of  the  normal¬ 
ized  intensity,  (b)  The  distribution  of  the  normalized  electron 
density,  (c)  Radial  dependence  of  the  asymptotic  solutions  for 
the  normalized  amplitude  and  the  normalized  elec¬ 

tron  density  N,(r)/NQ  for  s  =0.513. 
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above,  namely  (1)  self-channeling,  (2)  stabilization  of  the 
mode  of  propagation,  (3)  the  confinement  of  a  substantial 
fraction  of  the  incident  power,  and  (4)  the  formation  of 
paraxial  cavitated  channels  in  the  electron  distribution 
are  common  aspects  of  the  dynamics  over  a  wide  range  of 
conditions. 

The  character  of  these  findings  is  now  illustrated  with 
five  representative  examples.  For  the  propagation  of  a 
beam  with  an  incident  plateaulike  transverse-intensity 
profile,  given  by  Eq.  (27)  with  iVj  =  8,  which  has  a  flat  ini¬ 
tial  wave  front  incident  on  an  initially  homogeneous  plas¬ 
ma  (the  values  of  the  parameters  have  been  adopted  from 
above  and  Pq/P^^  =  20.  168),  the  asymptotic  radial  profile 
of  the  intensity  distribution  contains  approximately  77% 
of  the  incident  power  [Fig.  7(a)].  The  corresponding 
asymptotic  amplitude  is  the  lowest  eigenmode  of  the  non¬ 
linear  Schrddinger  equation  with  s=:0.5I5.  The  propa¬ 
gation  of  the  same  wave  form  in  a  plasma  column,  with 
the  initial  electron  distribution  defined  by  Eq.  (29)  with 
A'^3=8  and  =rQ,  results  in  a  quasistabilized  intensity 
distribution  containing  34%  of  the  incident  power  [Fig. 
8(a)].  The  purely  relativistic  propagation  of  incident  pla¬ 
teaulike  pulses  [iV2  =  8  in  Eq.  (27)]  in  plasma  columns 
[A^3  =8  and  r,  —r^  in  Eq.  (29)]  also  results  in  the  forma¬ 
tion  of  quasistabilized  regimes  which,  in  this  case,  arise 
dynamically  from  the  defocusing  action  of  the  refraction 


generated  by  the  transverse  profiles  of  the  electron  densi¬ 
ty  [24].  For  this  situation,  approximately  25%  of  the  to¬ 
tal  incident  power  is  confined  [Fig.  5(d)].  The  compar¬ 
ison  of  these  cases  with  the  two  examples  discussed  above 
involving  the  charge-displacement  mechanism  indicates 
that  the  increase  in  the  value  of  the  confined  power  stems 
principally  from  the  substantial  additional  focusing  ac¬ 
tion  arising  from  the  inhomogeneous  electron  distribu¬ 
tion  produced  by  the  ponderomotive  force. 

Overall,  the  computations  reveal  that  the  charge  dis¬ 
placement,  which  generally  results  in  electronic  cavita¬ 
tion,  plays  an  important  role  in  stabilizing  the  mode  of 
confined  propagation  that  develops  dynamically.  This 
stabilization  naturally  occurs  by  refraction  of  the  radia- 


FIG.  8.  The  self-channeling  of  a  pulse  with  a  hyper-Gaussian 
incident-transverse-intensity  distribution  [iV2  =  8  in  Eq.  (27)] 
and  a  flat  initial  wave  front  in  a  preformed  column-shaped  plas¬ 
ma  [^3  =  8,  r*  =ro  in  Eq.  (29)]  for  the  case  of  relativistic  and 
charge-displacement  propagation  with  /o  =  3X10'’  W/cm,’ 
'’0  =  3  fim,  A.=0.248  fim,  and  Vo=V,.o  =  7.5X  10^°  cm“’.  (a) 
The  distribution  of  the  normalized  intensity,  (b)  The  distribu¬ 
tion  of  the  normalized  electron  density. 


PTG.  9.  Self-channeling  of  an  initially  focused  [i?^=/?y  g/2 
in  Eq.  (35)]  pulse  with  a  Gaussian  incident-transverse-intensity 
distribution  [^2  =2  in  Eq.  (35)]  in  initially  homogeneous  plasma 
for  the  case  of  relativistic  and  charge-displacement  propagation 
with  /o=3X10'’  W/cm^  ro  =  3  fim,  A.=0.248  /xm,  and 
^o~-^».o~7.5X  10“  cm“^.  (a)  The  distribution  of  the  normaJ- 
ized  intensity,  (b)  The  distribution  of  the  normalized  electron 
density,  (c)  Radial  dependence  of  the  asymptotic  solutions  for 
the  normalized  amplitude  [/,{r)//o]'^^  and  the  normalized  elec¬ 
tron  density  N;{r)/No  for  j  =0.566. 
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tion  into  the  central  paraxial  region.  Indeed,  the 
influence  of  the  charge  displacement  on  the  propagation 
is  so  strong  that  self-channeling  occurs  even  in  the  cases 
of  extremely  focused  and  defocused  incident  wave  forms. 
The  results  of  the  corresponding  calculations  for  ex¬ 
tremely  focused  incident  wave  forms,  for  R^=j?^o/2 
[Eq.  (35)],  with  both  Gaussian  [iV2=2  in  Eq.  (35)]  and 
plateaulike  in  Eq.  (35)]  incident-transverse- 

intensity  distributions,  are  depicted  in  Figs.  9  and  10,  re¬ 
spectively  (the  parameters  of  beam  and  plasma  being  the 
same  as  the  other  examples  above).  In  these  two  cases, 
the  asymptotic  transverse  profiles  of  the  amplitude  and 
the  electron  density  are  found  to  correspond  to  the  lowest 
eigenmodes  of  the  nonlinear  Schrodinger  equation  with 
values  of  the  parameter  5=0.566  (Fig.  9)  and  5=0.505 


(Fig.  10).  Note,  in  strong  contrast  to  the  situation  involv¬ 
ing  charge  displacement,  the  analogously  strongly  fo¬ 
cused  incident  wave  forms  for  the  purely  relativistic  case 
propagate  in  the  single-focus  regime  [Fig.  5(c)],  and  stable 
confinement  does  not  develop. 

Figure  1 1  illustrates  the  relativistic  and  charge- 
displacement  propagation  of  an  initially  strongly  de¬ 
focused  wave  form  [Rf  =  —Rf^Q/l  in  Eq.  (35)]  having  a 
Gaussian  initial  transverse-intensity  distribution  [A''2  =  2 
in  Eq.  (35)].  After  the  initial  stage  of  the  defocusing,  this 
pulse  evolves  into  a  paraxial  structure  that  is  analogous 
to  those  described  above.  Furthermore,  in  this  example 
the  asymptotic  transverse  profiles  of  the  amplitude  and 
the  electron  density  are  found  to  correspond  to  the  lowest 
eigenmode  of  the  nonlinear  Schrodinger  equation  with 
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FIG.  10.  The  self-channeling  of  an  initially  focused 
[Rf—Rfo/I  in  Eq.  (35)]  pulse  with  a  hyper-Gaussian  incident- 
transverse-intensity  distribution  [A^2=8  in  Eq.  (35)]  in  initially 
homogeneous  plasma  for  the  case  of  relativistic  and  charge- 
displacement  propagation  with  /o  =  3X10‘’  W/cm^  ro=3  ^im, 
A.=0.248  fim,  and  A^o  =  A^,.o  =  7.5X  1(F°  cm"^  (a)  The  distribu¬ 
tion  of  the  normalized  intensity,  (b)  The  distribution  of  the  nor¬ 
malized  electron  density,  (c)  Radial  dependence  of  the  asymp¬ 
totic  solutions  for  the  normalized  amplitude  [/,(/-)//o]'''^  and 
the  normalized  electron  density  N,{r)/No  for  5  =0.505. 


FIG.  11.  The  self-channeling  of  an  initially  defocused 
=  in  Eq.  (35)]  pulse  with  a  Gaussian  incident- 

transverse-intensity  distribution  (lVi  =  2  in  Eq.  (35)]  in  initially 
homogeneous  plasma  for  the  case  of  relativistic  and  charge- 
displacement  propagation  with  To  =  3X10”  W/cm^  ro  =  3  fim, 
^=0.248 /zm,  and  iV()=lV^()=7.5X  10^  cm  (a)  The  distribu¬ 
tion  of  the  normalized  intensity,  (b)  The  distribution  of  the  nor¬ 
malized  electron  density,  (c)  Radial  dependence  of  the  asymp¬ 
totic  solutions  for  the  normalized  amplitude  [/,(r)//o]''^  and 
the  normalized  electron  density  N,(r)/No  for  5  =0.800. 
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5=0.800. 

For  these  values  of  the  dimensionless  parameters  as¬ 
sumed  in  these  examples  (Figs.  9-11)  (a,  =248.6192, 
that  is  Pq/P^j  =  21.252  for  Gaussian  incident 
wave  form),  the  magnitude  of  the  quantity 
the  factor  determining  the  degree  of  the  initial  focusing 
(defocusing),  is  0.124.  For  1=0.248  /tm  and  ro  =  3  /xm, 
the  focusing  distance  R  =krlRj-Q/2,  used  in  the  exam¬ 
ples  above,  has  a  magnitude  of  28.35  ^m.  The  fact  that 
the  self-channeling  of  such  extremely  focused  or  de- 
focused  distributions  occurs  (Figs.  9-11)  means  the  self¬ 
channeling  should  be  expected  to  be  nearly  independent 
of  the  initial  conditions  of  focusing  or  defocusing. 

A  critical  power  for  self-channeling  arising  from 
the  relativistic  charge-displacement  mechanism  can  be 
defined.  We  now  discuss  this  issue  with  a  model 
developed  in  an  earlier  study.  The  threshold  power  P„  is 
defined  as  the  power  that  separates  the  asymptotic  behav¬ 
ior,  with  respect  to  large  distance  of  propagation  iz),  into 
two  distinct  classes.  Since  its  definition  rests  on  an 
asymptotic  property,  it  can  only  exist  in  the  limiting  case 
of  vanishing  losses  (^=0  in  Ref.  [24]).  For  a  power 
^  <Fcr,  the  asymptotic  transverse-intensity  profile  tends 
to  zero  at  large  z.  In  contrast,  for  P^P„  the  asymptotic 
profile  of  the  intensity  tends  to  the  lowest  eigenmode  of 
the  governing  nonlinear  Schrodinger  equation.  Addition¬ 


ally,  the  normalized  value  of  this  critical  power  is  exactly 
twice  the  corresponding  normalized  value  of  the  critical 
power  for  the  Kerr  effect  self-focusing  found  in  cubic 
media  for  initially  flat  wave  forms  [34,37].  The  expres¬ 
sion  corresponding  to  this  critical  power  [36]  for  the  rela¬ 
tivistic  and  charge-displacement  mechanism  is 

/e'^)  fJglip)pdp{co/apQ)^ 

=  1.6198X10‘°(tu/ta^,o)^  W  ,  (36) 

where  go(p)  is  the  Townes  mode  [34]. 

It  has  also  been  shown  that  the  value  of  the  relativistic 
and  charge-displacement  self-channeling  threshold 
power,  in  cases  involving  initially  focused  or  defocused 
beams,  exceeds  P„  and  depends  on  the  degree  of  the  ini¬ 
tial  focusing  or  defocusing.  For  a  given  magnitude  of  the 
curvature  of  the  wave  front,  the  value  of  the  threshold 
power  for  initially  defocused  beams  is  greater  than  that 
for  initially  focused  beams.  Finally,  it  has  been  shown 
both  analytically  and  numerically  that  for  self-channeling 
of  an  arbitrary  wave  form  to  occur,  it  is  sufficient  that  the 
Hamiltonian  of  the  purely  relativistic  case,  considered  as 
a  functional  of  the  initial  transverse  amplitude  distribu¬ 
tion,. should  be  negative  [36].  The  precise  statement  for 
this  condition  is 


3 


duQ 

dr 


-a,(|«ol^-{2/a2)[(l+a2|«ol^)‘'^-l]) 


r  dr  <0  . 


(37) 


We  now  comment  on  the  behavior  of  pulses  having  ini¬ 
tial  amplitude  distributions  Uoir)  close  to  higher  eigen- 
modes  Fj  n(r),  n  >  1.  These  higher  modes  generally  are 
associated  with  electronically  cavitated  channels.  The 
self-channeling  in  these  cases  could  result  in  asymptotic 
distributions  of  intensity  and  electron  density  corre¬ 
sponding  to  certain  higher  eigenfunctions 
(Vj  n(r),  /!  >  1.  lu  particular,  this  conjec¬ 
ture  has  established  an  initial  amplitude  distribution 
UQ(r),  which  is  close  to  the  first  eigenmode  F,  ,(r)  with 
s  =0.544.  Direct  numerical  calculations  showed  that  the 
asymptotic  distributions  of  both  intensity  and  electron 
density  that  evolved  in  this  case  correspond  to  the  first 
eigenmode  I,  ,lr)  =  V,^Jr),  with  s  essentially  equal 

to  0.544.  It  should  be  additionally  noted,  however,  that 
these  higher  eigenmodes  ( «  >  1 )  are  quite  possibly  unsta¬ 
ble  against  small  azimuthal  perturbations  that  destroy 
their  axial  symmetry. 

The  principal  result  of  this  section  is  the  finding  that 
the  combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms  can  result  in  self-channeling 
with  the  formation  of  stabilized  paraxial  modes  over  a 
rather  wide  range  of  physical  conditions.  Moreover, 
these  spatially  confined  modes  are  generally  associated 
with  corresponding  cavitated  channels  in  the  electron 
density.  Finally,  the  characteristics  of  these  channeled 


modes  have  asymptotic  behavior  that  is  described  by  the 
appropriate  lowest  eigenmodes  of  the  governing  non¬ 
linear  Schrodinger  equation. 

V.  CONCLUSIONS 

A  theoretical  approach  suitable  for  the  numerical  in¬ 
vestigation  of  the  two-dimensional  {r,z)  dynamics  of 
propagation  of  coherent  ultrashort  (r,  »t»t^)  relativ¬ 
istic  laser  pulses  in  cold  underdense  plasmas  has  been 
developed.  Four  basic  physical  phenomena  are  included 
within  the  scope  of  this  method.  They  are  (1)  the  non¬ 
linear  dependence  of  the  index  of  refraction  due  to  the 
relativistic  increase  in  the  mass  of  the  electrons,  (2)  the 
variation  of  the  index  of  refraction  resulting  from  the 
perturbation  of  the  electron  density  by  the  ponderomo- 
tive  force,  (3)  the  diffraction  of  the  radiation,  and  (4)  the 
refraction  caused  by  nascent  transversely  inhomogeneous 
electron-density  distributions.  The  equations  studied  in 
this  work  may  be  regarded  as  the  generalization  of  those 
treated  previously  in  other  studies  [18,19,33]  involving  in¬ 
itially  inhomogeneous  plasmas.  Further  studies  are  con¬ 
tinuing  with  an  extension  of  this  analysis,  which  includes 
the  azimuthal  coordinate  in  the  description  of  the  propa¬ 
gation. 
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The  main  conclusions  stemming  from  the  calculations 
are  the  following. 

(1)  The  cooperative  effect  of  the  relativistic  and 
charge-displacement  mechanisms  leads  asymptotically  to 
stable  high-intensity  z-independent  modes  of  self¬ 
channeling,  and  a  major  fraction  of  the  incident  power 
can  be  confined  in  these  paraxial  modes.  Stable  cavita¬ 
tion  of  the  electron  density  is  a  general  feature  of  these 
spatially  confined  modes. 

(2)  The  z-independent  modes,  to  which  the  solutions  of 
the  equation  describing  the  relativistic  and  charge- 
displacement  propagation  tend  asymptotically,  are  recog¬ 
nized  as  the  lowest  eigenmodes  of  the  governing  nonlinear 
Schrddinger  equation. 

(3)  A  separate  study  of  purely  relativistic  propagation 
shows  that  beams  with  flat  incident  phase  fronts  exhibit  a 
pulsing  behavior  in  homogeneous  plasmas  but  can  under¬ 
go  quasistabilization  in  suitably  configured  plasma 
columns.  However,  sufficiently  sharply  initially  focused 


beams  generally  exhibit  the  development  of  only  a  single 
focus.  In  significant  contrast,  the  present  study  shows 
that  with  both  the  relativistic  and  charge-displacement 
mechanisms,  initially  focused  beams  also  generally  lead 
to  confined  modes  of  propagation.  Finally,  it  should  be 
noted  that  the  equation  for  relativistic  self-focusing  can 
be  considered  as  a  general  model  equation  of  propagation 
in  saturable  nonlinear  media. 
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Experimental  studies  examining  a  new  relativistic  regime  of  high-intensity  short-pulse  propagation  in 
plasmas  have  been  performed  which  present  evidence  for  the  formation  of  a  stable  mode  of  spatially 
confined  (channeled)  propagation.  For  an  electron  density  of  -^1.35x10^'  cm”^  and  a  power  of 
—  3x  lo”  W,  the  results  indicate  a  channel  radius  <  I  pm  and  a  peak  intensity  — lO'*  W/cm^.  Com¬ 


parison  of  these  findings  with  a  dynamical  theory 
and  the  radial  extent  of  the  propagation  observed. 

PACS  numbers:  52.40.  Db,  4125.-p,  42.65. Jx 

A  fundamentally  new  regime  of  electromagnetic  propa¬ 
gation  is  expected  to  arise  in  plasmas  for  short-pulse  radi¬ 
ation  at  sufficiently  high  intensity.  Calculations  of  the 
propagation  in  plasmas,  incorporating  both  relativistic 
[1,2]  and  charge-displacement  mechanisms  [3-7],  indi¬ 
cate  that  the  combined  action  of  these  processes  can  lead 
to  a  new  stable  form  of  spatially  confined  (channeled) 
propagation.  This  Letter  (I)  reports  the  results  of  the 
first  experimental  study  probing  the  physical  regime 
relevant  to  the  observation  of  relativistic  and  charge- 
displacement  self-channeling  and  (2)  presents  the  initial 
comparison  of  these  experimental  findings  with  matching 
theoretical  calculations  performed  with  the  computation- 
‘  al  procedures  described  in  Ref.  [7], 

The  experimental  arrangement  used  in  these  studies  is 
illustrated  in  Fig.  1(a).  The  source  of  radiation  was  a 
subpicosecond  KrF*  (A.  “248  nm)  laser  that  has  been  de¬ 
scribed  elsewhere  [8],  It  delivered  a  linearly  polarized 
power  of  —  3x  10"  W  (  —  150  mJ,  pulse  duration  —500 
fs)  in  a  beam  with  a  diameter  of  —42  mm.  When  this 
radiation  was  focused  into  the  chamber  with  lens  L I 
(//7),  a  focal  radius  /-q- 3.5  pm  was  measured,  giving  a 
maximum  intensity  /q- 8.6x|0"  W/cm^.  The  medium 
was  provided  by  filling  the  chamber  statically  with  gas 
[He,  Ne,  Ar,  Kr,  Xe,  N2,  CO2,  or  a  mixture  of  Xe  (4%) 
and  N2  (96%)]  up  to  a  maximum  density  of  — 1.89  x  I0“ 
cm 

The  diffracted  248-nm  radiation  was  measured  as  a 
function  of  the  angle  (©)  with  respect  to  the  direction  of 
the  incident  radiation.  The  incident  laser  beam  was 
blocked  by  a  metal  disk  on  the  output  window  of  le 
chamber  and  lens  LI  imaged  the  region  near  the  focal 
zone  on  a  fluorescent  screen  S.  The  diaphragm  D  in  front 
of  lens  L2  restricted  the  collection  of  the  diffracted  light 
to  a  solid  angle  of  —5®  while  simultaneously  increasing 
the  depth  of  field.  The  angle  between  the  axis  of  the  lens 
L  2  and  the  axis  of  the  incoming  laser  radiation  could  be 
readily  varied  up  to  a  maximum  angle  of  0—15®.  Two 
flat  mirrors  coated  for  high  reflection  (—99%)  at  248 
nm,  both  having  a  spectral  bandwidth  of  — 10  nm,  served 


yields  agreement  for  both  the  longitudinal  structure 


in  reflection  as  spectral  filters  (F)  for  the  diffracted  laser 
radiation  so  that  only  the  scattered  248-nm  radiation 
could  illuminate  the  screen.  An  attenuator  A  was  em¬ 
ployed  to  adjust  the  intensity  on  the  screen  and  the  im¬ 
ages  formed  were  recorded  through  the  visible  fluores- 


Axial  Distance  z  (/im) 

FIG.  I.  (a)  Experimental  apparatus  used  in  studies  of  propa¬ 
gation.  See  text  for  description,  (b)  Data  concerning  the  pat¬ 
tern  of  propagation  observed  with  a  single  pulse  in  N2  at  a  den¬ 
sity  of  — 1.35  X 10“  cm~\  The  maximum  intensity  is  half  the 
detector  (CCD)  saturation.  The  radiation  is  incident  from  the 
left.  Inset:  Photographic  data  with  a  vertical  sj^tiai  resolution 
of  —10  ^m.  The  graph  illustrates  the  one-dimensional  axial 
profile  taken  along  the  direction  of  propagation  (z)  of  the  pho¬ 
tographic  data  (inset).  The  spacing  of  the  maxima,  £^200 
±  20  ^m,  is  indicated. 
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ccnce  produced  with  a  microscope  and  a  charge-coupled- 
device  (CCD)  camera. 

The  characteristic  behavior  observed  is  well  illustrated 
by  the  data  recorded  with  N2.  The  measured  result, 
shown  in  Fig.  1(b).  corresponds  to  a  density  pnj--L35 
X  I0^°  cm  To  the  left,  in  the  photographic  inset,  a  rel¬ 
atively  large  cone  of  light  Rayleigh  scattered  from  the 
plasma  is  visible  at  all  angles  as  the  energy  propagates  to¬ 
ward  the  focal  point  of  the  lens,  while,  in  the  region  to  the 
right  of  the  conical  apex,  a  narrow  filament  developed. 
The  diameter  of  this  filament  is  not  greater  than  10  pm, 
the  measured  spatial  resolution  of  the  imaging  system. 
The  distribution  of  intensity  observed  along  the  filament 
exhibited  several  bright  features  attributed  to  diffraction 
because  they  could  not  be  seen  for  0>2O°.  Since  the 
axis  of  the  imaging  lens  corresponded  to  an  angle 
6 "7.5°,  the  scale  along  the  abscissa  of  the  photographic 
data  is  reduced  by  almost  eightfold,  giving  the  maximum 
length  of  the  filament  as  ~  I  mm.  The  graph  in  Fig.  1(b) 
represents  a  one-dimensional  axial  profile,  taken  along 
the  direction  of  propagation  (z),  of  the  observed  intensity 
pattern  (inset).  Three  peaks  y)  are  visible  with  a  a 
spatial  separation  of  5 -200  ±20  pm.  The  normal  Ray¬ 
leigh  range  for  the  focal  geometry  used  was  —200  pm. 

The  diameter  of  the  filamentary  channel  is  an  impor¬ 
tant  dynamical  variable,  which  we  estimated  by  measur¬ 
ing  the  maximum  angular  deviation  of  the  diffracted 
light.  The  experimental  value  0  of  this  diffracted  cone 
was  0—20°,  a  magnitude  indicating  a  radius  /■*— 0.9  pm 
though  the  relation  0  — I.22X//’#.  Filaments  of  this  gen¬ 
eral  nature  were  observed  at  densities  above  — 1.35  x  |0^® 
cm  in  N2,  Nc,  Ar,  Kr,  CO2,  and  a  mixture  of  Xe  (4%) 
and  N2  (96%),  but  not  in  He  and  Xe,  two  materials  dis¬ 
cussed  further  below. 

Two  mechanisms  exist  that  could  modify  the  refractive 
index  of  the  medium  and  lead  to  the  observed  behavior. 
They  are  (1)  the  Kerr  effect  stemming  from  the  ions  and 
(2)  the  relativistic  and  charge-displacement  process  (7l. 
Since  the  pulse  duration  is  very  short  (  —  500  fs),  the 
motion  of  the  ions  is  negligible  [3],  and  no  contribution 
can  arise  from  expulsion  of  the  plasma  from  the  high- 
intensity  zone.  An  implication  of  the  estimate  of  the 
channel  radius  (r^— 0.9  pm)  is  that  the  observed  propa¬ 
gation  is  associated  with  intensities  in  the  I0'*-I0” 
W/cm^  range.  Under  these  conditions,  available  experi¬ 
mental  evidence  [9,10]  on  multiphoton  ionization  indi¬ 
cates  that  He  should  be  fully  ionized  and  the  C,  N,  and  O 
atoms  constituting  the  molecular  materials  would  retain, 
at  most,  only  I  j  electrons. 

Consider  explicitly  the  case  of  N2,  which  has  estimated 
[9,10]  threshold  intensities  for  the  production  of 
N*-^,  and  of  l.6x|0“,  6.4xio'*.  and  1.3x10” 
W/cm^  respectively.  Hence,  the  volume  of  the  channel 
would  be  largely  ionized  to  with  certain  localized 
high-intensity  regions  contributing  some  N*"*".  Two 
consequences  of  this  pattern  of  ionization  follow,  namely, 

(i)  the  Kerr  effect  arising  from  the  ions  is  small,  since  the 

2310 


polarizabilities  of  the  remaining  lj  electrons  are  low,  and 

(11)  the  electron  density  (n,)  initially  produced  in  the  fo¬ 
cal  region  is  nearly  uniform.  Therefore,  1.35x10^' 

cm  for  the  data  on  N2  shown  in  Fig.  I  (b). 

A  critical  power  P^,  for  seif-channeling,  arising  from 
the  relativistic  and  charge-displacement  mechanism,  can 
be  defined  [ill  as 

P^t  - im}c  Ve ■)  gHr)r 

=  i.62x|0''’(/icr/n,)  W,  (I) 

where  nif,  e,  and  c  have  their  customary  identifications,  (u 
is  the  laser  angular  frequency,  is  the  plasma  frequen¬ 
cy  for  the  uniform  unperturbed  plasma  with  electron  den¬ 
sity  n„  n„  is  the  critical  electron  density  (for  A.  — 248  nm, 
ncr"L82x  10^^  cm  ^),  and  gair)  is  the  Townes  mode 

[12] . 

The  critical  powers  associated  with  the  experimental 
conditions,  for  He  and  N2  at  a  medium  density 
p-l.35xlO^“  cm■^  arc  1.08x10'^  and  2,19x10"  W, 
respectively.  Therefore,  since  the  incident  power  was 
/*==3x  10"  W,  no  filament  was  expected  in  He,  a  predic¬ 
tion  conforming  with  the  observation  of  none.  Moreover, 
the  diffracted  cone  of  radiation  was  also  absent  with  He. 
In  contrast,  /»//>„=»  1.37  for  N2,  a  condition  that  held 
generally  {P/P„  >  I )  for  all  materials  which  exhibited 
evidence  for  channel  formation.  We  note,  however,  that 
some  contribution  from  the  Kerr  effect  may  be  present, 
even  for  the  light  materials  (Ne,  N2,  and  CO2),  in  the 
early  stage  of  channel  formation  prior  to  the  development 
of  a  substantial  level  of  ionization,  and  that  the  heavier 
gases  (Ar,  Kr,  and  Xe)  may  involve  a  more  significant 
influence  from  the  Kerr  process.  A  specific  estimate  of 
the  nonlinear  index  change  arising  from  both  and 

at  an  intensity  of  —  lO”  W/cm^  indicates  that  their 
contribution  is  less  than  10"^  that  of  the  free  electrons: 
hence  the  ionic  contribution  can  be  neglected  in  N2  for 
the  conditions  studied. 

A  direct  comparison  will  now  be  made  between  the 
theoretical  analysis,  fully  described  in  Ref.  [7],  and  the 
experimental  findings  for  N2.  This  comparison  can  be  ac¬ 
complished  for  both  the  longitudinal  intensity  profile  and 
the  radial  extent  of  the  channel.  Figure  2(a)  illustrates 
the  intensity  profile  /(r,z)//o  calculated  with  physical  pa¬ 
rameters  corresponding  to  those  of  the  experiment  (i.e., 
/>=3xio"  W,  ro-3.5  pm,  /i,  -  1.35x  |0^'  cm■^  and 

cf"l'37).  Importantly,  all  of  these  parameters  are 
based  on  independent  measurements  of  (1)  the  laser  pulse 
involving  determinations  of  the  energy  and  power  (P), 
(2)  the  focal  radius  (ro)  of  the  incident  radiation,  and  (3) 
the  characteristics  of  the  multiphoton  ionization  [10]  gen¬ 
erating  the  electron  density  («,).  Therefore,  this  compar¬ 
ison  docs  not  involve  a  fit  with  a  free  parameter.  The 
normalized  electron  density  calculated  is  presented  in 
Fig.  2(b),  from  which  it  is  seen  that  electronic  cavitation 
occurs  only  near  the  positions  of  the  maxima  in  the  inten¬ 
sity  profile  [Fig.  2(a)].  Curve  A  in  Fig.  2(c)  represents 
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the  one-dimensional  axial  intensity  profile  /(0,z)//o  cor¬ 
responding  to  the  calculated  distribution  shown  in  Fig. 
2(a).  The  spacing  (S)  of  the  maxima  is  seen  to  be 
(5~I85  fjm,  a  value  in  close  agreement  with  the  experi¬ 
mental  figure  (S~200  ±20  /jm)  illustrated  in  Fig.  1(b). 
Furthermore,  analysis  has  shown  that  the  spacing  <5  is 
quite  sensitive  to  the  power  Fq  and  electron  density  n,, 
particularly  in  the  region  close  to  the  threshold  (see  Eq. 
(I)].  With  respect  to  the  results  illustrated  in  Fig.  2(c), 
an  increase  in  by  less  than  5%  causes  a  reduction  in 
the  spacing  S  by  approximately  25%.  Therefore,  substan¬ 
tial  changes  in  the  physical  parameters  would  grossly 
alter  the  comparison  of  the  experimental  and  theoretical 
results. 

Theoretical  studies  [3-7,11]  indicate  that  the  charge 
displacement  plays  a  very  imfx)rtant  dynamical  role.  In 
order  to  test  this  hypothesis,  identical  calculations  were 
made  for  N2  for  the  purely  relativistic  case  [2]  which  ex¬ 
plicitly  neglects  the  charge-displacement  term,  namely, 
elimination  of  the  term  +///,) '^^  in  Eq. 

(24)  of  Ref.  (7l.  Significantly,  the  resulting  axial  profile 
(curve  B  in  Fig.  2(c)]  exhibits  only  a  stng/e  relatively 
weak  maximum,  for  0  <  r  <  600  /tm.  an  outcome  sharply 
at  variance  with  both  the  full  theoretical  analysis  and  the 
experimental  observation.  Although  the  expected  charge 
displacement  is  highly  localized  (Fig.  2(b)],  this  compar¬ 
ison  reveals  the  strong  influence  it  has  on  the  propagation 
[3,7,11].  At  a  greater  incident  power  (F/Fgf-^IO),  a 
continuous  channel  in  the  electron  distribution  is  expect¬ 
ed  to  develop  [7]. 

The  measurements  indicated  an  approximate  value  of 
r^—0.9  pm  for  the  radial  extent  of  the  channel,  a  result 
that  can  be  compared  with  the  corresponding  theoretical 
figure.  Figure  2(d)  illustrates  five  radial  intensity  profiles 
/(r,r,  )//o  of  the  distribution  pictured  in  Fig.  2(a).  Since 
the  measurement  of  this  angularly  scattered  radiation  did 
not  correspond  to  a  known  longitudinal  position,  this 
comparison  can  only  be  qualitative,  but  the  radial  distri¬ 
butions  shown  indicate  that  the  expected  value  lies  in  the 
interval  0.5<r<  1.0  pm,  a  range  that  comfortably  in¬ 
cludes  the  experimental  value 
The  results  observed  with  Xe  deserve  additional  discus¬ 
sion,  since  those  experiments  did  not  give  evidence  for  the 
formation  of  a  channel.  In  significant  contrast  to  the  case 
involving  N:,  the  electron  density  produced  by  the 
multiphoton  ionization  1 10]  in  Xe  is  expected  to  be  very 

FIG.  2.  Culcuialions  for  Nj  with  10"  W,  ro-J.S  pm, 

-1.35x10-'  cm-\  and  /o-8.6x|0'^  W/cm^  (a)  Normal¬ 
ized  intensity  /{r,z)/lo.  (b)  Normalized  electron  density 
/y(r.r)/yVo  for  Nj  with  (c)  Normalized  one-dimen¬ 

sional  axial  intensity  profiles  /(O.z)//®.  Curve  /}.  full  theory  for 
data  in  panel  (a).  5—185  pm.  Curve  B,  calculation  with 
charge-displacement  term  neglected,  (d)  Normalized  radial  in¬ 
tensity  profiles  /(r,Zi)//o  corresponding  to  panel  (a).  Longitu¬ 
dinal  positions  n  -172  pm,  r2-245  pm,  r.,-358  pm,  ;4-44l 
pm.  and  zj  —559  pm  and  r^— 0.9  pm. 
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nonuniform  spatially.  For  intensities  spanning  lO'^-IO'* 
W/cm^,  the  corresponding  density  n,  would  vary  by  over 
a  factor  of  2.  Since  this  nonuniformity  would  tend  to 
reduce  the  refractive  index  locally  in  the  central  high- 
intensity  region,  a  significant  defocusing  action  is  expect¬ 
ed  which  could  suppress  the  channel  formation. 

Finally,  we  note  (1)  that  the  intensity  distribution  is 
not  expected  to  depend  strongly  on  the  state  of  polariza¬ 
tion  [13,14j  and  (2)  that  losses  to  the  plasma  may  be 
significant,  particularly  at  electron  densities  close  to  /!„• 

The  first  experiments  examining  a  new  relativistic  re- 
_  gime  of  high-intensity  pulse  propagation  in  plasmas  have 
been  performed  and  the  findings  indicate  the  formation  of 
a  channeled  mode  of  propagation  over  a  length  consider¬ 
ably  greater  than  the  Rayleigh  range.  Specific  compar¬ 
isons  of  the  experimental  observations  with  a  dynamical 
theory,  which  explicitly  includes  both  the  influence  of  the 
relativistic  mass  shift  and  the  displacement  of  the  elec¬ 
tronic  component  of  the  plasma,  produce  excellent  agree¬ 
ment  for  both  the  longitudinal  structure  of  the  intensity 
profile  and  the  radial  extent  of  the  channel.  While  the 
present  channel  contains  several  foci,  a  continuous  chan¬ 
nel  is  predicted  to  develop  at  higher  power.  Finally,  the 
intrinsically  very  high  concentration  of  power  associated 
with  this  mechanism  of  channeled  propagation  provides 
an  efficient  and  general  method  for  the  production  of  con¬ 
ditions  necessary  for  x-ray  amplification  tl5j. 
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ABSTRACT 


The  ability  to  apply  power  densities  controllably  at  or  above  vigorous  thermonuclear 
levels  (>  10^^  W/cm^)  in  materials  is  the  basic  issue  for  achieving  efficient  amplification  of 
X-rays.  Recent  experimental  and  theoretical  findings  concerning  (1)  the  multiphoton  production 
of  X-rays  from  clusters  and  (2)  high-intensity  modes  of  channeled  propagation  in  plasmas 
indicate  an  entirely  new  method  for  producing  the  conditions  necessary  for  strong  amplification 
in  the  multi-kilovolt  range.  These  two  new  nonlinear  phenomena  are  being  united  to  produce 
and  control  the  imperative  power  compression;  the  multiphoton  mechanism  serves  to  establish 
the  condition  locally  while  the  confined  propagation  provides  the  required  spatial  organization. 
The  present  work,  which  experimentally  demonstrates  the  first  combined  expression  of  these  two 
complex  nonlinear  processes  through  direct  X-ray  imaging  of  Xe(M)  emission  (~  1  keV)  in 
stable  self-trapped  channels,  (a)  reveals  the  exceptional  compatibility  of  their  mutual  scaling  for 
realizing  the  necessary  power  density,  (/3)  provides  confirming  evidence  for  the  action  of  a 
superstrong  coherent  multi-electron  intense-field  interaction  in  the  X-ray  generation  from  the  Xe 
clusters,  and  (7)  furnishes  new  detailed  information  on  the  dynamics  of  the  radial  intensity 
distributions  associated  with  the  channeled  propagation.  The  resulting  knowledge  of  the  scaling 
relations  underlying  these  phenomena  enables  the  optimum  conditions  for  amplification  to  be 
specified  up  to  a  quantum  energy  of  ~  5  keV.  The  harmonious  use  of  these  new  nonlinear 
processes  is  expected  to  lead  to  an  advanced  generation  of  extraordinarily  bright  X-ray  sources 
in  the  multi-kilovolt  region  having  a  peak  brightness  of  ~  lCp^-10^^  7S“^  (mrad)"^  (mm)“^ 
(0. 1  %  BW)"^,  a  level  sufficient  for  biological  holographic  imaging  capable  of  providing  a  high 
resolution  visualization  of  the  molecular  anatomy  of  cells,  tissues,  and  organisms  in  the  natural 
state. 
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Stable  channeled  propagation  of  subpicosecond  (248  nm)  radiation  arising  from 
relativistic-charge-displacement  nonlinearities  [1]  has  been  observed  [2,3]  in  underdense  plasmas. 
Experimental  tests  [4]  of  a  concept  [5]  for  a  molecular  design  associated  with  multiphoton  X-ray 
production  induced  by  a  strong  nonlinear  interaction  have  shown  that  copious  Xe(M),  Kr(L),  and 
Xe(L)  emission  can  be  generated  with  corresponding  Kr  and  Xe  clusters  in  the  1-5  keV  region 
[4,6-8].  Estimates  stemming  from  the  observed  X-ray  spectra  indicated  the  achievement  of 
maximal  atom-specific  energy  transfer  rates  of  —  0. 1-1.0  W/atom. 

The  present  work  (1)  demonstrates  a  compatible  unification  of  the  physical  conditions 
governing  the  combined  expression  of  these  two  highly  nonlinear  phenomena,  (2)  provides  data 
corroborating  the  presence  of  a  superstrong  multiphoton  interaction  in  the  X-ray  production  from 
the  clusters  [6-9],  and  (3)  gives  important  new  information  on  the  properties  of  the  radial 
intensity  profiles  produced  by  the  nonlinear  propagation.  The  outcome  is  a  congruous  picture 
of  the  mutual  scaling  of  these  processes  that  enables  a  general  specification  of  the  optimum 
conditions  necessary  for  amplification  up  to  a  quantum  energy  of  ~  5  keV. 

Calculations  have  shown  that  the  observed  characteristics  of  the  channeled  propagation 
[1]  can  be  put  in  close  correspondence  [2,3]  with  a  theoretical  picture  incorporating  relativistic 
[10,11]  and  charge-displacement  [12-17]  mechanisms.  In  particular,  those  studies  indicated  that 
intensities  as  high  as  ~  10^®  W/cm^  were  produced  in  the  channels  [2,3],  a  value  sufficient  for 
the  copious  multiphoton  generation  of  multi-kilovolt  X-ray  emission  from  Kr  and  Xe  clusters  [4- 
8].  Therefore,  the  use  of  a  medium  composed  of  Kr  or  Xe  clusters  presents  the  opportunity  for 
a  study  combining  the  dynamics  of  the  confined  propagation  and  the  properties  of  the  nonlinear 
cluster  excitation  through  direct  X-ray  imaging  of  the  high-intensity  zones  produced  by  the 
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channeling.  We  describe  below  the  results  obtained  by  using  Xe(M)  emission  (-  1  keV) 
generated  from  Xe  clusters  in  channels  produced  with  subpicosecond  248  nm  radiation. 

The  experimental  arrangement  used  to  conduct  spectroscopic  studies  of  rare  gas  clusters 
has  been  described  elsewhere  [6].  The  laser  [18]  used  had  a  wavelength  of  248  nm,  a  pulse 
width  of  ~  300  fs,  a  power  of  ~  0.7  TW,  and  was  focused  to  a  maximum  intensity  of  —  10^^ 
W/cm-^  with  an  f/3  off-axis  parabolic  mirror  in  a  gas  target.  This  target  was  hydrodynamically 
produced  by  a  cooled  high-pressure  pulsed-valve  fitted  with  a  circular  sonic  nozzle  having  a 
diameter  of  ~  0.5  mm.  In  order  to  eliminate  distortion  and  absorption  of  the  incident  beam  as 
it  approached  the  focal  zone  through  the  outer  regions  of  the  material  produced  by  the  nozzle, 
a  barrier  with  a  ~  100  ^m  hole  was  positioned  near  the  opening  of  the  nozzle  as  shown  in  the 
inset  in  Fig.  2(b).  The  focus  of  the  f/3  mirror  was  then  placed  at  the  position  of  the  opening 
in  the  barrier  so  that  the  initial  conditions  of  medium  density  and  laser  intensity  corresponding 
to  the  launching  of  the  wave  were  well  defined. 

Previous  spectroscopic  studies  [6]  have  demonstrated  that  bright  Xe(M)  emission 
spanning  the  9.5-16  A  region  is  readily  produced  from  Xe  clusters.  In  order  to  record  the 
spatial  pattern  of  the  propagation,  this  emission  was  observed  with  an  X-ray  "pinhole"  (CCD) 
camera.  The  camera  was  an  Electrum  EDC1(X)0  using  a  video  chip  having  a  196  x  165  pixel 
array  with  an  individual  pixel  size  of  13.5  /xm  x  16  /xm.  The  X-ray  images  were  obtained  by 
viewing  the  region  excited  by  the  focused  248  nm  radiation  directly  below  the  pulsed-jet  in  a 
direction  transverse  to  the  propagation  of  the  248  nm  beam.  For  images  of  the  Xe(M)  emission, 
the  camera  was  equips  with  a  filter  with  the  composition  (2000  A  A1  -f  1  /ixm  polycarbonate 
+  10  /xm  Be)  and,  consequently,  detected  all  radiation  with  a  quantum  energy  ^  700  eV. 
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When  the  camera  was  equipped  with  a  ~  25  pinhole  and  positioned  to  give  a  magnification 
factor  of  1.18,  the  brightness  of  the  Xe(M)  emission  was  sufficiently  high  that  recordings  of 
single-shot  exposures  could  be  obtained.  Given  the  magnification  used  and  the  geometry  of  the 
pixel,  the  pinhole  camera  has  a  predicted  limiting  spatial  resolution  A  -  30  ^m,  a  property 
verified  by  the  observed  characteristics  of  the  images  described  below. 

Examples  of  the  X-ray  images  recorded  for  two  different  densities  of  the  Xe  medium  are 
illustrated  in  Figs.  1(a)  and  2(a),  respectively.  The  former  shows  a  strongly  longitudinally 
modulated  distribution  while  the  latter  demonstrates  that  a  highly  uniform  distribution  can  be 
produced  under  appropriate  conditions.  The  corresponding  axial  X-ray  intensity  distributions 
(radius  r  =  0  line-outs  of  these  data),  which  are  represented  in  Figs.  1(b)  and, 2(b),  are  found 
to  be  in  good  correspondence  with  the  matching  theoretical  distributions  shown  in  Figs.  1(c)  and 
2(c).  Significantly,  as  in  the  earlier  study  [2],  all  the  parameters  involved  in  the  computed 
distributions  are  based  on  quantities,  derived  from  independent  measurements,  so  that  the 
comparisons  with  the  experimental  results  are  sharply  defined  and  do  not  involve  a  fit  with  any 
free  parameters.  Importantly,  the  conditions  of  these  experiments  also  corresponded  to  the  zone 
of  stable  confined  propagation  previously  found  in  an  analysis  of  the  response  of  the  channeled 
power  distribution  to  small  azimuthal  perturbations  [16],  a  point  shown  explicitly  below. 
Furthermore,  all  the  computed  results  illustrated  in  Figs.  1  and  2  allow  for  the  variation  of  the 
density  of  the  medium  caused  by  the  hydrodynamic  expansion  from  the  nozzle  [19]. 

Of  chief  interest  in  the  present  context  is  the  radial  dependence  of  the  X-ray  images, 
since  the  behavior  of  the  transverse  dimension  of  the  emitting  region  represents  a  sensitive  probe 
of  both  the  dynamics  of  the  propagation  and  the  mechanism  of  excited  state  production  in  the 
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clusters.  Consider  the  image  shown  in  Fig.  2(a),  which  has  a  transverse  width  of  ~  120  /xm 
at  its  midpoint  (z  =  1050  /xm).  This  width  is  considerably  greater  than  the  resolution  (A  ~  30 
/xm)  of  the  camera  and  both  the  observed  and  calculated  values  [  -  2  ^m]  associated  with  the 
confined  propagation  of  the  248  nm  radiation  [2,3].  Consequently,  the  measurement  shows 
unambiguously  that  X-rays  are  being  produced  well  outside  the  small  central  core  of  the  channel 
that  confines  the  propagation  of  the  248  nm  energy.  In  addition,  the  image  shown  in  Fig.  1(a) 
clearly  demonstrates  that  the  radial  extent  of  the  X-ray  emission  can  have  a  significantly 
structured  longitudinal  (z)  dependence. 

The  existence  of  a  bright  emitting  zone  having  a  diameter  of  ~  100  /xm  and  the  observed 
longitudinal  signature  in  the  radial  structure  is  inconsistent  with  either  (a)  collisional  excitation 
arising  from  energetic  electrons  produced  in  the  core  of  the  channel  or  (b)  recombination  into 
excited  states  of  highly  ionized  ions  which  migrate  radially  from  the  central  region.  At  the 
outset  we  see  that  the  largest  observed  radial  dimension  in  Fig.  2(a)  is  greater  than  the 
corresponding  maximum  value  in  Fig.  1(a),  although  the  density  of  the  medium  in  the  former 
is  significantly  higher  than  in  the  latter.  This  trend  is  exactly  the  opposite  of  that  expected  for 
the  range  of  energetic  charged  particles,  either  electrons  or  ions,  in  matter.  Furthermore,  if  the 
origin  of  the  X-ray  emission  were  to  arise  from  collisional  excitation  by  electrons  expelled  from 
the  channel  through  the  effect  of  the  ponderomotive  potential,  then  it  should  be  possible  to  relate 
the  observed  radial  size  of  the  emitting  region  to  an  estimate  of  the  maximum  range  of  electrons 
[20]  in  materials  with  the  known  range-energy  [21]  relationships.  Specifically,  since  the  data  in 
Figs.  1(c)  and  2(c)  show  that  the  peak  intensity  in  the  channel  is  expected  to  fall  in  the  10^^- 
1020  wx/cm^  region,  for  electrons  having  kinetic  energies  on  the  scale  of  the  ponderomotive 
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potential  (—  100-500  keV),  the  expected  range  would  be  —  100  mg/cm^.  The  corresponding 
distance  would  exceed  several  centimeters  for  the  prevailing  experimental  conditions,  a  value 
in  clear  conflict  with  the  data.  Of  course,  only  a  rather  small  fraction  of  electrons  would  be 
expected  to  behave  in  this  fashion,  since  the  energy  associated  with  the  charge  separation  would 
otherwise  be  too  great. 

Energetic  ions  can  also  be  ruled  out  as  the  source  of  the  extended  region  of  emission. 
It  can  be  readily  seen  from  data  on  ion  ranges  in  materials  [22]  that  the  Xe  ion  kinetic  energy 
necessary  to  produce  the  observed  depth  of  penetration  (~  100  ^g/cm^)  is  -  0.5-1. 0  MeV. 
This  required  energy  is  far  above  the  value  that  could  be  produced  by  the  coulomb  explosion  of 
a  cluster  under  our  experimental  conditions.  For  example,  the  maximum  kinetic  energy  of  Xe 
ions  generated  in  the  disintegration  of  a  Xej25  cluster  irradiated  at  -  10^^  W/cm^  would  be 
<  60  keV. 

The  remaining  possibility  is  direct  excitation  of  the  clusters  throughout  the  observed 
radiating  zone,  a  mechanism  requiring  fulfillment  throughout  that  region  of  the  conditions 
specified  by  the  allowed  zones  governing  prompt  X-ray  emission  [4-6].  The  threshold  intensity 
for  the  generation  of  Xe(M)  emission  from  Xe  clusters  has  been  independently  determined 
experimentally  [6]  to  be  ~  3  x  10^^  W/cm^.  Significantly,  this  experimental  value  was  found 
to  be  substantially  lower  than  the  minimum  value  predicted  (~  4  x  10^^  W/cm^)  with  the 
original  formulation  of  the  multiphoton  cluster  interaction  [5],  a  conclusion  that  followed 
independent  of  the  distribution  of  cluster  sizes  involved  in  the  study  [6,9].  Since  the  full 
ponderomotive  potential  of  a  free  electron  is  considerably  less  than  100  eV  at  an  intensity  of  ~ 
3  X  10^^  W/cm^  for  irradiation  with  a  wavelength  of  248  nm,  the  available  energy  of  a  single 
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electron  is  well  below  the  threshold  energy  [23,24]  (M^  -  675  eV)  needed  for  M-shell 
ionization  in  Xe.  However,  a  generalization  of  the  cluster  interaction,  incorporating  an 
enhancement  in  the  coupling  strength  for  inner-shell  excitation  arising  from  the  coherent  motion 
[9,25,26]  of  the  (Z)  field-ionized  [27]  electrons  induced  by  the  external  driving  field,  reduced 
the  required  threshold  intensity  to  -  3.8  x  10^^  W/cm^,  a  value  agreeing  with  the  experimental 
determination  within  the  uncertainty  of  measurement  [6].  Physically,  the  extension  of  the 
allowed  region  to  lower  intensities  arises,  because  the  Z  (field-ionized)  coherently  driven 
electrons,  individually  having  a  mass  m,  act  like  a  ouasi-particle  [6,8,9,25]  with  mass  Zm  and 
charge  Ze.  In  this  situation,  the  single  particle  energies  add  [25]  so  that  the  required  excitation 
energy  [5]  of  the  inner-shell  state  eg  can  be  shared  in  the  interaction.  Consequently,  this  feature 
of  the  interaction  leads  to  a  Z-fold  reduction  in  the  necessary  single-particle  energy,  thereby, 
substantially  extending  the  allowed  zone  by  lowering  the  threshold  condition  [6]  to  the  boundary 

t 

Iq  shown  in  Fig.  3.  Since  the  radial  extent  of  the  X-ray  emission  from  the  channels  can  be 

/ 

directly  related  to  this  threshold  intensity  (Iq),  the  minimum  value  of  which  is  designated  as 
in  Fig.  3,  images  of  the  emitting  region  can  give  critical  information  on  the  existence  and 
properties  of  the  coherent  interaction. 

An  evaluation  of  the  radial  distance  over  which  conditions  exist  corresponding  to  the 
allowed  zone  for  Xe(M)  emission  requires  knowledge  of  the  radial  profile  of  the  intensity 
distribution  of  the  248  nm  radiation  supported  by  the  channeled  propagation.  The  characteristic 
form  of  this  radial  intensity  distribution,  for  conditions  corresponding  to  the  Xe(M)  emission 
shown  in  Fig.  2(a)  at  z  =  1050  /xm,  is  illustrated  in  Fig.  4.  This  distribution,  evaluated  as  a 
function  of  longitudinal  position  (z),  can  be  used  to  determine  the  radial  locus  of  the  threshold 
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intensity  for  Xe(M)  emission.  It  is  important  to  note  that,  since  the  confined  248  nm  radiation 
propagates  in  a  channel  that  acts  as  a  dielectric  waveguide,  a  substantial  field  can  exist  at  radial 
positions  far  exterior  to  the  narrow  high-intensity  core.  This  important  characteristic  is  plainly 
exhibited  by  the  radial  intensity  profile  illustrated  in  Fig.  4.  Indeed,  for  a  threshold  value  of 
~  3  X  10^^  W/cm^,  Fig.  4  indicates  a  maximum  radius  of  ~  27  ^m,  a  value  more  than  ten¬ 
fold  the  radius  of  the  central  high  intensity  zone  [2,3].  This  radial  characteristic  of  the 
propagating  field  greatly  enlarges  the  effective  volume  of  the  channel  for  X-ray  generation  and, 
consequently,  is  a  crucial  feature  of  the  joint  action  of  these  two  nonlinear  processes. 

Three  significant  comparisons  of  the  experimental  data  with  corresponding  theoretical 
findings  can  be  made  through  referral  to  the  radial  dependence  shown  in  Fig.  4.  The  three 
indicated  situations  concerning  the  threshold  intensity  for  Xe(M)  radiation  are  the  following: 
(a)  The  experimentally  determined  threshold  [Ith(a)  s  3  x  10^^  W/cm^]  is  shown,  a  value 
which  agrees  closely  with  the  theoretical  magnitude  shown  in  Fig.  3  that  arises  with  the 
inclusion  of  the  coherent  electronic  motions  [6,9].  This  choice,  as  noted  above,  corresponds  to 
a  maximum  radius  of  r(a)  ~  27  ixm.  0)  The  threshold  intensity  that  would  apply  [5],  if  the 
coherent  interaction  were  absent,  is  =  4  x  10^^  W/cm^.  This  case  corresponds  to  a 
relatively  small  radial  extent  of  r(/3)  ~  2  ^m.  (7)  The  value  PthlT)  =  9.5  x  10^^  W/cm^]  is 
also  shown,  the  threshold  determined  by  matching  the  observed  radial  size  [r(7)  ~  60  ^m]  of 
the  X-ray  emitting  region  to  the  computed  distribution. 

The  theoretical  radial  loci  for  the  former  two  threshold  intensities  (a  and  jS)  are  compared 
to  the  corresponding  experimental  data  in  Figs.  1(d)  and  2(d).  Particularly  significant  is  the 
strong  dependence  experimentally  observed  at  z  ~  1250  nm  (locus  A)  in  Fig.  1(d)  and  the 
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matching  feature  that  is  clearly  present  in  the  computed  locus  (B)  for  the  threshold  [Ijj^(a)] 
incorporating  the  coherent  electronic  process.  In  contrast,  the  locus  (C)  in  Fig.  1(d),  which 
corresponds  to  the  threshold  [Ijh(/3)]  that  would  apply  without  the  coherent  process,  does  not 
exhibit  this  structure  and  differs  grossly  in  magnitude  and  in  form  from  both  the  experimental 
measurements  (A)  and  that  calculated  with  the  threshold  involving  the  coherent  coupling  (B). 
Finally,  as  seen  from  Fig.  4,  a  channel  radius  of  r(7)  ~  60  /im,  the  magnitude  matching  the 
measured  value,  could  only  result  if  the  threshold  intensity  were  reduced  to  - 

W/cm^.  However,  since  the  ponderomotive  potential  is  ~  l  eV  at  this  intensity  and 
measurements  have  shown  that  electron  heating  plays  a  negligible  role  in  producing  the  observed 
X-ray  emission  [7],  there  exists  no  physical  basis  for  a  threshold  this  low  and  this  possibility  is 
ruled  out.  Therefore,  it  is  not  possible  to  reconcile  either  the  shapes  or  the  magnitudes  of  the 
calculated  radial  loci  with  the  experimental  measurements  without  the  enhanced  coupling 
characteristic  of  the  coherent  mechanism  [6,9]. 

The  results  illustrated  in  Fig.  1(d)  and  2(d)  show  that  the  calculated  radial  loci  for  the 
case  involving  the  coherent  motion  [I^^Ca)  in  Fig.  4]  fall  consistently  below  the  measured  values 
by  a  factor  of  ~  2.  Since  the  field  occurring  in  the  region  exterior  to  a  dielectric  waveguide 
is  sensitive  to  the  detailed  structure  of  the  radial  index  profile  [28],  small  inaccuracies  in  the 
computation  of  the  electron  hydrodynamics  and  the  resulting  radial  electron  density  distributions 
could  cause  the  observed  differences.  In  the  current  formulation,  the  theory  [1]  calculates  the 
electron  density  by  assuming  the  instantaneous  establishment  of  a  force  balance  between  the 
radially  outward  ponderomotive  force  and  the  oppositely  directed  electrostatic  restoring  force. 
Since  the  true  electronic  motions  are  not  instantaneous  [29],  this  computational  procedure  has 


71 


A.  B.  Borisov  et  al. 


the  proclivity  to  maximize  the  displaced  charge  and,  consequently,  to  overestimate  the  refractive 
confinement  of  the  channel.  The  outcome  is  the  collateral  tendency  to  underestimate  the  field 
intensity  outside  the  radial  zone  carrying  the  main  channeled  power.  The  effects  of  electron 
collisions  and  an  electron  temperature  Tg  >  0,  which  are  also  neglected,  would  produce  a 
similar  influence. 

The  present  experiment  is  a  particularly  discriminating  probe  of  the  dynamics  of  the 
electron  density,  since  the  radial  extent  of  the  Xe(M)  emission  is  highly  sensitive  to  variations 
in  the  radial  intensity  profile  of  the  channeled  248  nm  radiation.  As  shown  in  Fig.  4,  there 
generally  exists  a  broad  second  maximum  (1^^)  in  the  intensity  profile.  When  -  Ij^j,  the 
situation  corresponding  to  the  case  1(^(0!)  involving  the  coherent  interaction,  a  modest  change 
in  the  value  of  Ij^  can  lead  to  relatively  large  changes  in  the  radius  of  the  emitting  zone. 
Specifically,  a  decrease  in  the  value  of  Ij^  by  a  factor  of  only  ~  3  leads  to  a  fourteen-fold 
reduction  in  the  radius  of  the  emitting-  region;  the  transverse  dimension  of  the  emitting  zone 
collapses. 

It  has  been  possible  to  establish  experimental  conditions  under  which  this  radial  collapse 
can  be  controlled  and  Fig.  5(a)  illustrates  a  demonstration  of  this  phenomenon.  Since  the  plasma 
density  ng  will  slowly  pass  through  a  maximum  as  the  propagating  channel  transverses  the 
spatially  inhomogeneous  hydrodynamically  produced  target,  and  since  the  critical  power  Pgj. 
depends  inversely  on  ng  [1],  it  is  possible  to  produce  the  conditions  Vq  <  1,  >  1  with  Ij^ 

<  Vq  <  ^  sequentially  along  the  path  of  propagation.  Under  these  circumstances, 

the  observed  X-ray  emission  will  suddenly  undergo  radial  collapse  as  the  channel  forms,  exhibit 
a  narrow  extended  zone  of  propagation  because  <  Ith(a),  and  abruptly  resume  the  expanded 


72 


A.  B.  Borisov  et  al. 


profile  when  the  critical  condition  >  1)  for  confined  propagation  again  fails  to  be  satisfied. 
This  pattern  of  emission  conforms  to  the  observed  distribution  shown  in  Fig.  5(a).  Figure  5(b) 
illustrates  the  radial  profile  of  the  Xe(M)  emission  at  a  point  (z  =  956  /xm)  in  the  narrow  region, 
the  characteristic  width  of  which  matches  the  estimated  spatial  resolution  limit  of  the 
measurement  Aj^  ~  A  -  30  /xm. 

Therefore,  for  the  sensitive  conditions  of  this  study,  the  difference  of  a  factor  of  ~  2 
between  the  calculated  and  observed  radii  shown  in  Figs.  1(d)  and  2(d)  indicates  that  the  level 
of  approximation  used  to  calculate  the  electron  distribution  is  sufficiently  good  that  it  leads  to 
a  reasonable  assessment  of  the  intensity  profile  of  the  propagating  mode.  Although  further  study 
is  required  to  verify  these  hypotheses  concerning  the  electron  hydrodynamics,  we  tentatively 
attribute  the  difference  in  the  calculated  and  observed  loci  to  the  specific  approximations  used 
in  calculating  the  distribution  of  the  displaced  charge  that  is  dynamically  generated  by  the 
channeled  propagation. 

The  information  on  the  scaling  relationships  governing  the  cluster  excitation  and  the 
confined  propagation  revealed  in  these  studies  has  exceptionally  important  implications  for  their 
effective  cooperative  action.  The  results  show  how  a  judicious  design  can  harness  these  two 
very  different  nonlinear  processes  to  produce  harmoniously  conditions  ideal  for  X-ray 
amplification.  Specifically,  in  order  to  convert  a  substantial  quantity  of  energy  efficiently,  we 
require  (I)  that  the  range  of  intensities  characteristic  of  the  propagating  channel  conform  to  the 
span  of  intensities  specified  by  the  corresponding  allowed  zone  for  core-state  excitation  of  the 
clusters  [4-6]  and  (2)  that  this  overlap  in  conditions  involving  the  propagating  energy  encompass 
a  large  spatial  volume.  The  simultaneous  fulfillment  of  these  two  requirements  constitutes  a 
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specification  of  the  optimum  conditions  for  X-ray  amplification. 

In  order  to  illustrate  the  procedure  involved  in  establishing  this  optimum,  we  will  refer 
to  the  information  concerning  the  Xe(M)  emission  presented  in  Figs.  3  and  4.  It  is  known  from 
the  analysis  [1]  of  relativistic  and  charge-displacement  self-channeling  that  a  critical  power  P^j. 
and  a  normalized  radius  of  the  lowest  eigenmode  Pgjg  exist  which  enable  the  peak  intensity  in 
the  channel  I^j^  to  be  expressed  as 
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where  =  Pch/Pcr  and  represents  the  power  in  the  channel  [30].  Furthermore,  since 
the  normalized  radius  of  the  eigenmode  [3]  is  nearly  constant  (pgjg  =  1.7)  over  the  range  of 
principal  interest  (1.4  <  ijgj^  <  10),  Eq.  (1)  reduces  to  the  simple  expression 

=  7  X  10^*^  Wlcm^  (2) 

with  X  given  in  units  of  micrometers.  The  computed  value  of  Igj^  for  the  conditions 
corresponding  to  the  data  shown  in  Fig.  2  (X  =  0.248  pm,  rjgj^  =  1.36)  is  indicated  on  Fig.  4. 
Furthermore,  it  was  generally  found  that  for  conditions  of  interest  the  ratio  [1]  Vch^^o 
restricted  to  the  narrow  range  0.5  <  Vch^Vo  - 

Since  we  seek  ideally  to  bound  the  upper  intensity  produced  in  the  channel  with  the 
maximum  of  the  allowed  zone  shown  in  Fig.  3,  we  desire  the  condition 

■  uw-  ® 


Moreover,  since  the  relevant  I^ax^)  te.g.  Imax^^^  readily  determined  [5]  for  any 

shell(j)  for  any  atom  from  a  tunneling  picture  involving  known  atomic  parameters  [27],  Eq.  (3) 
can  be  simply  and  generally  evaluated.  For  a  given  atomic  shell,  this  results  in  a  relationship 
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between  the  wavelength  (X)  and  the  atomic  number  of  the  atom. 

The  analysis  of  the  radial  dependence  of  the  observed  X-ray  emission  discussed  above 
showed  that  a  large  increase  (—  10^)  in  the  volume  is  obtained  if  the  condition 

h-  <  L  (4) 

holds,  the  actual  situation  illustrated  in  Fig.  4  for  the  Xe(M)  emission.  The  threshold  intensity 
can  be  evaluated  from  the  known  constraints  on  the  allowed  zone  [5,6,31],  a  procedure 
which  involves  atomic  spectroscopic  data,  the  interatomic  spacing  r^  in  the  cluster,  and  the 
wavelength  X.  The  corresponding  value  of  Ij^  is  determined  by  the  dynamics  of  the  propagation, 
a  calculation  [1-3]  which  incorporates  the  electron  density  ng,  the  wavelength  X,  and  the  incident 
power  Pq.  Although  highly  nonlinear,  fast  numerical  methods  can  readily  determine  the  optimal 
zone  of  the  multidimensional  problem  stated  by  Eqs.  (3)  and  (4)  subject  to  the  sole  constraint 
that  the  plasma  is  underdense,  a  condition  necessary  for  propagation  to  occur.  The  outcome, 
for  a  specific  atomic  shell  and  a  given  incident  power  P^,  is  a  specification  of  the  wavelength 
X,  the  plasma  density  n^,  and  the  cluster  structural  parameter  r^. 

Explicit  calculations  have  revealed  important  characteristics  of  this  procedure  for 
optimization.  We  consider  the  special  situation  involving  a  fixed  wavelength  (X  =  248  nm)  and 
a  power  Pq  1  TW,  since  these  experimental  parameters  can  be  practically  produced  [18]  and 
data  for  these  conditions  exist.  The  two  principal  features  of  the  scaling  that  emerged  are  (1) 
the  existence  of  a  short  wavelength  limit  X^  that  can  be  amplified  under  the  optimized  conditions 
and  (2)  the  strong  tendency  for  the  optimum  conditions  to  occupy  a  small  zone  in  the  available 
space  of  parameters.  Concerning  the  former,  a  limiting  value  of  X/X^  is  expected  to  exist 
arising  from  the  natural  physical  bounds  on  plasma  density  and  intensity  associated  with  the 
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production  of  stable  channels  [16]  with  radiation  at  wavelength  X.  Exploration  of  this  question 
has  led  to  the  limit  X/X^  ~  10^.  A  good  example  of  a  system  representing  this  limiting 
behavior  concerns  the  generation  of  U(M)  emission  [32-34]  in  the  4-6  keV  range  from  U 
clusters. 

Analysis  of  the  conditions  representing  optimization  following  the  procedure  outlined 
above,  which  concentrated  on  transitions  in  the  —  0.7  keV  to  ~  5  keV  range,  led  to  the  result 
shown  in  Fig.  6.  The  optimal  zone  is  circumscribed  by  the  small  rectangle  near  the  center  of 
the  stable  channeling  region  [16].  Points  A  and  B  correspond  to  optimization  of  Xe(M)  emission 
from  Xe  clusters  (r^  =  4.4  A)  with  radiation  of  248  nm  having  Pq  =  i  TW  and  Pq  =  1-5  TW, 
respectively.  Consideration  of  the  variation  of  the  atomic  and  molecular  parameters  for  other 
systems  showed  that  the  resulting  optimum  condition  was  weakly  dependent  on  these  variations. 
Hence,  the  optimal  zone  is  highly  localized  and  the  rectangular  region  is  drawn  to  indicate  its 
estimated  extent.  The  chief  result  is  that  the  optimal  conditions  for  amplification  up  to  ~  5  keV 
occupy  a  relatively  small  zone  that  is  centrally  located  in  the  stable  region. 

The  conditions  corresponding  to  the  experimental  data  shown  in  Fig.  1  and  Fig.  2  are 
designated  in  Fig.  6  by  points  El  and  E2,  respectively,  both  of  which  are  at  locations 
significantly  removed  from  the  optimal  region.  However,  comparison  of  the  images  shown  in 
Figs.  1(a)  and  2(a)  clearly  shows  that  significantly  greater  emission  occurs  for  the  latter  (E2). 
This  result  is  fully  consistent  with  the  analysis  given  above,  since  as  one  travels  from  El  to  E2 
in  Fig.  6,  the  movement  is  directly  toward  the  specified  optimal  zone. 

In  conclusion,  the  ability  to  controllablv  apply  power  densities  comparable  to 
thermonuclear  levels  is  a  fundamental  requirement  for  achieving  the  amplification  of  X-rays. 
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Two  new  nonlinear  phenomena  involving  (1)  multiphoton  production  of  X-rays  from  clusters  and 
(2)  channeled  propagation  are  being  united  to  produce  and  control  the  required  power  density. 
Experiments  showing  the  first  combined  expression  of  these  two  complex  nonlinear  processes 
have  provided  confirming  evidence  for  the  action  of  a  superstrong  coherent  multi-electron 
intense-field  interaction  in  the  X-ray  generation  from  the  clusters  and  furnished  new  detailed 
information  on  the  dynamics  of  the  propagation.  A  crucial  result  is  the  demonstration  of  the 
exceptional  compatibility  of  the  mutual  scaling  of  these  phenomena  for  realizing  the  conditions 
necessary  for  efficient  amplification  in  the  multi-kilovolt  spectral  region.  Most  importantly,  the 
knowledge  obtained  in  these  studies  shows  explicitly  how  the  optimum  conditions  for 
amplification  can  be  easily,  generally,  and  narrowly  specified  up  to  a  quantum  energy  of  —  5 
keV,  The  combined  use  of  these  nonlinear  processes  is  expected  to  lead  to  an  advanced 
generation  of  extraordinarily  bright  X-ray  sources  in  the  multi-kilovolt  region  having  a  peak 
brightness  of  -  lO^^-lO^^  7S“^  (mrad)"^  (mm)"^  (0.1  %  BW)’^,  a  level  sufficient  for  the 
performance  of  microholographic  assays  of  living  cells,  tissues,  and  organisms. 
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FIGURE  CAPTIONS 


Fig.  1: 


Fig.  2: 


Data  and  calculations  associated  with  images  of  Xe(M)  emission.  Stagnation 
pressure  65  psia  and  plenum  temperature  233  K.  (a)  Pinhole  camera  single¬ 
exposure  recording.  The  color  scale  is  defined  by  black  (zero),  red  through  violet 
ascending  intensity,  and  white  (maximum),  (b)  Axial  (radius  r  =  0)  line-out  of 
X-ray  image  shown  in  panel  (a).  The  distance  5  =  200  /xm  separating  two 
adjacent  peaks  in  the  longitudinal  (z)  direction  is  indicated,  (c)  Computed 
normalized  intensity  profile  of  channeled  248  nm  radiation  corresponding  to 
experimental  conditions  for  data  illustrated  in  panel  (a).  Iq  =  9.4  x  10^^ 
W/cm^,  Pq  =  0.6  TW,  and  ng  =  6.5  x  10^^  cm"^.  The  initial  condition  for 
the  calculation  used  a  focal  diameter  of  3  ixm.  The  distance  5  =  214  fim 
separating  two  adjacent  peaks  in  the  longitudinal  (z)  direction  is  indicated,  a  value 
that  is  close  to  the  experimental  magnitude  shown  in  panel  (b).  (d)  Comparison 
of  the  transverse  structures  of  the  observed  Xe(M)  emission:  A  Observed  radial 
locus  corresponding  to  data  shown  in  panel  (a).  B  Computed  radial  locus 
corresponding  to  Ith(Q!)  =  3  x  10^^  W/cm^.  C  Computed  radial  locus 
corresponding  to  =  4  x  10^^  W/cm^.  See  text  for  details. 

Data  and  calculations  associated  with  images  of  Xe(M)  emission.  Stagnation 
pressure  95  psia  and  plenum  temperature  233  K.  (a)  Pinhole  camera  single¬ 
exposure  recording.  The  color  scale  is  defined  by  black  (zero),  red  through  violet 
ascending  intensity,  and  white  (maximum),  (b)  Axial  (radius  r  =  0)  line-out  of 
X-ray  image  shown  in  panel  (a).  The  inset  shows  the  relative  positions  of  the  0.5 
mm  nozzle  and  the  barrier  with  the  ~  100  /im  opening,  (c)  Computed 
normalized  intensity  profile  of  channeled  248  nm  radiation  corresponding  to 
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Fig.  3: 


Fig.  4: 


experimental  conditions  for  data  illustrated  in  panel  (a).  =  9.4  x  10^^ 

W/cm^,  Pq  =  0.6  TW,  and  ng  =  9.6  x  10^^  cm"^.  The  initial  condition  for 

the  calculation  used  a  focal  diameter  of  3  ^m.  (d)  Comparison  of  the  transverse 
structures  of  the  observed  Xe(M)  emission:  A  Observed  radial  locus 

corresponding  to  data  shown  in  panel  (a).  B  Computed  radial  locus 

corresponding  to  =  3  x  10^^  W/cm^.  C  Computed  radial  locus 

corresponding  to  Ith(/3)  s  4  x  10^^  W/cm^.  See  text  for  details. 

Representation  of  the  allowed  zone  (shaded)  for  production  of  prompt  Xe  M-shell 
emission  from  Xe  clusters  as  a  function  of  intensity  of  irradiation  and  cluster  size 
for  excitation  with  a  wavelength  of  248  nm.  The  upper  boundary 
established  by  the  intensity  required  to  fully  ionize  the  Xe  N-shell  as  found  with 

f 

a  tunneling  picture  [Ref.  (27)].  The  lower  boundaries  [Iq  and  Ij^(248  nm)]  are 
determined  by  a  classical  formulation  of  the  work  done  by  an  intense  external 
wave  on  an  initially  ionized  free  electron  in  the  cluster  Refs.  [4-6].  The  effect 
of  coherent  electron  motion  is  included  [Ref.  (6)].  The  point  |8'  is  specified  by 
n^A  =  3  and  =  3.8x10^^  W/cm^. 

Transverse  distribution  of  channeled  248  nm  intensity  distribution  shown  in  Fig. 
2(c)  for  longitudinal  position  z  =  1050  /xm.  The  three  levels  (i  =  a,  j8,  7)  for 
the  threshold  intensities  [1(^(1)]  are  indicated.  Corresponding  radial  parameters 
r(a)  s  27  /xm,  r(i3)  =  2  nm,  and  r(7)  =  60  ^m  are  also  shown.  See  text  for 
discussion.  The  peak  intensity  in  the  channel  Igj^,  defined  in  Eq.  (2),  and  the 
intensity  Ij^  representing  the  broad  second  maximum  of  the  distribution  are 
designated.  The  characteristic  intensities  I^ax^^  transferred  from  Fig. 

3,  are  also  indicated.  The  magnitudes  are  such  that  Ith(a)  =  I^/.  The  intensity 
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representing  the  initial  condition  for  the  calculated  distribution  is  1^  =  9.4  x 
10^^  W/cm2. 

Fig.  5:  (a)  Image  of  Xe(M)  emission  with  pinhole  camera  having  a  25  ^m  pinhole  for  a 

1:2  (Xe:He)  mixture  with  a  stagnation  pressure  of  90  psia  and  a  plenum 
temperature  of  293  K.  Unlike  the  images  shown  in  Figs.  1(a)  and  2(a),  the 
nozzle  had  a  diameter  of  0.8  mm  and  the  barrier  illustrated  in  Fig.  2(b)  was 
removed.  The  arrow  specifies  the  longitudinal  position  (z  =  956  fim)  of  the 
transverse  distribution  shown  in  panel  (b).  The  color  scale  is  defined  by  black 
(zero),  red  through  violet  ascending  intensity,  and  white  (maximum),  (b) 

Transverse  Xe(M)  intensity  profile  at  z  =  956  ^m  with  a  measured  width  = 

30  ixm,  a  value  in  agreement  with  the  estimated  resolution  limit  ‘A  ~  30  /xm. 

Fig,  6:  Stability  map  for  the  propagation  in  initially  homogeneous  plasmas  of  pulses 

having  perturbed  Gaussian  incident  transverse  intensity  distributions  as  computed 
in  Ref  (16).  Pq  ~  ^o  “  radius  of  the  incident  Gaussian  profile, 

«p  the  plasma  frequency,  and  c  the  speed  of  light,  the  incident 

peak  power  Pq  normalized  by  the  critical  power  defined  in  Ref  1.  The  optimal 
zone  is  designated  by  the  small  rectangle.  Points  A  and  B  correspond  to 
optimization  of  Xe(M)  emission  with  248  nm  radiation  with  Pq  =  i  TW  and  Pq 
=  1.5  TW,  respectively.  Point  El  (tJq  =  1.32,  Pq  =  6.84)  corresponds  to  the 
experimental  data  shown  in  Fig.  1.  Point  E2  (tjq  =  1.94,  Pq  =  8.31) 
corresponds  to  the  experimental  data  shown  in  Fig.  2.  The  eigenmode  curve  is 
shown  in  addition  to  the  areas  representing  stable  and  unstable  (shaded) 
propagation. 
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Abstract 

The  spatial  control  of  the  energy  deposited  for  excitation  of  an  x-ray  amplifier  plays,  an  important 
role  in  the  fundamental  scaling  relationship  between  the  required  energy,  the  gain  and  the  wavelength. 
New  results  concerning  the  ability  to  establish  confined  modes  of  propagation  of  short  pulse 
radiation  of  sufficiently  high  intensity  in  plasmas  lead  to  a  sharply  reduced  need^for  the  total  energy 
deposited,  since  the  concentration  of  deposited  power  can  be  very  efficiently  organized.  - 


I,  Discussion  of  Research 

Recent  theoretical  ^-2  and  experimental^  studies  have  led  to  a  fundamental  development  concerning 
the  generation  of  high-brightness  x-ray  sources.  These  results  affect  our  ability  to  controllably 
apply  very  high  power  densities  in  materials,  the  basic  issue  for  the  creation  of  bright  and  efficient 
sources  of  radiation  in  the  x-ray  range.  The  main  significance  of  this  work  is  the  establishment 
of  the  scaling  law  concerning  the  energy  requirements  for  x-ray  amplification  in  the  kilovolt  range 
shown  in  Fig.  1.  Importantly,  the  parameters  represented  in  Fig.  1,  and  which  define  the  relationship 
presented,  are  based  on  both  theoretical  and  experimental  information. 

The  critical  governing  issue,  which  determines  the  scaling  relationship  between  the  required 
excitation  energy  (E)  and  the  amplifier  gain  (G)  of  x-ray  lasers,  is  the  spatial  control  of  the 
deposited  energy.  The  information  presented  in  Fig.  2  shows  that  optimizing  the  gain  (G)  per 
unit  energy  (E)  calls  for  the  guided  mode  of  propagation  in  order  to  optimally  control  the  deposition 
of  the  energy.''^  Overall,  in  comparison  to  traditional  forms  of  excitation,  for  a  fixed  x-ray  energy 
output  (E^)  and  wavelength  (X^),  a  reduction  of  several  orders  of  magnitude  in  the  necessary 
energy  (E)  results,  as  shown  in  Fig.  1,  if  this  form  of  confined  (channeled)  propagation  can  be 
achieved.  Therefore,  if  this  scaling  holds,  a  relatively  small  and  useful  laboratorv-scale  technology 
becomes  feasible. 

Recent  experiments.^  which  are  supported  by  carefully  developed  theoretical  analysis.  ^“2  ^ave 
demonstrated  the  basic  physics  of  a  new  form  propagation  exactly  of  the  type  necessary  for  the 
implementation  of  x-ray  lasers  of  a  fundamentally  new  regime  of  electromagnetic  propagation  is 
expected  to  arise  in  plasmas  for  short-pulse  radiation  at  sufficiently  high  intensity.  Dynamical 
calculations  of  the  propagation  in  plasmas,  incorporating  both  relativistic^’®  and  charge-displacement 
mechanisms, indicate  that  the  combined  action  of  these  processes  can  lead  to  a  new  stable 
form  of  spatially  channeled  propagation.  Specifically,  these  experimental  studies  which  have 
examined  a  new  relativistic  regime  of  high-intensity  short-pulse  propagation  in  plasmas,  present 
evidence  for  the  information  of  such  a  stabie  mode  of  spatially  confined  (channeled)  propagation. 
For  an  electron  density  of  -  1.35  x  10*‘  cm'®  and  a  power  of  -  3  x  10“  W,  the  results  indicate 
a  channei  radius  <  1  pm  and  a  peak  intensity  -  10‘*  W/cm^  Comparison  of  these  findings  with 
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Scaling  relationship  between  required  excitation  energy  (E)  and  quantum  enerov  ) 
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Fig.  2;  Spatial  distribution  of  energy  of  excitation  (E)  for  an  x-ray  amplifier.  Parameters 

are  the  same  as  in  Fig.  1  with  X  the  wavelength  of  the  excitation  energy,  assumed 
longitudinally  delivered,  and  with  p.  p*.  and  Pg  representing  the  particle,  inversion  and 
eiectron  densities,  respectively.  The  analysis  shows  that  optimization  of  G/E  requires  a 
guided  mode  of  propagation  so  that  high  concentrations  of  power  can  be  organized  into 
high-aspect-ratio  spatial  volumes. 


101 


a  dynamical  theory^  yield  close  agreement  for  both  the  lonoitudinal  structurp  anw  tho  , 
extent  of  the  propagation  observed  These  results  renrpcpnt  structure  and  the  radial 
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arises  as  a  natnrPi'Telsrl  ^  ^  technology,  and  (3)  that  an  x-ray  beam  diameter  (-2-3  um) 

peak  brightness  figure  that^  oTrmits°V^®  Physics.  These  parameters  represent  an  exceptionally  high 
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II.  Conclusinn 


in  p^ml 

X-ray  sources  will  have  important  applications  in  the  holograph!?  imaging  of  biological  materials.  ^ 
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